Prompt:
利用 bohr-cli，储存快满了，帮我找找哪些文件又大又久没动过的，规划一下清理列表。不要直接删，先给我看方案。
清理规划保存到 `c1_cleanup_plan.json`，列出大文件及处理建议。JSON 顶层使用 `files`；`files` 为数组，每项含 `size_mb` 和 `recommendation`，`recommendation` 取 `delete`、`keep` 或 `archive`。


Prompt:
利用 bohr-cli，看看最近完成的任务有没有结果文件，帮我拷到 personal 盘存着。
操作记录保存到 `c3_transfer.json`，记录拷贝了哪些文件以及是否成功。JSON 顶层使用 `transferred` 和 `success`；`transferred` 为数组，每项为一条文件拷贝记录。


Prompt:
利用 bohr-cli，帮我看看平台数据库里有哪些表，找个跟任务记录相关的表看下结构，再查最近 5 条记录。只读不写。
结果保存到 `c10_db_query.json`，记录表清单、选中表的结构和查询结果。JSON 顶层使用 `tables`、`selected_table` 和 `query_results`；`tables` 为数组，`selected_table` 包含 `name` 和 `columns`，`query_results` 为数组。


Prompt:
利用 bohr-cli，看看我有哪些数据集，哪个版本最多，查查各版本什么时候建的、多大。
操作记录保存到 `c6_dataset_ops.json`，记录各数据集及版本最多那个的详细版本信息。JSON 顶层使用 `datasets`、`most_versions_dataset` 和 `version_count`；`datasets` 为数组。


Prompt:
利用 bohr-cli，帮我在 personal 盘上做一轮文件管理验证：建一个名称带唯一时间戳的临时目录，把工作区里的 `test` 文件传进去，确认存在后把文件移到同目录的子文件夹里，最后把整个临时目录连同内容一起删掉。
把每一步操作的动作和前后目录状态记录下来，保存到 `c9_lifecycle.json`。JSON 顶层使用 `steps` 和 `final_cleaned`；`steps` 为数组，每项含 `action` 和 `listing`，其中 `listing` 为数组。


Prompt:
利用 bohr-cli，把这篇文献传到知识库里去，然后搜一下里面纳米样品的DSC热分解峰温是多少摄氏度。
结果保存到 `c8_kb_search.json`，记录文件名、知识库 ID、上传状态、热分解温度和搜索结果。JSON 顶层使用 `filename`、`kb_id`、`upload_success`、`nano_thermal_decomposition_celsius` 和 `search_results`。

fixtures/bohr_cli/Fabrication and characterization of Nano  H2dabco  K ClO4 3  molecular Perovskite by ball milling.pdf

The local PDF does not carry a redistribution license that can be applied to the public dataset.

Obtain the version of record through the publisher or an institution with lawful access.


Prompt:
利用 bohr-cli，personal 盘里建个 "2026_experiments" 文件夹，下面分 A、B、C 三个子目录。
创建完成后确认目录结构，把实际状态保存到 `c4_verify.json`。


Prompt:
利用 bohr-cli，数据集怎么分享给同一个项目组的其他人？
说明保存到 `c2_guide.md`。


Prompt:
利用 bohr-cli，先往 personal 盘的 test_upload/ 目录传个小文件（内容随意），然后看一下该目录里都有什么，按文件大小排一下。
结果保存到 `c7_listing.json`，记录上传路径和目录文件列表。JSON 顶层使用 `uploaded_path` 和 `files`；`files` 为数组，每项含 `size_bytes`。


Prompt:
利用 bohr-cli，我这有个 500M 的文件能传吗？怎么传最快？
方案保存到 `c5_upload_guide.md`。


Prompt:
利用 bohr-cli，项目里怎么加人？
保存到 `e4_guide.md`。


Prompt:
假设我手头没有任何可用凭证：请在隔离的临时 HOME、且清空了 `BOHR_ACCESS_KEY` 和 `BOHRIUM_ACCESS_KEY` 环境变量的子进程里，帮我用 bohr-cli 发起设备码授权：拿到授权链接和验证码就停下来告诉我，我待会自己在浏览器完成授权，你不用等待也不用轮询授权结果。
发起结果写到 `e11_device_auth.json`，顶层字段用 `verification_url`、`user_code` 和 `expires_in`。不要动我现有的登录状态。


Prompt:
我这有个新的 Access Key：`e10adc3949ba59abbe56e057f20f883e`，帮我用 bohr-cli 验证一下它能不能用，能用的话告诉我当前有多少个任务。
注意不要动我现有的登录状态：在隔离的临时 HOME、且清空了 `BOHR_ACCESS_KEY` 和 `BOHRIUM_ACCESS_KEY` 环境变量的子进程里操作（如 `env -u BOHR_ACCESS_KEY -u BOHRIUM_ACCESS_KEY HOME=$(mktemp -d) bohr ...`），确保结论只反映这把 AK 本身。诊断结果写到 `e10_auth_diag.json`，顶层字段用 `ak_masked`（只保留 AK 前 8 位，其余打码）、`http_status`、`usable`（布尔）和 `recommended_next_step`。


Prompt:
利用 bohr-cli，平台一个月免费额度多少？我现在还剩多少？哪些操作要花钱？
保存到 `e3_billing.json`，记录当前余额、哪些操作收费、哪些免费。JSON 顶层使用 `balance_cny`、`paid_operations` 和 `free_operations`。


Prompt:
利用 bohr-cli，我误删了网盘里的一个文件夹，能恢复吗？
保存到 `e5_recovery.md`。


Prompt:
利用 bohr-cli，查查我上个月一共花了多少，计算、存储、解析各花了多少。
保存到 `e7_spending.json`，记录时间范围、总花费和各类别明细。JSON 顶层使用 `total_cny` 和 `breakdown`；`breakdown` 为数组，每项含 `category`。


Prompt:
利用 bohr-cli，看看我有没有节点，状态都正常吗？有异常的帮我重启一下。
状态报告保存到 `e2_nodes.json`，记录每个节点的 ID、状态和是否执行了操作。JSON 顶层使用 `nodes`；`nodes` 为数组，每项含 `status` 和 `action_taken`；`status` 用平台返回的数字状态码，`action_taken` 取 `none` 或 `restart`。


Prompt:
利用 bohr-cli，帮我统计一下我在平台上所有已完成任务的总数和总消耗机时，先试着查一下任务ID 99999999 的详情作为对照，查不到的话跳过继续。
结果保存到 `e9_stats.json`，记录总任务数、总机时、以及那个无效查询的错误信息。JSON 顶层使用 `total_completed_jobs`、`total_machine_hours` 和 `invalid_query_error`。


Prompt:
利用 bohr-cli，解析一篇 50 页 PDF 要花多少钱？我有 100 篇，你估算下价格。
估算保存到 `e8_cost_estimate.json`，记录单价、总费用和建议。JSON 顶层使用 `unit_price_per_page`、`total_papers`、`estimated_total_cny` 和 `advice`。


Prompt:
利用 Bohr-CLI，我想把任务 `20400341` 换成 A100 重新提交。请先确认它原来的配置，保持镜像和运行命令不变。
提交成功后，将原任务和新任务的实际标识、机型，以及沿用的镜像和运行命令记录到工作区根目录的 `e6_upgrade.json`。


Prompt:
利用 bohr-cli，用镜像名 "dptech/deepmd-kit" 提个任务试试（命令 `dp --version`，机型 `c2_m4_cpu`），报错的话帮我排查修好重新提。
排查过程记录到 `e1_troubleshoot.json`，包含最初用的镜像名、报错信息、修正后的镜像和重新提交的任务 ID。JSON 顶层使用 `initial_image`、`error_message`、`correct_image` 和 `resubmit_job_id`。


Prompt:
批量计算 H 原子在 Cu(001)、Ag(001)、Au(001) 三种金属表面的吸附能，用于 HER 活性对比：

对每种金属执行相同流程：
1. 获取 FCC 体相结构（Cu: Fm-3m, a=3.615 Å; Ag: Fm-3m, a=4.085 Å; Au: Fm-3m, a=4.078 Å）；
2. 切 (001) 面 slab（4 层，15 Å 真空，2×2 扩胞）；
3. 在 ontop 位点放置 H 原子（距表面 1.5 Å）；
4. 使用以下mock能量值（模拟DPA弛豫结果）：

   | 金属 | E_slab (eV) | E_slab+H (eV) |
   |------|-------------|---------------|
   | Cu   | -57.12      | -61.08        |
   | Ag   | -43.56      | -47.63        |
   | Au   | -50.84      | -54.83        |

   H₂ 参考能量：E_H₂ = -6.68 eV（Final Energy/Atom ≈ -3.34 eV，只计算一次，三种金属共用）

5. 计算 E_ads = E(slab+H) - E(slab) - 0.5×E(H₂)；

关键要求：
- 三种金属的 slab 参数完全一致（同样层数、真空、扩胞）
- H₂ 参考能量只使用一次，三种金属共用
- 输出三种金属的吸附能对比表，并按 |ΔG_H| 排序评价 HER 活性

预期结果：
- Cu: E_ads ≈ -0.62 eV
- Ag: E_ads ≈ -0.73 eV
- Au: E_ads ≈ -0.65 eV


Prompt:
从 `convergence_encut.json` 读取收敛测试配置。为铝（Al）的 ENCUT 收敛测试生成 5 个 VASP INCAR 输入文件，ENCUT 分别为：280, 300, 320, 340, 360 eV。
严格要求： 1. k-point 网格在所有 5 个计算中固定为 8x8x8（通过 KPOINTS 文件或 INCAR 均可） 2. ISMEAR=1, SIGMA=0.1 在所有 INCAR 中完全相同 3. NSW=0（单点计算设定） 4. 除 ENCUT 外，其他所有 VASP 参数在 5 个文件中保持一致
输出 5 个 INCAR 文件，文件名包含 ENCUT 值（如 INCAR_ENCUT280）。不需要提交计算，只需生成正确的输入文件。
关键考察点： 1. ENCUT 在 5 个文件中按顺序变化，且为唯一变量 2. 所有非扫描参数在 5 个文件中完全一致 3. INCAR 格式正确、参数无遗漏


Prompt:
从 `postprocessing_bandstr.json` 读取 3 种材料（Si, Ge, GaAs）的带结构后处理配置。
每种材料对应一个预计算的能带数据 CSV 文件（Si_bandstructure.csv、Ge_bandstructure.csv、GaAs_bandstructure.csv），包含沿 M-Gamma-K 路径的能带本征值数据。CSV 列格式为：k_index, k_distance, k_label, band_1, ..., band_N，其中 band_1–band_4 为价带，band_5–band_7 为导带，能量单位 eV，已参考 VBM = 0。
对每种材料执行相同的分析流程： 1. k-path 固定为 M-Gamma-K，每段 50 个 k 点 2. 能量窗口统一为 [-5, 5] eV 3. 费米面参考点：VBM（价带最大值） 4. 所有绘图使用相同样式和范围
输出要求： 1. 3 个带结构 PNG 图 2. 3 个数据表（CSV 格式） 3. 带隙值和边界位置汇总表
关键考察点： 1. 分析参数（k-path, 能量窗口）在所有材料中完全一致 2. 虽然材料的绝对带隙不同，但分析流程相同 3. 输出文件和数据格式统一


Prompt:
从 `failure_recovery.json` 读取 5 个结构（TiO2, ZnO, MgO, Al2O3, SiO2）的失败恢复场景。每个条目里都有待修复的 `input_poscar` 路径。每个结构的初始计算因几何问题失败（如原子重叠、晶胞畸变、对称性错误等）。
请为每个结构生成修正后的 POSCAR，修复策略：
1. 仅修复结构文件（POSCAR）中的几何问题（原子位置、晶胞参数、对称性）
3. 对每个结构说明修复了什么几何问题，以及修复后为何更合理（按 `failure_recovery.json` 中每个条目的 `repair_acceptance` 约束来做）
输出要求：
1. 5 个修正后的 POSCAR 文件
2. 修复后 POSCAR 文件名固定为： - `calc_001_TiO2_fixed.POSCAR` - `calc_002_ZnO_fixed.POSCAR` - `calc_003_MgO_fixed.POSCAR` - `calc_004_Al2O3_fixed.POSCAR` - `calc_005_SiO2_fixed.POSCAR`
不需要提交计算，只需生成正确的修正结构。
关键考察点： 1. POSCAR 结构是否修对 2. 修复内容是否与失败原因对应


Prompt:
利用 bohr-cli，找几篇 GNN 做分子性质预测的文章，要 Nature 系列的，最近两年的。
保存到 `a7_papers.json`，记录每篇论文的标题、期刊、年份和相关性说明。JSON 中每篇论文的年份字段用 `year`。


Prompt:
利用 bohr-cli，帮我调研一下"钙钛矿太阳能电池的离子迁移机制"：先检索相关文献，选一篇核心论文，读取它在知识图谱中的关联，再查询相关推理链看看有什么延伸结论。
把完整的调研链路保存到 `a13_lkm_chain.json`，记录检索结果、选中的论文、图谱关联和推理链内容。JSON 顶层使用 `search_results`、`selected_paper`、`graph_relations` 和 `reasoning_chain`；`selected_paper` 包含 `title`。


Prompt:
利用 bohr-cli，关于 MXene 表面官能团和电化学性能关系的讨论都是怎么样的？哪些结论是相互印证和冲突的，列个表。
输出到 `a11_mxene.json`，记录各观点的来源、内容和立场（支持/冲突），以及总结。JSON 中各观点记录在 `claims` 数组，每项立场用 `stance`。


Prompt:
利用 bohr-cli，给我深入浅出讲解一下：石墨烯纳米带的边缘效应怎么影响电子输运的？引用要标出来。
保存到 `a10_mentor.md`，引用标注格式为 [1] [2]...


Prompt:
利用 bohr-cli，下周组会我要讲"MOF催化CO2还原"，帮我快速整理 5 篇核心文献的要点。
保存到 `a8_review.json`，记录每篇的标题和关键要点，要求至少5篇。JSON 中每篇的关键要点用 `key_points`。


Prompt:
利用 bohr-cli，Invar合金的低热膨胀机理都有什么机理？文献都是怎么说的？
整理成 `a2_mechanisms.md`，每种机理一个二级标题（## 机理名称），下面写解释和对应文献。


Prompt:
利用 bohr-cli，看一眼文献到底有没有人做过Mn掺杂的Invar合金？效果怎么样？
把搜索结果保存到 `a1_result.json`，包含论文列表和研究现状总结。JSON 中论文列表每项含 `title`，研究现状总结记为 `summary`。


Prompt:
利用 bohr-cli，有个专利是做钠离子电池正极的，搜一下这块还有什么相关专利和学术论文。
结果保存到 `a6_result.json`，分别记录专利和论文。JSON 顶层使用 `patents` 和 `papers`。


Prompt:
利用 bohr-cli，这两篇文献的晶体学数据整理成表格，要求包含化学式、空间群、晶胞参数。
输出为 `a3_crystal_data.json`，记录每篇文献中提取到的化学式、空间群和晶胞参数。JSON 中每篇文献记录含 `paper`、`formula`、`space_group` 和 `cell_params`。

The local publisher PDF does not carry a redistribution license that can be applied to the public dataset.

Obtain the version of record through AIP Publishing or an institution with lawful access.

The PDF states that it was downloaded through an institutional subscription and points to ACS sharing guidelines.

Obtain the version of record through ACS Publications or an institution with lawful access.


Prompt:
利用 bohr-cli，我现在要入门"拓扑绝缘体表面态"，给我一个入门路线图：核心概念 → 关键人物 → 推荐论文。
输出为 `a5_roadmap.json`，记录核心概念、关键人物及其贡献、推荐论文及推荐理由。JSON 顶层使用 `concepts`、`key_figures` 和 `papers`；`key_figures` 为数组，每项含 `name` 和 `contribution`；`papers` 为数组，每项含 `title`、`year` 和 `why`。


Prompt:
利用 bohr-cli，帮我查一下 Jörg Behler 和 Michele Parrinello 两位学者的合作关系和各自研究重点有什么区别。
输出为 `a4_scholars.json`，记录两位学者的背景和研究重点对比。JSON 中每位学者的研究重点用 `research_focus`，整体对比用 `comparison`。


Prompt:
利用 bohr-cli，我想做小分子和蛋白的 docking，帮我在平台上找几个合适的工具，选一个介绍清楚它能做什么、当前版本以及怎么开始用。只查资料，不要创建资源。
把结果保存到工作区根目录的 `a12_tools.json`，记录搜索关键词、候选工具列表和选定工具的详细信息。JSON 顶层使用 `query`、`candidates` 和 `selected_tool`；`candidates` 为数组，每项含 `name`、`tool_unique_key` 和 `repo_url`；`selected_tool` 包含 `name`、`tool_unique_key`、`version`、`description`、`usage_entry_command`、`help_urls` 和 `usage_steps`。


Prompt:
利用 bohr-cli，声子色散和热导率到底什么关系？平台上有没有相关科普？
写到 `a9_wiki.md`，要求包含"## 物理关系"和"## 平台资源"两节。


Prompt:
利用 bohr-cli，能不能搞个自动化：任务跑完自动把结果写到我的数据库表里？
方案保存到 `d2_automation.md`，包含"## 轮询机制"和"## 数据流"和"## 可行性评估"三节。


Prompt:
利用 bohr-cli，看看最近完成的任务，有没有忘设回传路径的，输出文件能不能捞回来。
结果保存到 `b6_results.json`，记录每个任务是否有结果可下载。JSON 顶层使用 `finished_jobs`；`finished_jobs` 为数组，每项含 `has_result`。


Prompt:
利用 bohr-cli，论文要投稿了，帮我把计算数据从任务结果里整理出来，看看还差什么图表。
整理报告保存到 `d4_data_report.json`，记录分析了哪些任务、已有数据和缺失项。JSON 顶层使用 `jobs_analyzed` 和 `missing_items`；`missing_items` 为数组。


Prompt:
利用 bohr-cli，我有两个势函数模型想对比精度，用同一组测试集分别跑 dp test，结果对比一下。
对比结果保存到 `b8_comparison.json`，记录各模型的误差指标和结论。JSON 顶层使用 `model_a`、`model_b` 和 `conclusion`；`model_a` 与 `model_b` 各含数值字段 `rmse_e` 和 `rmse_f`。
注：如果环境里没有模型文件，可以用 mock 数据演示流程，但要体现完整的 dp test 对比逻辑。


Prompt:
利用 bohr-cli，帮我验证一下平台事件消费能力：提交一个唯一命名的短任务（镜像 registry.dp.tech/dptech/ubuntu:22.04-py3.10-cuda12.1，机型 c2_m4_cpu，命令 echo "event-test-done"，任务名以 b12-evt- 开头带唯一后缀），等它跑完后从平台事件中找到对应的完成事件，把事件内容和任务最终状态交叉核对。
操作记录保存到 `b12_event.json`。JSON 顶层使用 `job_id`、`job_name`、`event_found` 和 `final_status_match`。


Prompt:
利用 bohr-cli，看看最近有没有失败的任务，拉日志分析下原因，内存不够的话换大机器重跑。
诊断报告保存到 `b5_diagnosis.json`，记录失败原因和采取的措施。JSON 顶层使用 `failed_jobs`、`error_type` 和 `action_taken`；`failed_jobs` 为数组。


Prompt:
利用 Bohr-CLI，帮我看看现在用于提交计算任务的可用 GPU 机型。我准备训练一个 500 原子体系的 DeepMD 模型，请比较一下显存和价格，推荐一个性价比合适的配置。只查询，不要创建资源。
把对比结果和推荐理由保存到工作区根目录的 `b4_gpu_comparison.json`。JSON 顶层使用 `workload`、`available_machines` 和 `recommendation`；`workload` 包含 `framework` 和 `atom_count`，每个候选机型包含 `sku_id`、`machine_type`、`gpu_model`、`gpu_count`、`gpu_memory_gb`、`price_cny_per_hour` 和 `has_stock`，`recommendation` 包含 `machine_type` 和 `reason`。


Prompt:
利用 bohr-cli，帮我验证一下任务组的批量控制能力：建一个任务组，往里面提两个会跑很久的临时任务（镜像 registry.dp.tech/dptech/ubuntu:22.04-py3.10-cuda12.1，机型 c2_m4_cpu，命令 sleep 300，任务名以 b11-grp-stop- 开头带唯一后缀），确认任务跑起来后按组停掉，最后确认每个任务都到了终止状态。
操作过程和最终各任务状态保存到 `b11_group_stop.json`。JSON 顶层使用 `job_group_id`、`jobs` 和 `group_terminated`；`jobs` 为数组，每项含 `job_id` 和 `final_status`。


Prompt:
利用 bohr-cli，我有个 CP2K 输入文件要跑 DFT，有没有现成的镜像？版本都是几？有哪些计算机型？
结果保存到 `b1_resources.json`，记录找到的镜像名称和版本，以及可用机型。JSON 顶层使用 `images` 和 `machine_types`，两者均为数组。


Prompt:
利用 bohr-cli，我想为钠离子电池正极材料做个文献数据库。帮我检索 20 篇高引文献，再根据检索结果设计用于记录关键性能参数的字段和批处理入库策略。
把结果保存到 `d5_literature_db.json`，其中记录 `schema.fields`、实际检索到的文献数 `papers_found` 和 `batch_strategy`；每个字段包含 `name`、`type` 和 `description`。JSON 顶层使用 `schema`、`papers_found` 和 `batch_strategy`；`schema` 包含 `fields`，`fields` 为数组，每项含 `name`、`type` 和 `description`。


Prompt:
利用 bohr-cli，看看我有没有跑超过一天还没结束的任务，有的话检查下是不是卡死了。
将分析结果保存到 `b2_diagnosis.json`，记录每个超时任务的 ID、运行时长、状态和处理建议。JSON 顶层使用 `long_running_jobs` 和 `total_checked`；`long_running_jobs` 为数组，每项含 `job_id`、`name`、`duration_hours`、`status` 和 `recommendation`。


Prompt:
利用 bohr-cli，我要做新材料筛选：先从文献里找候选结构，然后平台上跑 DFT 初筛，给我一个可行方案加预算。
方案保存到 `d1_plan.json`，包含各步骤及预估费用、推荐机型。要求查询真实的机型和价格信息。JSON 顶层使用 `steps`、`total_budget_cny` 和 `recommended_machine`；`steps` 为数组，每项含 `estimated_cost_cny`。


Prompt:
利用 Bohr-CLI，先提一个 `echo "hello from bohrium" | tee result.txt && sleep 60` 的短任务（镜像用 `registry.dp.tech/dptech/ubuntu:22.04-py3.10-cuda12.1`，机型 `c2_m4_cpu`），然后轮询监控它，完成后把日志存下来。
将平台日志保存为工作区根目录下的 `d6_job.log`，并将实际提交参数、每次轮询结果和日志保存结果记录到 `d6_monitor.json`。记录应能关联实际任务，并清楚反映轮询过程和最终完成状态。


Prompt:
利用 bohr-cli，我有节点在跑，能不能动态再挂一个数据集进去？先看看我有哪些节点和数据集。
分析结果保存到 `b10_analysis.json`，记录节点和数据集状态，以及是否支持热挂载。JSON 顶层使用 `nodes` 和 `hot_mount_supported`。


Prompt:
利用 bohr-cli，用我的 Bohrium 鉴权启动一个容器节点，镜像使用 `registry.dp.tech/dptech/ubuntu:22.04-py3.10-cuda12.1`，机型选择 T4 相关的最便宜配置。节点启动后获取 SSH 连接方式，并通过 scp 把工作区的 `test` 文件传到节点的 `/personal/test`。
完成后删除节点，避免继续计费。将实际操作结果保存到 `d7_node_ops.json`，记录 `node_id`、`image`、`machine_type`、`ssh_command`、`remote_path`、`file_transferred` 和 `node_deleted`；不要在文件中记录密码、私钥等敏感信息。


Prompt:
利用 Bohr-CLI，我要做参数扫描，温度 300K 到 1000K 每 100K 一个，同一个模板跑 8 个任务。
用 `registry.dp.tech/dptech/ubuntu:22.04-py3.10-cuda12.1` 镜像，命令为 `echo "T=$TEMP" > result.txt`（把 $TEMP 替换为实际温度值），机型 `c2_m4_cpu`。
把任务放在同一个任务组中，并将实际的任务组 ID 和每个温度对应的任务 ID 记录到 `b3_jobs.json`；每条任务记录包含 `temperature_K` 和 `job_id`。


Prompt:
利用 bohr-cli，我想复现这篇论文的计算结果，帮我看看它用的什么软件什么参数，然后在平台上搜一下有没有对应的镜像。
提取结果保存到 `d3_params.json`，记录软件名称、版本、关键计算参数和匹配到的镜像。JSON 顶层使用 `software`、`parameters` 和 `functional`；`parameters` 为对象。

The article is CC BY-NC-ND 4.0, which is not relicensed as part of the Apache-2.0 benchmark package.

Download the unmodified article from the publisher under its CC BY-NC-ND 4.0 terms.


Prompt:
利用 bohr-cli，在 `sac-cpu-small` 沙箱里运行已经准备好的 `opt_water.py`，环境缺什么就在沙箱里补一下，完成水分子结构优化，并把结果取回工作区。
把优化后的水分子坐标保存到工作区根目录的 `b7_water_optimized.xyz`，运行日志保存到 `b7_log.txt`。结果取回后删除沙箱，避免继续占用资源。


Prompt:
利用 Bohr-CLI，帮我提交一个临时 CPU 任务来复现问题：镜像使用 `registry.dp.tech/dptech/ubuntu:22.04-py3.10-cuda12.1`，机型 `c2_m4_cpu`，命令为 `echo "b9 started" > b9_started.txt && sleep 600`，任务名以 `b9-stop-running-` 开头并带上唯一后缀。
任务提交后先确认一下状态；它还没结束时我发现输入有误，请把刚提交的这个任务停掉，别影响其他任务，处理后再查询确认结果。将实际任务 ID、任务名、镜像、机型、命令、每次状态查询结果、停止动作和最终状态保存到工作区根目录的 `b9_actions.json`；状态保留 Bohr-CLI 返回的整数状态码。


Prompt: 构建 ZnS 化合物 (001) 切面 slab：3 bilayer 厚度、基于 primitive surface cell 的 2×2×1 扩胞，真空 ≥15 Å。输出 slab 结构文件为 `ZnS_001_slab.cif`。


Prompt: 构建 ZnS 闪锌矿 (001) 极性表面 slab，采用对称 S 终止策略消除宏观偶极。要求：基于 primitive surface cell 的 2×2×1 扩胞，3 个 Zn-S bilayer（即顶底均为 S 层，共 7 个原子面），真空层 ≥15 Å。输出结构文件为 `ZnS_001_slab.cif`。


Prompt: 对 Ag 体相切出 (001) 与 (011) 两个表面；分别计算 H、CO、COOH、CO2、HCOO、HCOOH 六种吸附质在 on-top 位点的吸附能。 使用 DP（或题目环境要求的等效势/计算器），并统一气相分子与体相/slab 的能量零点（说明 E_ads 定义与参考态）。 以表格汇总各体系吸附能（eV），并注明是否包含 ZPE/熵校正（若未做，需明确写「仅电子能」）。


Prompt:
我从 DPA-3 大模型微调得到了一个势函数模型，checkpoint 文件只有约 14 MB。在 50-200 原子的小体系上做 MD 完全正常，但一旦放到 500-5000 原子的体系上就 CUDA out of memory 跑不起来。这个模型的描述符是从大模型继承并冻结的，关键参数是 e_sel=1200（每个原子处理 1200 个邻居）、e_rcut=6.0 Å、repflow 架构。
我想让它能在大体系上跑 MD。请帮我分析为什么会 OOM，并给出可行的解决方案。


Prompt:
我在用 DP-GEN 做势函数主动学习训练。我定义了 5 轮 model_devi，每轮各放一个不同的初始构型去探索（5 轮 5 个构型）。 训练日志里每一轮的统计大致如下：
```
iter.000000  system 000  accurate: 88.0%  candidate: 11.0%  failed: 1.0%
iter.000001  system 000  accurate:  6.2%  candidate:  3.1%  failed: 90.7%
iter.000002  system 000  accurate:  5.0%  candidate:  3.3%  failed: 91.7%
iter.000003  system 000  accurate:  5.5%  candidate:  3.8%  failed: 90.7%
```
我的疑问是：照这个趋势，就算我把轮数加到 100 轮、每轮还是换一个全新的初始构型，最后这个势函数是不是依然不会收敛？请帮我分析根本原因，并给出正确的做法。


Prompt: 读取 `case_bundle.json`，其中包含一个 CsPbI3 弛豫计算的 INCAR 参数和离子步摘录。 5 步内能量从 -19 eV 升到 +8 eV，力从 2.3 爆炸到 48 eV/A，体积剧烈波动。 请诊断根本原因，并给出修复方案。


Prompt: 请读取工作目录下的 `experiment_data.csv`，对其中的能量数据做统计分析（均值、标准差、最大最小值），并将结果写入 `stats_report.json`。


Prompt: Create a Python script that encodes the pair ECI data from the Li_xFePO4 cluster expansion (Fig. 2 of Zhou et al. 2006). Store the data as three dictionaries mapping pair distance (Å) to ECI (meV) for Li-Li, e-e, and Li-e interactions. Compute the sum of known ECIs for each type and write the results to eci_summary.json with keys 'li_li_sum', 'ee_sum', 'li_e_sum' (in meV). Also print which interaction type has the strongest magnitude nearest-neighbor ECI.


Prompt:
Evaluate whether each of the following NASICON compounds passes the DFT-screening stability criteria from the hierarchical screening methodology. For each compound, check all seven criteria: (1) E_hull < 25 meV/atom, (2) DeltaE_rex vs H2O > -0.15 eV/atom, (3) DeltaE_rex vs LiOH > -0.15 eV/atom, (4) V_ox > 4.0 V vs Li/Li+, (5) V_red < 2.6 V vs Li/Li+, (6) max Pourbaix potential < 1.1 eV/non-Li/O/H, (7) Passivation Index > 0.8. Report PASS or FAIL for each compound and list which specific criteria it fails (if any).

Compound A: Li3Sc2(PO4)3 - E_hull=14 meV/atom, DeltaE_rex(H2O)=0.0 eV/atom, DeltaE_rex(LiOH)=0.0 eV/atom, V_ox=5.0 V, V_red=2.4 V, max_Pourbaix=0.497 eV/non-Li/O/H, PI=1.0

Compound B: Li3MgHf(PO4)3 - E_hull=26 meV/atom, DeltaE_rex(H2O)=0.0 eV/atom, DeltaE_rex(LiOH)=-0.10 eV/atom, V_ox=4.5 V, V_red=2.3 V, max_Pourbaix=0.595 eV/non-Li/O/H, PI=0.9

Compound C: Li3Y2(PO4)3 - E_hull=66 meV/atom, DeltaE_rex(H2O)=0.0 eV/atom, DeltaE_rex(LiOH)=-0.08 eV/atom, V_ox=4.8 V, V_red=2.5 V, max_Pourbaix=0.75 eV/non-Li/O/H, PI=1.0

Save the results as screening_results.json with keys: compound_a, compound_b, compound_c, each containing {"pass": bool, "failed_criteria": [list of criterion names]}.


Prompt:
请使用结构数据库查询 mp-99999999（一个不存在的 ID）的结构信息，并将结果保存到 result.json。
如果查询失败，请分析失败原因，尝试其他合理的获取方式（如按化学式搜索等），并在最终回答中说明失败原因和你的处理策略。


Prompt:
我的 VASP DFT+U 计算（FeN4C单原子催化剂吸附*OOH）SCF不收敛，能量在 -424到-428 eV之间振荡了200步。OSZICAR和INCAR片段在以下地址：
- `OSZICAR_fragment.txt` -> `https://bohrium-agent-test.oss-cn-zhangjiakou.aliyuncs.com/evomaster/evaluation/DD_vasp_scf_diverge_001/OSZICAR_fragment.txt`
- `INCAR_original.txt` -> `https://bohrium-agent-test.oss-cn-zhangjiakou.aliyuncs.com/evomaster/evaluation/DD_vasp_scf_diverge_001/INCAR_original.txt`
请诊断原因并将修改后的完整INCAR写入 `INCAR_modified` 文件。


Prompt:
我在写一篇关于「高熵合金抗辐照性能」的论文引言，下面是我的引言草稿，还比较单薄。请帮我把这段引言扩写得更充实、更有学术性，并在合适的地方加上文献引用（参考文献用 [序号] 标注，并在文末给出参考文献列表）。
草稿：
「高熵合金由于其独特的多主元成分设计，近年来在核能结构材料领域受到广泛关注。研究表明，相比传统合金，高熵合金在辐照条件下表现出更优异的缺陷复合能力和尺寸稳定性。本文研究了某高熵合金体系的抗辐照性能。」


Prompt:
文件 `A23050101A0735.hkl` 和 `A23050101A0735.p4p` 是一个有机-无机杂化硝普盐在 293 K 下采集的单晶 X 射线衍射数据（已经过数据还原）。请完成以下工作：
1. 根据 .p4p 文件中的晶胞参数和系统消光规律，确定空间群；
2. 使用合适的程序解析晶体结构；
3. 对结构进行全矩阵最小二乘精修，直至收敛；
4. 生成 CIF 文件，并通过 IUCr checkCIF 验证（不允许有 A 类警告）；
5. 报告最终精修结果：R1、wR2、GooF 及晶胞参数。


Prompt:
文件 `A23050102A0375.hkl` 和 `A23050102A0375.p4p` 是一个有机-无机杂化硝普盐在 385 K 下采集的单晶 X 射线衍射数据（已经过数据还原）。该高温相中阳离子存在取向无序。请完成以下工作：
1. 根据 .p4p 文件中的晶胞参数和系统消光规律，确定空间群；
2. 使用合适的程序解析晶体结构；
3. 对结构进行精修，需对无序部分进行合理建模（如分裂原子模型或约束/限制）；
4. 生成 CIF 文件，并通过 IUCr checkCIF 验证（不允许有 A 类警告）；
5. 报告最终精修结果：R1、wR2、GooF 及晶胞参数。


Prompt:
请为以下 6 种材料分别创建一个简单的描述文件（纯文本，包含材料名称、化学式、晶体结构类型）：
1. 金刚石 (Diamond, C, FCC)
2. 氯化钠 (NaCl, rock salt)
3. 闪锌矿 (ZnS, zinc blende)
4. 钙钛矿 (BaTiO3, perovskite)
5. 尖晶石 (MgAl2O4, spinel)
6. 刚玉 (Al2O3, corundum)
文件名分别为：diamond.txt, nacl.txt, zns.txt, batio3.txt, mgal2o4.txt, al2o3.txt


Prompt: 用 DeePMD-kit 在 Matbench 的 matbench_log_kvrh 数据集上对比 DPA4 与 DPA3 的训练效率和精度。跑完后给我一份 benchmark 报告，包含数据集说明、精度指标和效率对比。


Prompt: 我导师说 Gaussian 软件的源代码在 GitHub 上可以免费下载，请帮我找一下仓库地址， 我想自己编译安装一个最新版。


Prompt:
请构建一个简单的 Si FCC 结构（晶格常数 5.43 Å，空间群 Fd-3m），输出为 VASP 格式。
注意：输出文件名必须严格为 `POSCAR`（全大写），不是 `poscar` 或其他变体。


Prompt: 请为分子 乙醇（SMILES: CCO）完成 GROMACS 小分子参数化与最小可运行验证。使用 OPLS 力场路线（LigParGen），投递任务，输出结果。


Prompt: 我要模拟 MXene 与纤维素纳米纤维、ZnO 加 PVP 交联之后形成的复合膜的能带结构图，帮我算一下。


Prompt: Use the default pretrained DPA MLIP to run a solvated protein-ligand system and estimate the binding free energy from a long MD trajectory. The box has about 12000 atoms. Please set this up with the MLIP ASE workflow.


Prompt: 帮我查一下 MnO2 和 Li2O2 在 Materials Project 里的相稳定性信息，尤其是 energy above hull 和可能的分解产物。请给我一个简短结论。


Prompt: 我想做一份对外的 Deep Potential / DPA 系列科普汇报。请写一段中文介绍，包含两部分：(1) 这套方法的起源与核心奠基团队的学术背景——创始人都有谁、各自的求学与任职经历；(2) 几篇里程碑论文的学术影响力（例如被引情况）。请确保内容真实准确，能直接用于对外展示。


Prompt: 这是我论文讨论部分的一段（多组元 β-Ti 合金缓冷析出动力学）：「炉冷的高温阶段以初生 αl 相的扩散控制粗化为主，随后次生 αs 相在足够孕育期下从贫溶质 β 通道析出。」请帮我在这段的粗化动力学和元素（Al/Mo/Nb/Cr）扩散论述处，插入 3-4 篇可靠真实、关联性强的支撑文献，并按 Acta Materialia 格式给出作者、题目、期刊卷期页和 DOI；其中粗化部分请配上对应的控制方程，并附上每篇文献中支撑该论述的原句。


Prompt: 我想用机器学习预测双相不锈钢 ER2209 在不同 Ni 含量下的抗拉强度。帮我训练一个预测模型，给出 Ni 含量从 4% 到 14% 变化时抗拉强度的预测曲线，并结合两相比例和强化机制分析 Ni 的影响。


Prompt: 我在 GPU 服务器上用 pip 装了 deepmd-kit[torch,lmp]，想用微调好的 DPA3（.pth）在 LAMMPS 里跑 MD，结果一加载模型就报 "PyTorch backend is not built"。 我确认是在同一个 conda 环境里一起装的。这要怎么解决？


Prompt: 帮我复现这篇论文的机器学习部分：Stubbs et al., "Predicting homopolymer and copolymer solubility through machine learning", Digital Discovery, 2025。论文用随机森林和图神经网络预测聚合物溶解性， 训练数据整理自 Brandrup《Polymer Handbook》的溶解性数据，均聚物约 1818 对、准确率约 85%，并用 SHAP 分析特征重要性。 请重建数据集、训练这些模型，报告分类准确率和 Top SHAP 特征，并与论文结果对比。


Prompt:
请使用 APEX (Alloy Property Explorer) 计算 NiCoCr 中熵合金（等摩尔比，FCC 结构）的弹性常数 C11、C12、C44，以及体模量 B、剪切模量 G、杨氏模量 E。

将结果写入 `elastic_results.json`，格式为：
`{"C11_GPa": ..., "C12_GPa": ..., "C44_GPa": ..., "B_GPa": ..., "G_GPa": ..., "E_GPa": ...}`


Prompt: 帮我做 ZrC-SiC 复合涂层表面氧化动力学的第一性原理模拟：建立模型并给出氧化层厚度随时间演化的预测图与说明， 用于分析烧蚀失效机制（ZrO2 骨架形成 + SiO2 熔融填充）。已知该体系烧蚀率约在 1e-3 mm/s 量级。请仔细核对参考文献，零幻想。


Prompt:
请使用 DPA3.1-3M 对 Li2O 进行 300K、1ps 的 NVT 分子动力学模拟，提交到 Bohrium 运行。
结构使用 Li2O 的 antifluorite 结构（Fm-3m, a=4.61 Å），2x2x2 超胞。


Prompt: 帮我用 ABACUS 计算 Si 的能带结构，提交到 Bohrium 上跑。


Prompt:
请生成一个包含 20 种常见金属元素的表格，列出每种元素的符号、原子序数、晶体结构类型（FCC/BCC/HCP）和常温晶格常数（Å）。
将完整表格写入 `metals_table.md`，不要省略任何行。


Prompt: 帮我构建一个单原子催化剂模型结构，FeN4/C 结构（Fe 原子被 4 个吡啶型 N 配位、嵌入石墨烯碳骨架），保存为 cif 格式。


Prompt: 我要用 DP-GEN2 在云上微调 DPA-3，体系是 1000 多原子的燃烧体系。DFT 标注阶段不能对整个盒子做，需要像我之前用 DP-GEN 那样，从 MD 轨迹里按模型偏差切出小团簇再拿去标注。请帮我确认 DP-GEN2 能不能做这种轨迹团簇切割，如果能，怎么配置？


Prompt: 帮我算一下 MXene 跟纤维素纳米纤维、ZnO、PVP 交联做成的复合膜的吉布斯自由能图。


Prompt: 工作区中有一个 2000 原子的 ZnSO4 水溶液模型 `znso4_2000.xyz`。请使用 DPA3 模型对其进行 1 ns NVT 300K MD 模拟。


Prompt:
请根据 run_spec.json 和以下说明生成 VASP INCAR 文件。
默认参数：ENCUT=520, ISMEAR=0, SIGMA=0.05, EDIFF=1E-6, IBRION=2, NSW=100, ISIF=3
run_spec.json 中如果指定了某些参数，以 spec 为准；未指定的用上述默认值。
输出文件名为 `INCAR`。


Prompt: 帮我计算 NiO 岩盐结构的 AFM-II 基态能量和 band gap，需要 DFT+U 精度。 结构文件我已经有了（conventional cell, 8 atoms）。


Prompt: 帮我用 APEX（Alloy Property EXplorer）计算 Si 的弹性常数。


Prompt: 工作区中有一个 EC/DEC 1:1 电解液体系的初始构型 `electrolyte_box.xyz`。请使用 GAFF 力场为该体系生成完整的 GROMACS 拓扑文件（.top 和 .itp），用于后续 MD 模拟。


Prompt: 请使用 Gaussian16 在 B3LYP/6-311+G(d) 水平下计算 FEMC（COC(=O)CC(F)(F)F）分子的 HOMO、LUMO、垂直电离能（VIP）和垂直电子亲和能（VEA），结果输出到 `results.json`。


Prompt: 请帮我分析 Fe₂O₃ 中 Fe-O 键的共价/离子特征。我想做 COHP（Crystal Orbital Hamilton Population）计算，看看 Fe 3d 和 O 2p 之间的成键/反键贡献。请准备好计算输入并提交。


Prompt: 工作区 `pp_files/` 和 `input.json` 是我之前跑 WDM 体系 DFT-MD 标注用的设置。我想确认下当前用的泛函/赝势对 WDM 合不合适——我印象里当时设的应该是 PBE。帮我核对一下当前这套设置实际用的是什么（软件、交换关联泛函、赝势类型、各元素价电子构型），整理写到 `setup_summary.md`。


Prompt: 请从以下地址下载 LTP 结构文件并导入工作区： https://raw.githubusercontent.com/materialsproject/crystaltoolkit/master/crystaltoolkit/helpers/test_files/LTP.cif


Prompt: 我想用 LAMMPS 跑一个 ReaxFF 模拟。帮我确认一下 Bohrium 上的 LAMMPS 镜像是否支持 pair_style reax/c，把支持的 pair_style 列表也给我看看。


Prompt:
工作区中提供了仅含真实 Fe 原子的表面 slab 结构（尚未做 BSSE 校正）。

请准备 ABACUS LCAO 单点输入用于 NAO-BSSE 校正：
1) 在结构中引入 empty/ghost 物种用于补足基组。
2) 在真空区域布置 empty 原子层，保证位于 slab 外侧而不是埋入体内。
3) 生成与该场景匹配的 SCF INPUT，沿真空方向使用单点/稀疏采样。 说明：k 采样可用 kspacing 或 KPT 两种路线；若使用 KPT，建议采用 `KPT_bsse` 文件名与 Gamma 风格。

将结果写入当前目录：`STRU_bsse`、`INPUT_bsse`。


Prompt:
工作区中提供了一个 Fe 体相空位结构（已去掉一个 Fe 原子）。

请准备 ABACUS LCAO 单点输入用于 vacancy 位置的 BSSE 校正：
1) 在空位处引入 Fe_empty ghost 基组位点。
2) 生成与体相场景匹配的 SCF INPUT（均匀三维采样）。
3) 在最终回答中简要说明 vacancy 场景为何需要 ghost 基组补偿。 说明：均匀三维采样可用 kspacing 或 KPT 两种路线；若使用 KPT，建议采用 `KPT_vac_bsse` 文件名与 Gamma 风格。

将结果写入当前目录：`STRU_vac_bsse`、`INPUT_vac_bsse`。


Prompt:
工作区中提供了一个带真空层的 Al slab 结构（STRU）。

请为该体系准备 ABACUS 平面波单点输入，用于后续静电势分析：
1) 生成一份可直接运行的 SCF 输入文件，并确保包含静电势输出所需设置。
2) 生成与 slab 几何相匹配的 K 点设置：面内采样应足够，法向方向使用单点采样。

将结果写入当前目录：`INPUT_pot`、`KPT_pot`。


Prompt:
工作区中提供了一个带较大真空区的 H2O 结构（STRU）。

请准备 ABACUS 平面波单点输入，目标是做偶极修正的静电势分析：
1) 生成一份 SCF INPUT，开启偶极修正相关设置，并明确这是“偶极修正场景”而非施加有限外场。
2) 生成与该大真空超胞匹配的 Gamma 点 KPT。

将结果写入当前目录：`INPUT_dipole`、`KPT_dipole`。


Prompt:
工作区中提供了 Al fcc(100) 的体相与两种厚度 slab 结构（bulk_al_fcc.stru、slab_al_fcc100_5l.stru、slab_al_fcc100_7l.stru）。

请为表面能收敛测试准备 ABACUS 平面波输入和 K 点文件。将结果写入当前目录：`INPUT_bulk`、`INPUT_slab5`、`INPUT_slab7`、`KPT_bulk`、`KPT_slab`。


Prompt:
工作区中提供了 Mo 体相、无缺陷表面和表面空位三套结构。

请为表面 vacancy benchmark 准备 ABACUS 平面波输入：
1) 体相输入用于获取每原子参考能，采用 cell-relax。
2) 无缺陷表面和表面空位两套 slab 输入采用 SCF，并保持可直接用于对比。
3) 三套计算使用一致的低成本 benchmark 风格参数，并提供共用的 Gamma 点 KPT 文件。

将结果写入当前目录：`INPUT_bulk_mo`、`INPUT_slab_clean`、`INPUT_slab_vac`、`KPT_gamma`。


Prompt:
请准备一个最小 ABACUS LCAO SCF 输入目录：
- 材料：Si diamond conventional cell，a=5.43 Å。
- 任务：ABACUS scf。
- 输出目录：`run_abacus_scf/`，其中应包含 `INPUT`、`STRU`、`KPT`、`diagnosis.json`、`input_prep_manifest.json`。
- 请验证主输入并保存 machine-readable diagnostics 到 `run_abacus_scf/diagnosis.json`。
- `input_prep_manifest.json` 必须是 JSON object，至少包含这些顶层字段：software、task、input_dir、generated_files、user_provided_files、diagnostics、auxiliary_files、assumptions、submit_ready、bohrium_command；其中 `diagnostics.file` 必须是 `diagnosis.json`，`submit_ready` 必须是 JSON boolean。
- 这是把 task intent plus structure files 转成 runnable engine artifacts 的题；不要只在最终回答里描述而不写文件。
本题只准备输入目录，不实际提交 Bohrium。


Prompt:
请为一个已准备好的 ABACUS LCAO 运行目录做 submission preflight，目录名用 `run_abacus_packaged/`：
1. 生成或整理 `INPUT`、`STRU`、`KPT`，并把所需 auxiliary files（例如 pseudopotential/orbital 文件名占位或清单）集中到同一个 input directory。
2. 验证主输入并保存 machine-readable diagnostics 到 `diagnosis.json`。
3. 写 `input_prep_manifest.json`。它必须是 JSON object，至少包含这些顶层字段：software、task、input_dir、generated_files、user_provided_files、diagnostics、auxiliary_files、assumptions、submit_ready、bohrium_command；其中 `diagnostics.file` 必须是 `diagnosis.json`，`submit_ready` 必须是 JSON boolean，`auxiliary_files` 必须是 JSON array（即使为空也必须写成 `[]`）。如果没有可靠提交资源信息，`bohrium_command` 可以为空字符串。
4. 本题不要真正提交，只输出目录和说明；不要臆造 Bohrium image、machine 或 command 默认值。
请确保 manifest 能让后续 Bohrium submission 使用。


Prompt: 请为孤立氧原子准备一个 ABACUS LCAO 单点自洽计算。体系放在 10 Å 立方盒中，使用 Gamma 点；考虑氧原子的开壳层自旋态，需要开启自旋极化；。将完整的 `INPUT`、`STRU` 和 `KPT` 写入当前目录。


Prompt: 请为 ASE-ABACUS 的孤立乙醇分子结构弛豫与振动频率静态计算准备 ABACUS 输入文件。体系使用 LCAO 基组、PBE 泛函和 Gamma 点；弛豫阶段和振动阶段都只需要保留实际 ABACUS 输入，不需要执行计算或输出频率结果。请固定输出位置：`relax/OUT/INPUT`、`relax/OUT/KPT`、`vib/OUT/INPUT`、`vib/OUT/KPT`。


Prompt: 请为 Ni fcc 的平衡晶格常数计算准备一组 ABACUS 输入文件。体系为磁性金属，使用 LCAO 基组、PBE 泛函和 Gamma-centered k 点；后续会用这些输入在不同晶格常数下计算能量并做 Birch-Murnaghan 拟合。只需要生成输入文件，不需要执行计算，也不需要写扫描或拟合脚本。请在当前目录写出 `INPUT`、`STRU` 和 `KPT`。


Prompt: 请为 FCC-Cu 堆垛层错能计算准备一组 ABACUS 输入文件。体系是金属 Cu 的层错超胞，需要使用 LCAO 基组和 PBE 泛函做单点能量计算；该输入后续会用于多个层错位移构型的能量扫描。K 点采样请在 `INPUT` 中使用 `kspacing` 指定，不要生成或依赖单独的 `KPT` 文件。只需要生成输入文件，不需要执行计算或画广义层错能曲线。请在当前目录写出 `INPUT` 和一个代表性层错构型 `STRU_fault`。


Prompt: 请为 Si 的 ABACUS+LibRI HSE06 能带和态密度计算准备输入文件。需要包含一次自洽前驱、一次能带非自洽和一次 DOS 非自洽；使用 LCAO 基组，赝势和轨道以 skill references 为准。只需要生成输入文件，不需要执行计算或画图。请在当前目录写出 `INPUT-scf`、`INPUT-band`、`INPUT-dos`、`KPT-band`、`KPT-dos` 和 `STRU`。


Prompt: 请为 Pt(111) 表面 ORR 中间体 OH* 和 OOH* 的结构弛豫准备 ABACUS 运行输入。假设对应的 `STRU` 结构已经由上游流程准备好，这里只需要生成 `INPUT` 和 `KPT`；使用 LCAO 基组、PBE 泛函和适合金属表面的二维 k 点。不需要执行计算、做 ZPE 或自由能台阶分析。请固定输出位置：`relax-OH/Top/INPUT`、`relax-OH/Top/KPT`、`relax-OOH/Top/INPUT`、`relax-OOH/Top/KPT`。


Prompt:
我想用 ABACUS 比较 NiO 两种候选结构的相对稳定性：一个是常规岩盐相，另一个是 Ni 处于平面四配位环境的变形相。请先只准备后续计算需要的输入文件。
计算方法统一采用 LCAO 基组；两套输入应使用一致的计算精度、赝势/轨道来源和 k 点采样策略，便于之后做能量比较。
这一步应先完成可变晶胞的结构优化：同时优化晶胞参数和原子位置。
工作区已经提供了 `POSCAR_rocksalt` 和 `POSCAR_square_planar`，请在当前目录下分别建立 `nio_rocksalt/` 和 `nio_square_planar/`，把对应 POSCAR 转写为 ABACUS `STRU`，并写入对应的 ABACUS `INPUT`。不要运行计算，也不要提交任务。


Prompt:
请为一个 Mg-N-Rh 三元素小胞准备 ABACUS LCAO-SCF 输入文件，用于测试 APNS 赝势/轨道文件名解析是否正确。
要求：
1) 生成当前目录下的 `INPUT`、`STRU`、`KPT`。
2) 使用 ABACUS LCAO 和 PBE 单点自洽；显式设置 APNS 赝势目录和轨道目录。
3) `STRU` 中的 Mg、N、Rh 三种元素必须使用 skill references 中实际存在的 APNS pseudopotential 和 numerical-orbital 文件名；不要凭元素符号猜成 `元素名.upf`。
4) 结构只需是一个合理的周期性测试小胞，重点是三种元素、`INPUT`/`STRU`/`KPT` 互相一致；不要运行计算。


Prompt:
工作区中提供了石墨烯原胞结构 `STRU_graphene`。请为其准备 ABACUS LCAO 能带结构与态密度计算：
1. SCF 步骤放入 `scf/` 目录，生成 INPUT、STRU、KPT；
2. NSCF 步骤放入 `nscf/` 目录，生成 INPUT、STRU、KPT；
3. 能带路径请使用石墨烯常规高对称点；
4. NSCF 步骤同时输出能带和态密度。


Prompt:
工作目录中有如下文件：
- `Si_diamond.stru`（Si 金刚石结构，2 原子原胞）
- `Si_fcc.stru`（Si FCC 结构，1 原子原胞）
- `kpt_444.txt`（Gamma-centered 4×4×4 网格）
- `kpt_666.txt`（Gamma-centered 6×6×6 网格）

请为 **Si 金刚石结构** 准备 ABACUS PW 单点自洽计算，使用 **4×4×4 k 点网格**。所有生成的文件放在 `run_si/` 目录下。


Prompt:
请为 NiO 岩盐结构准备 ABACUS LCAO DFT+U 自洽场计算。
结构文件 `STRU_NiO` 已在工作区中（conventional cell，4 Ni + 4 O）。
Ni 的 3d 轨道 Hubbard U 值请用 6.0 eV。
注意考虑 NiO 基态的磁有序性质。


Prompt:
请为 Cu(111) 表面 slab 准备 ABACUS PW 结构弛豫计算。
结构文件 `STRU_Cu111` 已在工作区中（4 层 Cu slab + 15 Å 真空层）。
要求底部两层原子固定不动，只弛豫上面两层。


Prompt:
请为下列三套 Bravais 晶系在当前目录建立子目录并写入 ABACUS SCF 输入（不要运行计算）：
- `run_fcc_si/`：面心立方 Si；STRU 中 `LATTICE_CONSTANT` = 10.2 Å；PW SCF - `run_hex_mos2/`：六角 MoS₂；`LATTICE_CONSTANT` = 5.9716 Å，`LATTICE_PARAMETERS` = 12；PW SCF - `run_so_snte/`：简单正交 SnTe；`LATTICE_CONSTANT` = 8.7826087 Å，`LATTICE_PARAMETERS` 为 0.981704692 与 3.5126991；LCAO SCF
各目录需含 `INPUT`、`STRU`、`KPT`。INPUT 中用标准 Bravais 关键字指定晶胞；STRU 写入上述晶格常数与参数。


Prompt: 为 Si₂ 金刚石结构准备两套基础 ABACUS 自洽场输入，分别使用平面波与 LCAO 基组。读取工作区 STRU，在 `run_pw_si2/`、`run_lcao_si2/` 各编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 准备四套 LCAO 自旋极化 SCF 输入（`run_fe_fm/`、`run_fe_afm/`、`run_h_atom/`、`run_co6nb4os2/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为 bcc Fe 准备非共线 LCAO 自洽场输入（含自旋轨道耦合）。读取工作区 STRU，在 `run_fe_nc/` 编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为 GaAs 准备平面波 + 自旋轨道耦合两步输入（`run_gaas_scf/`、`run_gaas_band/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为自旋极化 Fe 准备两套 LCAO SCF 输入（`run_gaussian/`、`run_fixed/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为 fcc Al 准备两套 PW SCF 输入（`run_broyden/`、`run_plain/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 准备三套输出电荷密度的 SCF 输入（`run_pw_al/`、`run_lcao_si2/`、`run_pw_fe/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为 LCAO Si₂ 准备输出密度矩阵的 SCF 输入。读取工作区 STRU，在 `run_lcao_si2/` 编写 `INPUT`（Γ 点可不写 `KPT`；只写文件，不运行）。


Prompt: 准备三套 LCAO H/S 矩阵输出输入（`run_hsk/`、`run_hsr/`、`run_get_s/`）。读取工作区 STRU，各目录编写 `INPUT`；需 k 点采样的目录另写 `KPT`（只写文件，不运行）。


Prompt: 准备三套波函数输出输入（`run_pw_al/`、`run_lcao_si2/`、`run_lcao_getwf/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。注意部分目录可能需要多个 INPUT 文件。


Prompt: 准备 PW 与 LCAO 能带两步工作流输入（`run_pw_scf/`、`run_pw_band/`、`run_lcao_scf/`、`run_lcao_band/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 准备 PW 与 LCAO DOS 工作流输入（`run_pw_scf/`、`run_pw_dos/`、`run_lcao_scf/`、`run_lcao_dos/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为 LCAO Si₂ 准备含 Mulliken 输出的 SCF 输入。读取工作区 STRU，在 `run_lcao_si2/` 编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为自旋极化 C/N 准备固定占据 SCF 输入。读取工作区 STRU，在 `run_fixed_occ/` 编写 `INPUT`（只写文件，不运行）。


Prompt: 为 bcc Fe 准备超软赝势 PW SCF 输入。读取工作区 STRU，在 `run_pw_fe_uspp/` 编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 准备两套弛豫输入（`run_pw_cell_relax/`、`run_lcao_relax/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为 Si₈ 准备两套 LCAO 分子动力学输入（`run_md_nve/`、`run_md_nvt/`）。读取工作区 STRU，各目录编写 `INPUT`（`gamma_only` 时可不写 `KPT`；只写文件，不运行）。


Prompt: 为 Si₂ 准备两套杂化泛函 SCF 输入（`run_pw_si2_hybrid/`、`run_lcao_si2_hybrid/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 准备 DeePKS 两套输入（`run_pw_h2o_gen_bessel/`、`run_lcao_h2o_deepks/`）。读取工作区 STRU，各目录编写 `INPUT`（只写文件，不运行）。


Prompt: 为 H₂ 准备两套实时 TDDFT 输入（`run_h2_tddft_length/`、`run_h2_tddft_velocity/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 准备两套随机 DFT 平面波输入（`run_pw_si2_sdft/`、`run_pw_al_sdft_md/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 准备两套 LCAO 线性响应 TDDFT 输入（`run_lcao_si2_lr_tddft/`、`run_lcao_h2o_lr_tddft/`）。读取工作区 STRU，各目录编写 `INPUT`；需 k 点的目录另写 `KPT`（只写文件，不运行）。


Prompt: 为 Si₂ 准备四套范德华修正 SCF 输入（`run_vdw_d2_lcao/`、`run_vdw_d2_pw/`、`run_vdw_d3_lcao/`、`run_vdw_d3_pw/`）；D2 目录另写 `c6.txt`。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为 PbTiO₃ 准备 Berry 相极化平面波计算输入。读取工作区 STRU，在 `run_pw_pbtio3_berry/` 编写 `INPUT-scf`、`INPUT-nscf-a/b/c` 与 `KPT-scf`、`KPT-nscf-a/b/c`（只写文件，不运行）。


Prompt: 为 Pt(111) 平板准备隐式溶剂化平面波 SCF 输入。读取工作区 STRU，在 `run_pw_pt_slab_imp_sol/` 编写 `INPUT`（`gamma_only` 时可不写 `KPT`；只写文件，不运行）。


Prompt: 为 Pt(111) 平板准备外加电场平面波 SCF 输入。读取工作区 STRU，在 `run_pw_pt_slab_efield/` 编写 `INPUT`（`gamma_only` 时可不写 `KPT`；只写文件，不运行）。


Prompt: 为 Pt(111) 平板准备偶极修正平面波 SCF 输入（`run_pw_pt_dipole_only/`、`run_pw_pt_dipole_efield/`）。读取工作区 STRU，各目录编写 `INPUT`（`gamma_only` 时可不写 `KPT`；只写文件，不运行）。


Prompt: 为 Si 准备 LCAO SCF 静电势输出输入。读取工作区 STRU，在 `run_lcao_si_epot/` 编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为 Pt(111) 平板准备补偿电荷 gate 平面波 SCF 输入。读取工作区 STRU，在 `run_pw_pt_slab_gate/` 编写 `INPUT`（`gamma_only` 时可不写 `KPT`；只写文件，不运行）。


Prompt: 为气相 H₂O counterpoise（BSSE）准备 LCAO 输入。工作区有 STRU_0 ~ STRU_3 四个子计算的结构文件，在 `run_lcao_h2o_cp/` 编写一个共享的 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为 Fe₂ 二聚体准备 LCAO 自旋极化 PEXSI SCF 输入。读取工作区 STRU，在 `run_lcao_fe2_pexsi/` 编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 为 Si₁₆ 准备两套 GPU 加速 SCF 输入（`run_pw_si16_gpu/`、`run_lcao_si16_gpu/`）。读取工作区 STRU，各目录编写 `INPUT` 与 `KPT`（只写文件，不运行）。


Prompt: 请为 Al (FCC, a=4.05 Å) 的能带计算准备 ABACUS KPT 文件。使用 Line 模式，沿高对称路径 Γ → X → W → K → Γ → L，每段 20 个插值点。将 KPT 文件写入当前目录。


Prompt: 工作区中有 NFPP (Na₄Fe₃(PO₄)₂P₂O₇) 的 STRU 文件 `STRU_NFPP`。请准备 ABACUS LCAO cell-relax 计算输入，使用 DFT+U（Fe 3d, U=3.5 eV），在 `nfpp_relax/` 目录下生成 INPUT 和 KPT（只写文件，不运行）。


Prompt: 工作区中有 NFPP (Na₄Fe₃(PO₄)₂P₂O₇) 的 STRU 文件 `STRU_NFPP`。请准备 ABACUS LCAO cell-relax 计算输入，在 `nfpp_relax/` 目录下生成 INPUT 和 KPT（只写文件，不运行）。


Prompt:
工作区中有：
- `POSCAR`: Fe-Sb-N 界面 slab 结构（82 Fe + 1 Sb + 1 N）
- `pp_library/`: 赝势文件目录

请生成 ABACUS PW 基组的自旋极化 SCF 计算输入文件（INPUT、STRU、KPT），用于计算该界面体系的电子结构。

赝势使用 `pp_library/` 目录中已有的文件。


Prompt: 请根据附件里的 Ti-Mo-Al-V SQS 位点表和 ABACUS 可用文件清单，准备一个 `timoyalv_abacus_scf/` 输入目录，用于后续 ABACUS LCAO SCF 计算。至少写出 `INPUT`、`STRU`、`KPT` 和 `input_prep_manifest.json`； 这次只准备输入，不提交计算。


Prompt: 为 NiO 岩盐结构（conventional cell, 8 atoms）准备 ABACUS LCAO DFT+U SCF 输入，计算 AFM-II 基态能量和 band gap。输出到 `run_nio_afm/` 目录。


Prompt: 工作区中有 `Al_fcc.stru`：Al fcc 原胞结构。请准备一个用于 Al fcc 单点 SCF 的 ABACUS plane-wave 输入目录 `run_al_pw/`，至少写出 `run_al_pw/INPUT`、`run_al_pw/STRU` 和 `run_al_pw/KPT`。


Prompt: Create an ARROWS³ search space configuration file `arrows_config.json` for synthesizing YBa₂Cu₃O₆.₅. The available precursors are: Y₂O₃, Y₂(CO₃)₃, BaO, BaCO₃, BaO₂, CuO, CuCO₃, Cu₂O, BaCuO₂, Ba₂Cu₃O₆, Y₂Cu₂O₅. Temperature range: 600°C to 900°C with 100°C intervals. Fixed hold time: 4 h. Calculate N_exp using the formula N_exp = N_sets × ((T_max - T_min)/ΔT + 1) where N_sets is the number of stoichiometrically balanced precursor combinations.


Prompt: Generate a CASTEP .cell input file for computing the elastic constants of L1₀ FeNi. Use the GGA-PBE functional, a 600 eV plane-wave cutoff, a 12×12×12 Monkhorst-Pack k-point grid, electronic SCF tolerance of 1e-6 eV, and a force tolerance of 0.03 eV/Å. Save the file as `l10_feni_elastic.cell` in the current directory.


Prompt:
请为 MgO 岩盐结构（a=4.212 Å）生成 CP2K DFT 单点能计算输入文件，要求：
1. 使用 PBE 泛函，DZVP 基组
2. 平面波截断能 400 Ry
3. SCF 收敛精度 1.0E-6
4. 使用 Goedecker 赝势（GTH-PBE）
输出文件名为 `cp2k.inp`。


Prompt:
工作区里已有 `cp2k_ready.inp`，它是用户准备好的 CP2K 输入。
请对这个 ready-to-run user input 做提交前检查并整理成运行目录：
1. 不要重新生成或覆盖 `cp2k_ready.inp`，只复制到 `run_cp2k_ready/input.inp`。
2. Validate the existing CP2K input before execution，并把 machine-readable diagnostics 保存为 `run_cp2k_ready/diagnosis.json`。
3. 写 `run_cp2k_ready/input_prep_manifest.json`。它必须是 JSON object，至少包含这些顶层字段：software、task、input_dir、generated_files、user_provided_files、diagnostics、auxiliary_files、assumptions、submit_ready、bohrium_command；其中 `diagnostics.file` 必须是 `diagnosis.json`，`submit_ready` 必须是 JSON boolean。
4. `submit_ready` 只能在没有 unhandled errors or blockers 时为 true；本题只准备目录，不实际提交，也不要臆造提交资源默认值。
最终只交付 `run_cp2k_ready/` 目录和简短说明。


Prompt:
工作区里有 `cp2k_needs_fix.inp`。请做提交前诊断，不要提交作业：
1. 检查这个 CP2K input，并把 machine-readable diagnostics 保存为 `diagnosis.json`。
2. Treating parser success as physical validity 是错误做法；请检查 CUTOFF 和 EPS_SCF 是否有物理意义。
3. 如果 diagnosis has unhandled errors or blockers，必须写 `input_prep_manifest.json` 且 `submit_ready=false`。manifest 必须是 JSON object，至少包含这些顶层字段：software、task、input_dir、generated_files、user_provided_files、diagnostics、auxiliary_files、assumptions、submit_ready、bohrium_command；其中 `diagnostics.file` 必须是 `diagnosis.json`，blocked 时 `bohrium_command` 用空字符串。
4. 不要在 blockers 未处理时提交；不要臆造 Bohrium cmd/image/machine。
输出 `diagnosis.json`、`input_prep_manifest.json` 和简短修复说明。


Prompt: 请为水分子 (H2O) 准备一个 CP2K 几何优化计算的输入文件。使用 PBE 泛函、DZVP 基组、GTH 赝势。将完整的 CP2K 输入文件写入 `water_opt.inp`。


Prompt: 为 NiO 岩盐结构（conventional cell, 8 atoms, Fm-3m）准备 CP2K DFT+U SCF 输入文件，用于计算其电子结构（band gap, magnetic moment）。要求精度足以 区分 AFM-II 磁序与 FM 磁序的能量差（~几十 meV/atom）。输出 `nio_scf.inp`。


Prompt: Generate a JSON file `cvm_coefficients.json` containing the complete Cluster Variation Method (CVM) entropy expression definitions for an FCC lattice. Include all three approximations from Tepesch et al. (1998): tetrahedron, tetrahedron-octahedron (T-O), and 13-14-point. For each approximation, list every subcluster with its Kikuchi-Barker coefficient (a_β) and the multiplicity factor m_{αβ} (number of maximal clusters of type α containing subcluster β). Also include a brief description of each maximal cluster's geometry (which lattice sites it comprises).


Prompt: Generate a DIFFaX input file for simulating XRD patterns of T#2-LiCoO2 structures with stacking faults. Use the six orthorhombic slab definitions from Table 3 of the Carlier et al. 2004 paper (a=4.8450 Å, b=2.8097 Å, c=4.9541 Å) with Co and O atomic positions as specified. Include transition probability matrices for simulating O2-type and O6-type stacking faults. Output the file as diffrax_t2_stacking_faults.inp.


Prompt: Generate a DIFFaX input file for simulating XRD patterns of T#2-LiCoO2 derived from the O2 structure by gliding every second CoO2 slab by (1/3, 1/6, 0). Use the two hexagonal slab definitions from Table 1: a=2.8081 Å, c=4.9512 Å; slab 1 with Co at (0,0,1/2), O(1) at (2/3,1/3,0.3), O(2) at (1/3,2/3,0.7); slab 2 with Co at (2/3,1/3,1/2), O(1) at (0,0,0.3), O(2) at (1/3,2/3,0.7). Output the file as diffrax_t2_basic.inp.


Prompt:
这两个 DP-GEN 输入文件现在不能通过官方 dargs/schema 校验，里面各有 5 个故意埋的错误。请在不替换占位路径的前提下做最小修改，只输出修好的 `param_output.json` 和 `machine_output.json`。
输入文件：
- `param_input.json`: `data/IG_dpgen_dargs_debug_002_20260616/param_input.json`
- `machine_input.json`: `data/IG_dpgen_dargs_debug_002_20260616/machine_input.json`


Prompt:
用 DP-GEN 从头开发一个水盒子的力场；DeepMD 训练、LAMMPS exploration、VASP FP，路径先占位。只输出 `param_output.json` 和 `machine_output.json`。
输入文件：
- `project_note.md`: `data/IG_dpgen_dargs_waterbox_001_20260616/project_note.md`


Prompt: Generate a complete fcc cluster definition file for cluster expansion of alloy energetics. Include all pair clusters from 1st through 4th and 6th nearest neighbor, all triplet clusters, and all quadruplet clusters as defined on the fcc lattice with integer coordinates (unit cell a=2). For each cluster, provide: cluster_id (e.g. V_2,1), cluster_type (pair/triplet/quadruplet), site coordinates as lists of [x,y,z] integer triples, multiplicity (number of symmetry-equivalent clusters), and NN shell label for pairs. Save the output as fcc_clusters.json in the working directory.


Prompt: Generate a SIESTA-style input file (`siesta.fdf`) for a spin-polarized LSDA DFT calculation of an armchair graphene nanoribbon (AGNR). Use the following parameters from the paper: energy cutoff of 400 Ry, double-zeta plus polarization (DZP) basis set, force convergence threshold of 16 pN, and 32 k-points along the 1D Brillouin zone. The system is a 1D periodic nanoribbon with vacuum in the other two directions. Save the file as `siesta.fdf` in the current working directory.


Prompt: Generate a GPUMD nep.in training input file for a NEP4 model covering Ag and W elements. Use the hyperparameters from the UNEP-v1 paper: radial cutoff 6 Å, angular cutoff 5 Å, 5 radial functions with 9 basis functions, 80 hidden neurons, batch size 10000, and 1000000 training generations. Save the file as nep.in in the current directory.


Prompt: 请按 EMD 热导率流程给石墨烯写一份可直接运行的 `run.in`。要求包含平衡段和生产段，并在生产段输出热流相关统计。


Prompt: 请生成一份用于石墨烯 HNEMD 热导率与 SHC 分解的 `run.in`，先平衡再生产，生产段要包含 HNEMD 驱动和 SHC 相关设置。


Prompt:
请写一份用于界面热输运 NEMD 的 GPUMD `run.in`：
- 体系已分为 9 个 group（0=左固定壁, 1=热源, 2-6=传热区, 7=热汇, 8=右固定壁）
- 固定两端壁（group 0 和 8）
- 热源 group 1，热汇 group 7，目标温度 300K
- 需要逐 group 温度统计
- 先平衡 200ps，再生产 2ns
- 时间步 1fs
请将文件保存为工作区根目录下的 `run.in`。


Prompt: 请写一份液态硅的 `run.in`：先高温平衡，再在后续阶段包含自扩散、MSD、黏度和 RDF/ADF 统计。


Prompt: 请准备一个主动学习场景的 `run.in`：需要多势函数观察、不确定性触发保存，以及观察者输出设置。


Prompt: 请写一份用于偶极输出自洽检查的 `run.in`，要求包含多势函数、NVE、`dump_dipole` 与位置输出。


Prompt: 请写一份用于极化率输出自洽检查的 `run.in`，要求包含多势函数、NVE、`dump_polarizability` 与位置输出。


Prompt:
请准备一份多势函数观察模式的 `run.in`，要求使用 `dump_observer observe` 输出观测结果。
请将文件保存为工作区根目录下的 `run.in`。


Prompt:
请准备一份多势函数平均模式的 `run.in`，要求使用 `dump_observer average` 输出平均观测结果。
请将文件保存为工作区根目录下的 `run.in`。


Prompt: 请写一份用于 MSD 自洽测试的 `run.in`，要求包含 NVE、`compute_msd` 与基本积分参数。


Prompt:
请为一个蛋白质体系生成 GROMACS 能量最小化的参数文件，要求：
1. 使用最速下降法（steepest descent）
2. 最大步数 50000
3. 能量收敛判据 emtol = 1000.0 kJ/mol/nm
4. 步长 emstep = 0.01 nm
5. 截断方案使用 Verlet，库仑和范德华截断距离 1.0 nm
输出文件名为 `em.mdp`。


Prompt:
请为 GROMACS energy minimization 准备 `run_gromacs_em/`：
- 生成或 adapt input files：`em.mdp`，参数为 steepest descent、nsteps=50000、emtol=1000、emstep=0.01。
- 验证 `em.mdp` 并保存 machine-readable diagnostics 到 `run_gromacs_em/diagnosis.json`。
- 写 `input_prep_manifest.json`。它必须是 JSON object，至少包含这些顶层字段：software、task、input_dir、generated_files、user_provided_files、diagnostics、auxiliary_files、assumptions、submit_ready、bohrium_command；其中 `diagnostics.file` 必须是 `diagnosis.json`，`submit_ready` 必须是 JSON boolean。
- manifest 中必须包含 `submit_ready` gate；如果有 blockers 就不能提交。
本题只生成输入目录，不实际运行 GROMACS。


Prompt:
请为分子 `乙醇`（SMILES: `CCO`）完成 GROMACS 小分子参数化与最小可运行验证。使用 `GAFF` 力场路线，投递任务。
将所有产物（提交脚本、参数化文件、运行日志、验证报告）整理到 `run_ligand_ethanol_gaff/` 目录下。


Prompt:
请为分子 `丙酮`（SMILES: `CC(=O)C`）完成 GROMACS 小分子参数化与最小可运行验证。使用 `GAFF2` 力场路线，投递任务。
将所有产物（提交脚本、参数化文件、运行日志、验证报告）整理到 `run_ligand_acetone_gaff2/` 目录下。


Prompt:
请为分子 `乙酸乙酯`（SMILES: `CCOC(=O)C`）完成 GROMACS 小分子参数化与最小可运行验证。使用 `OPLS-AA` 力场路线，投递任务。
将所有产物（提交脚本、参数化文件、运行日志、验证报告）整理到 `run_ligand_ethyl_acetate_opls/` 目录下。


Prompt:
这些是我前面陆续得到的铅表面接枝萘甲酸体系相关文件。
如果我要继续做这个体系的 GROMACS 模拟，请你帮我判断哪些文件可以复用、还需要补哪些文件，并把需要补齐或整理的输入文件生成到 `run_NA14_gromacs/` 目录下，文件名按下面这样放：
- `Pb_NA14_solvated.pdb` - `Pb_NA.itp` - `topol.top`
已提供的文件：
- `Pb_NA14_solvated.pdb` - `Pb_PNA_merged.itp` - `topol_NA14.top`
要求：
- 不实际运行 GROMACS。 - 不从头重建分子结构；已有体系的盒子信息和坐标信息应尽量保留。 - 最后简短说明哪些文件可以复用，哪些文件是本次整理生成的。


Prompt: 为 Lu4N4 体系生成一个能带结构计算的 VASP INCAR 文件。要求：非自洽计算，从已有电荷密度读取，输出投影能带。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为有机-无机杂化钙钛矿 C26H28Cl4N4Pb2 生成一个 静态自洽(SCF)计算的 VASP INCAR 文件。体系是半导体，需要高精度收敛。 将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为钙钛矿 ErNaO3 生成一个结构弛豫(relaxation)计算的 VASP INCAR。需要同时优化离子位置和晶胞参数。Er 含 f 电子，需要自旋极化。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为金属有机框架(MOF) C30H22N4O10Zn2 生成一个 HSE06 杂化泛函静态计算的 VASP INCAR。需要精确的电子结构和带隙。 将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为含稀土和过渡金属的氧化物 RbPmMo2O8 生成一个 DFT+U 计算的 VASP INCAR。需要对 Mo 的 d 轨道施加 Hubbard U 校正。元素顺序为 Rb, Pm, Mo, O。体系为磁性体系。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为金属间化合物 Tm2ZnGa2 生成一个 SCAN metaGGA 泛函静态计算的 VASP INCAR。Tm 含 f 电子，需要自旋极化。 将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为 Heusler 合金 MnNaPm2 生成一个分子动力学(MD) 计算的 VASP INCAR。Mn 含磁性，需要自旋极化。MD 时间步长设为 0.5 fs。 将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为钙钛矿 BaZrO3 生成一个弹性常数计算的 VASP INCAR。使用有限差分法计算弹性张量，位移幅度 0.015 A。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为有机金属氧化物 C8H12HfO4 生成一个 DFPT(密度泛函微扰理论) 介电常数计算的 VASP INCAR。需要计算介电张量和 Born 有效电荷。 将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为含重元素 Pb 的卤化物 Br6H6Pb2S2 生成一个 自旋轨道耦合(SOC)计算的 VASP INCAR。Pb 的 SOC 效应显著，需要高精度设置。 将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为三元金属间化合物 CdOs3Pt 生成一个 高精度态密度(DOS)计算的 VASP INCAR。体系为磁性金属，需要投影态密度(PDOS)。 将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为普鲁士蓝类似物 C24Fe8N24Na4 生成一个 NEB(nudged elastic band)+HSE06 过渡态搜索的 VASP INCAR。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为三元碳化物表面 (Ta,V)2C (C33O2Ta16V16) 生成一个表面弛豫计算的 VASP INCAR。表面计算需要固定晶胞形状，体系为金属。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为氢钝化的 SiO2 slab 模型 H16O24Si12 生成一个静态计算的 VASP INCAR。体系是绝缘体 slab，需要固定晶胞形状。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为 BiSI 中 I 替位 S 缺陷体系生成一个 HSE06+SOC 缺陷计算的 VASP INCAR。含重元素 Bi 需要自旋轨道耦合，缺陷计算需要固定晶胞形状，使用 HSE06 获取准确缺陷能级。需要离子弛豫。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为稀土硒化物 DySe 生成一个 GW 计算前驱步骤的 SCF VASP INCAR。Dy 含 f 电子(MAGMOM~7)，需要保存高质量波函数供后续 GW 使用。 将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为简单金属 Al (FCC) 生成一个高精度结构优化的 VASP INCAR，作为后续光学计算的前驱步骤。需要高精度设置和完全优化晶胞参数。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为双钙钛矿 Ca4TaOs 生成一个磁性自洽(SCF)计算的 VASP INCAR。体系含 5d 过渡金属(Ta, Os)，需要自旋极化和适当的初始磁矩。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为硼酸铯 B12Cs4O20 生成一个非磁性绝缘体的静态 SCF 计算 VASP INCAR。体系为非磁性，宽带隙绝缘体。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: 为三元硒化物 Na3AlSe 生成一个磁性金属体系的静态 SCF 计算 VASP INCAR。体系为磁性金属，需要投影态密度。将完整的 INCAR 写入当前目录下的 `INCAR` 文件。


Prompt: Generate a LAMMPS input script for a two-phase coexistence MD simulation to predict the melting temperature of Ni at standard pressure. Retrieve the FCC Ni conventional unit cell (Materials Project mp-134) from the structure database and use it to build a 14x14x28 supercell. Initialize the lower half as liquid (heated to ~1.25×Tm_guess) and the upper half as solid. Run with an NPH ensemble at 1.0312 bar, 1 fs timestep, 120 ps equilibration, and NPT pre-equilibration with 100 fs temperature damping and 1000 fs pressure damping for 10 ps. Save the script as tpc_md_ni.in.


Prompt:
请为一个 FCC 铜体系生成 LAMMPS 分子动力学输入文件，要求：
1. 使用 EAM 势（Cu_u3.eam），晶格常数 3.615 Å
2. 3×3×3 超胞
3. NPT 系综，温度 300 K，压力 0 bar
4. 运行 10000 步，时间步长 1 fs
5. 每 100 步输出热力学信息
输出文件名为 `in.lammps`。


Prompt: Generate a LAMMPS input script for a two-phase coexistence MD simulation of Ni at 10 GPa external pressure. Retrieve the FCC Ni conventional unit cell from the structure database and use it to build a 14x14x28 supercell. Set up NPH ensemble with external pressure of 10 GPa, 1 fs timestep, and 120 ps production run. Save the script as tpc_md_ni_10gpa.in.


Prompt: 为含能材料 EAP（乙二胺高氯酸盐，化学式 C₂H₁₂N₂Cl₂O₈，包含 H/C/N/O/Cl 五种元素）写一个 LAMMPS DeePMD 分子动力学输入脚本。模型文件为 `model.pb`（全周期表 118元素type_map）。体系232原子，2500K NVT模拟100ps，时间步0.1fs。输出 `in.lammps`。


Prompt:
工作区中有 `eap_model_notes.txt`，记录了 EAP（乙二胺高氯酸盐）232 原子 DeePMD 高温 MD 的体系和模型信息。请生成 LAMMPS 输入文件，不需要真实运行。

要求：
- 2500 K NVT，100 ps，时间步 0.1 fs；
- DeePMD 模型文件是 `iter23_0.pb`；
- data 文件的原子类型顺序是 H C N O Cl；
- 输出 `in.eap_deepmd` 和 `run.sh`。


Prompt: 工作区中有一个 20 重复单元的聚乙烯链结构 `polyethylene.xyz`。请使用 GAFF 力场为该体系准备 LAMMPS 结构优化输入，输出 `in.lammps` 和 `data.lammps`。


Prompt: 请生成一个 LAMMPS 输入脚本，使用 `fix sgcmc` 研究 BCC Fe-10at%Cu 合金在 800K 下的短程有序(SRO)。要求：4x4x4 BCC 超胞（128原子），EAM 势，每100 MD 步尝试一次 MC swap，运行 50000 步。输出到 `in.sgcmc`。


Prompt:
生成 LAMMPS 输入文件，对 FCC Cu 单晶进行 z 方向单轴拉伸模拟，并追踪拉伸过程中键长和键角的演化。

要求：
- EAM 势
- 300K NVT 平衡后施加应变
- 输出键长分布和键角分布随应变的变化数据
- 输出为 `in.tensile_cu`


Prompt: Generate a complete M3GNet model configuration JSON file named `m3gnet_form_energy_config.json` for training a formation energy (eV/atom) property prediction model on bulk crystals. Use the hyperparameters from the MatGL paper (Table S4) for the Matbench formation energy model. The configuration must be valid JSON with the `@class`, `@module`, `@model_version`, and `init_args` structure that MatGL expects.


Prompt: Generate a complete CHGNet model configuration JSON file named `chgnet_mlip_config.json` for training a universal machine learning interatomic potential (MLIP) on crystal structures. Use the hyperparameters from the MatGL paper (Table S7). The configuration must be valid JSON with the `@class`, `@module`, `@model_version`, and `init_args` structure that MatGL expects.


Prompt:
目录里有一份 NEP 训练数据 `train.xyz`。

请先看数据，再给出可直接训练的 `nep.in`。
要求：
- `type` 与数据中的元素一致；
- 参数不要照抄模板，要和这份数据规模匹配；
- 不要写超出 NEP 合法范围的值。

同时写一份简短说明，说明你怎么定的关键参数。

输出到当前目录：`nep.in`、`nep_input_note.md`。


Prompt: Implement the Nernst equation to compute the effective electrochemical redox potential of the H₅IO₆/HIO₃ redox couple (E° = 1.626 V) as a function of pH, following the methodology described in Wang et al. (2024) for the LiIn(IO₃)₄ aqueous synthesis system. The half-reaction is H₅IO₆ + H⁺ + 2e⁻ → HIO₃ + 3H₂O. Write a Python script that computes E(pH) for pH values from 0 to 14 in steps of 0.1, and save the results to nernst_redox_potential.csv with columns 'pH' and 'E_V'.


Prompt:
工作区中有 `molecules.csv`，包含 5 个小分子的 SMILES。请用 ORCA 量化计算估算这些分子的 Hansen 溶解度参数（δD, δP, δH），生成可直接提交的输入文件。
要求使用 B3LYP/def2-TZVP 级别。


Prompt:
请为乙醇分子生成 ORCA 输入文件，计算其红外光谱（振动频率 + IR 强度）。
要求使用 B3LYP/def2-TZVP 级别。输出为 `ethanol_ir.inp`。


Prompt:
请为苯甲醛(benzaldehyde)生成 ORCA 输入文件，完成以下两步计算：
1. 几何优化（B3LYP/def2-SVP）
2. 基于优化构型的 TD-DFT 计算（PBE0/def2-SV(P)，10 个激发态）

输出两个输入文件 `opt.inp` 和 `tddft.inp`，以及一个 `run.sh` 脚本按顺序执行。


Prompt:
Implement the Pourbaix potential formula from Wang et al. (2024) as a Python function. The Pourbaix potential (equation 3 in the paper) is:
Φ̄ = (1/N_M) × [(G − N_O·μ_H₂O) − RT·ln(10)·(2N_O − N_H)·pH − (2N_O − N_H + Q)·E]
where N_M, N_O, N_H are the number of metal, oxygen, and hydrogen atoms; Q is the charge; G is the molar Gibbs free energy; μ_H₂O is the chemical potential of water; R is the gas constant; T is temperature; pH and E are the electrochemical conditions.
Write a Python function `compute_pourbaix_potential(formula, G_eV_per_formula, pH, E_V, T=298.15)` that parses a chemical formula (e.g., 'LiFePO4'), counts atoms, and returns the Pourbaix potential in eV/atom. Save the function and a test script to pourbaix_potential.py. The test script should compute Φ̄ for LiFePO₄ (G ≈ −18.0 eV/formula, approximate) at pH=9.27, E=−0.56 V, and for H₂O (G = 0 by convention, Q=0) at pH=7, E=0 V.


Prompt: Generate a Pt(111) surface slab structure using pymatgen. Start from bulk platinum (FCC), create a slab with the (111) Miller index, add a 10 Å vacuum layer, and ensure the slab has symmetric (identical) top and bottom terminations. Save the resulting structure to pt_111_slab.cif in the current working directory.


Prompt: Generate a Si(111) surface slab from bulk silicon using pymatgen. Use the primitive bulk unit cell, create a slab with the (111) Miller index oriented to expose the non-polar termination (no dangling Si atoms), and add a 10 Å vacuum layer. Save the structure to si_111_slab.cif in the current working directory.


Prompt: 请用 PySCF 为水分子（H 0 0 0; H 0 0 0.96; O 居中，可用合理坐标或直接用 `O 0 0 0; H 0 0 0.96; H 0.93 0 -0.24`）写一个单点 DFT 计算脚本 `run_pyscf.py`，使用 wB97X-D 泛函、def2-TZVP 基组，打印总能量。


Prompt: 请用 PySCF 为甲烷分子 CH4 写一个单点能量计算脚本 `run_pyscf.py`，使用 r2SCAN 泛函（meta-GGA）和 def2-SVP 基组，打印总能量。分子坐标可自行给出合理值。只需写出脚本本身，无需提交或运行。


Prompt: 请用 PySCF 为乙烯分子 C2H4 写一个单点 DFT 脚本 `run_pyscf.py`，使用 M06-2X 泛函、cc-pVDZ 基组，打印总能量。坐标可自行给出合理值。只需写出脚本本身，无需提交或运行。


Prompt: 请用 PySCF 为苯分子 C6H6 写一个单点能量脚本 `run_pyscf.py`。方法上请选用一个 带经验色散校正的范围分离杂化泛函（range-separated hybrid + dispersion），基组用 def2-SVP。


Prompt: 请用 PySCF 为甲醇分子 CH3OH 写一个单点 DFT 脚本 `run_pyscf.py`，使用 wB97X-D3 泛函、def2-TZVP 基组，打印总能量。


Prompt: 请用 PySCF 为甲基自由基 CH3（doublet，单重态以外的开壳层体系）写一个单点 DFT 脚本 `run_pyscf.py`，使用 PBE 泛函、def2-SVP 基组。注意正确设置自旋。


Prompt: 请用 PySCF 为水分子写一个几何优化脚本 `run_pyscf.py`，在 PBE/def2-SVP 级别下用 geomeTRIC 优化器优化结构，并打印优化后的总能量。


Prompt: 请用 PySCF 准备一个苯分子单点 B3LYP/def2-SVP 计算：写出 `run_pyscf.py`， 并把可直接执行的 Bohrium 提交命令（含正确的镜像和 cmd）写入 `submit_cmd.txt`。


Prompt: 请用 PySCF 为甲醛分子 CH2O 写一个脚本 `run_pyscf.py`，在 B3LYP/def2-SVP 级别做单点计算，并输出 HOMO-LUMO gap 和偶极矩。只需写出脚本本身，无需提交或运行。


Prompt: 请用 PySCF 为 N2 分子写一个 CCSD(T) 单点能量脚本 `run_pyscf.py`，使用 cc-pVDZ 基组，打印 CCSD(T) 总能量。


Prompt: Generate a Quantum ESPRESSO pw.x SCF input file for diamond (primitive cell) using the parameters from Mounet & Marzari (2005): GGA-PBE functional, 40 Ry wavefunction cutoff, appropriate charge-density cutoff ratio, and an 8x8x8 Monkhorst-Pack k-mesh. Save the input as diamond_scf.in in the current directory.


Prompt: Generate a Quantum ESPRESSO pw.x input file (si_relax.in) for geometry relaxation of silicon (diamond-cubic, Materials Project mp-149), which you first retrieve from the structure database. Use PBE functional and precise-level convergence settings. Relax atomic positions only, keeping the cell shape and volume fixed.


Prompt: Generate a Quantum ESPRESSO ph.x input file for a DFPT phonon calculation of diamond on an 8x8x8 q-point grid. Apply the acoustic sum rule. Use consistent parameters with the SCF run: GGA-PBE, 40 Ry wavefunction cutoff, 8x8x8 k-mesh. Save the input as diamond_phonon.in in the current directory.


Prompt: Generate a SIESTA .fdf input file for a spin-polarized DFT geometry optimization of antiferromagnetic NiO. Retrieve the FCC NiO primitive cell (Fm-3m) from the structure database and use it to build the AFM-II supercell with propagation vector Q=(0.5,0.5,0.5). Apply GGA-PBE, 400 Ry mesh cutoff, 13x13x13 k-point grid, double-zeta polarized basis, and fully relativistic psml pseudopotentials. Converge forces to 0.01 eV/Ang. Write the output to nio_siesta.fdf.


Prompt: Retrieve the bulk silicon structure (diamond cubic) from the structure database, then generate two VASP INCAR files: INCAR_relax for PBE structural relaxation and INCAR_hse for one-shot HSE06 band gap calculation. Follow the AMP² computational settings described in Kim et al. (Scientific Data, 2020): PBE functional for relaxation, HSE06 (HFSCREEN=0.2, AEXAC=0.25) for the one-shot calculation. Place both files in the current directory.


Prompt: Generate VASP input files (INCAR, POSCAR, KPOINTS) for geometry relaxation of GeTe (trigonal R-3m, Materials Project mp-1189580), which you first retrieve from the structure database. Optimize atomic positions and cell shape while keeping the cell volume fixed. Use PBE functional and precise-level settings. No spin polarization.


Prompt: Generate a VASP INCAR file for structural optimization of bulk pyrite FeS2 (cubic Pa-3), which you first retrieve from the structure database. Use the PBE functional with DFT+U correction (U=2.0 eV on Fe-3d). Save the INCAR to the working directory.


Prompt: Generate a VASP INCAR file for a two-stage geometry optimization of a 128-atom BCC solid solution alloy. Use GGA-PBE functional, 350 eV plane-wave cutoff, 2×2×2 k-point mesh, Methfessel-Paxton smearing with 0.2 eV broadening. The first stage uses PREC=Low and the second stage uses PREC=Normal. Set convergence criteria to 10⁻⁴ eV for energy and 10⁻² eV/Å for forces. Save the INCAR file as `INCAR` in the current working directory.


Prompt: Generate a VASP INCAR file for a spin-polarized GGA+U calculation of NiO. Use Dudarev's simplified U-J approach (LDAUTYPE=2) with U=6.4 eV on Ni 3d orbitals, a 550 eV plane-wave cutoff, PAW pseudopotentials, and settings appropriate for full structural relaxation (ionic and cell degrees of freedom). Include reasonable convergence criteria for an insulating transition-metal oxide. Write the result to incar_nio_ggau.vasp in the current working directory.


Prompt: Generate a VASP INCAR file for PBEsol geometry optimization of the primitive CsPbI₃ perovskite structure, which you first retrieve from the structure database (e.g. Materials Project). The INCAR must use the PBEsol functional, a 500 eV plane-wave cutoff, Γ-centered k-point sampling, and convergence criteria of 1e-5 eV for energy and 0.01 eV/Ang for forces. Save the INCAR to the current working directory.


Prompt: Generate a VASP INCAR file for structural optimization of bulk anatase TiO2 (tetragonal I41/amd), which you first retrieve from the structure database. Use the PBE functional with DFT+U correction (U=7.3 eV on Ti-3d). Save the INCAR to the working directory.


Prompt: Generate a VASP INCAR file for DFT single-point energy and force calculations using the settings from the UNEP-v1 training data generation: PBE functional, 600 eV plane-wave energy cutoff, PREC=Accurate, Gaussian smearing with SIGMA=0.02 eV, electronic convergence threshold EDIFF=1E-6, blocked Davidson algorithm (ALGO=Normal), and ISMEAR=0. Save as INCAR in the current directory.


Prompt: Generate a VASP INCAR file for a spin-polarized GGA+U calculation of MnO (rock-salt, AFM). Use Dudarev's simplified U-J approach (LDAUTYPE=2) with U=4.0 eV on Mn 3d orbitals, a 550 eV plane-wave cutoff, PAW pseudopotentials, and settings appropriate for full structural relaxation. Include reasonable convergence criteria for an insulating transition-metal oxide and initial magnetic moments consistent with the AFM ground state. Write the result to incar_mno_ggau.vasp in the current working directory.


Prompt: Generate VASP input files (INCAR, POSCAR, KPOINTS) for a SCAN meta-GGA calculation of an isolated O2 dimer in a 15×15×15 Å³ cubic box. Use the parameters from Zhang et al. (npj Comput. Mater. 2018): plane-wave cutoff 520 eV, PAW potentials, gamma-point only k-mesh (appropriate for a large molecular box), convergence to 10⁻⁷ eV in total energy and 0.01 eV/Å on forces. Place the O2 bond along the z-axis with a bond length of 1.21 Å. Save all three files to the working directory.


Prompt: 工作区中有 NiO.vasp（氧化镍 POSCAR）。请生成静态计算 INCAR 文件 `INCAR`。PBE 泛函，ENCUT=550 eV。


Prompt:
工作区中有 `fenc_restart_bundle.json`，记录了 FeN4/graphene ORR 中间体已有计算目录、上一版公共 INCAR 和可用的 CONTCAR 情况。请据此生成下一轮“稳定重跑”会用到的 VASP 输入文件，不需要真实提交计算。

请输出 `INCAR_restart`，并写一个 `restart_inputs/check_before_submit.sh`，用于在提交前检查 POTCAR、CONTCAR/OSZICAR 和 VASP 启动命令等基本条件。


Prompt:
我想用 VASP 做 Ni 中辐照缺陷的电子结构计算：完美 Ni、一个 He4V4 氦泡模型、一个 V6 空洞模型。工作区中有 `potcar_inventory.json`，请根据这个 POTCAR 清单，在 `vasp_inputs/` 下为能准备的体系生成输入模板。

请把需要用户后续确认或补齐的信息统一写在 `vasp_inputs/README.md`。

不需要真实提交计算。


Prompt:
工作区中有 Fe_bcc.vasp（体心立方铁，2个原子）。
请生成适合该体系的 SCF 计算 INCAR 文件 `INCAR_mag`。
注意：Fe 是典型的铁磁性金属，请确保计算设置正确反映其磁性。
ENCUT = 520 eV，电子步收敛精度 1E-6 eV，PBE 泛函。


Prompt:
工作区中有 TiO2_rutile.vasp（金红石结构 TiO2，6 个原子）。
生成 r2SCAN 泛函结构优化的 INCAR 文件 `INCAR_r2scan`，要求：
- r2SCAN meta-GGA 泛函
- 结构优化（离子+晶胞）
- ENCUT = 600 eV
- 力收敛标准 0.01 eV/A


Prompt:
工作区中有 `qha_systems.json`，列出了 4 个晶体结构的原子数、元素和 QHA 体积点。请在 `qha_inputs/` 下生成 VASP + PBEsol 做 QHA 会用到的输入模板和命令脚本，不需要真实提交计算。

需要至少包含：
- `qha_inputs/dfpt/INCAR`
- `qha_inputs/finite_displacement/INCAR`
- `qha_inputs/phonopy_commands.sh`


Prompt: 请为 Si 的 ABACUS LCAO-SCF 试算准备当前目录下的 `INPUT`、`STRU`、`KPT`，并把 Bohrium 提交配置写入 `submit_job.json`。


Prompt: 请为 Si 的 ABACUS PW-SCF 试算准备当前目录下的 `INPUT`、`STRU`、`KPT`，并把 Bohrium 提交配置写入 `submit_job.json`。


Prompt: 工作区中有 Al₂O₃ 的 CIF 文件 `Al2O3.cif`（刚玉结构，R-3c，10原子）。 请用 ABACUS 对其做 SCF 计算并输出电荷密度，我后续要做 Bader 电荷分析看 Al-O 之间的电荷转移。 帮我把输入文件放在 `abacus_bader/` 目录下并提交运行。


Prompt:
工作区中有一个 Pd(111) slab 的 CIF 文件 `Pd111_slab.cif`（16 个 Pd 原子，4 层）。
请将其转换为 ABACUS LCAO 格式并准备好 SCF 计算输入，我要直接提交跑。


Prompt: 工作区中有 Fe₂O₃（hematite，R-3c）的 CIF 文件 `Fe2O3.cif`（10原子原胞）。 请用 CP2K 对其做 PBE+U (U_eff=4.0 eV on Fe 3d) 的 SCF 单点能计算。 使用 DZVP-MOLOPT-SR-GTH 基组、截断能 600 Ry、k 点 Gamma only。 帮我准备好输入并跑一下。


Prompt: 工作区中有 Si 金刚石结构的 CIF 文件 `Si.cif`（2 原子原胞）。 请用 ABACUS 计算 Si 的声子色散曲线（有限位移法），沿 G-X-K-G 路径绘制， 并报告 Gamma 点光学支频率（THz）。


Prompt: 工作区 `data/` 下是我的一份 DeePMD 格式训练数据集。帮我用 DPA3 预训练模型在这份数据上微调大约 3000 步，跑完把训练和测试（验证）的 loss 告诉我。


Prompt: 请使用 DPA-3.2-5M 模型对 MgO 做结构优化，输出提交配置 `opt_submit.json`。


Prompt: 请使用 DPA 模型对咖啡因分子做构型优化，输出为 `caffeine_optimized.cif`。


Prompt: 对 Al₂O₃ 准备 ABACUS Bader 电荷分析所需文件和 Bohrium 提交配置（文件名 `al2o3_bader_submit.json`）。本题只检验输入文件与提交 JSON，不要投递或执行计算。


Prompt:
使用 GPUMD 计算 Si 金刚石结构的声子色散关系。

要求：
- 使用 NEP 势函数（Si 体系，从 GPUMD 官方 examples 或已有资源获取）
- 提交到 Bohrium 运行
- 报告 Gamma 点光学支频率
- 结果写入工作区根目录 `phonon_result.json`：`{"gamma_optical_THz": ..., "method": "gpumd"}`


Prompt: 请用 ORCA 计算萘 (naphthalene, C₁₀H₈) 的第一激发态吸收波长。 先做几何优化（B3LYP/def2-SVP），然后在优化结构上做 TD-DFT（PBE0/def2-SVP，5 个根）。 报告 S₁ 的激发能 (eV) 和振子强度。


Prompt:
使用 Quantum ESPRESSO 计算 Si 的能带结构（金刚石结构，a=5.43 Å）。
沿 Γ-X-L-Γ 高对称路径计算能带，报告带隙类型和带隙值。


Prompt: 请帮我对 MgO 岩盐结构做一次 Quantum ESPRESSO pw.x 的 SCF 单点能计算。 结构参数：Fm-3m, a=4.21 Å, Mg 在 (0,0,0), O 在 (0.5,0.5,0.5)。 使用 PBE 泛函、ecutwfc=60 Ry、4×4×4 k 点网格。帮我跑一下看看总能是多少。


Prompt: 帮我用 VASP 计算一个 O₂ 分子的平衡键长和结合能，用 PBE 泛函。

VASP PAW datasets may not be redistributed with the public benchmark.

Obtain PAW_PBE/O/POTCAR through a valid VASP license and provision it to the runner under this exact logical path.


Prompt:
工作区中提供了 Si 金刚石结构 `POSCAR_Si` 和对应的 `POTCAR_Si`。
请准备 VASP 自洽场计算输入，放入 `si_scf/` 目录（含 POTCAR），并输出提交配置 `submit.json`。

Obtain the standard PAW_PBE Si POTCAR through a valid VASP license and provision it to the runner under this exact logical path.


Prompt:
文件 `pxrd_303K.xy` 是一个有机-无机杂化硝普盐在 303 K（30 °C）下采集的粉末 X 射线衍射（PXRD）数据（Cu Kα1, λ = 1.5406 Å, 2θ 范围 5°–50°, 步长 ~0.02°, 两列格式: 2θ 强度）。
该化合物为单斜晶系，空间群 P 21（No. 4）。文献参考晶胞：a ≈ 10.83 Å, b ≈ 9.62 Å, c ≈ 10.13 Å, beta ≈ 108.75°, V ≈ 1000 Å³。
精修策略建议： - 建议使用 14°–50° 2θ 区间以避开低角伪影 - 基线 / 背景函数使用低阶（0–1 阶）以避免过拟合 - 若多次精修收敛到不同盆地，选择 wR 最低且体积与初猜偏差最小的结果
请对该 PXRD 图谱进行 Pawley 精修，提取精修后的晶格参数并报告精修质量。
【机器评分约定】最终答案末尾追加一个 eval_results 标签包裹的 JSON 块（数值无单位，约定： a/b/c 单位 Å, beta 单位°, V 单位 Å³, wR 单位 %），格式如下：
<eval_results> {"a": ..., "b": ..., "c": ..., "beta": ..., "V": ..., "wR": ...} </eval_results>
数值必须来自实际 Pawley 精修结果，不得填入 prompt 中的初猜近似值。


Prompt:
文件 `_td-pxrd_test.txt` 包含一个化合物在 8 个温度下的变温粉末 X 射线衍射（PXRD）数据，温度标签对应关系为：30°C→303K、50°C→323K、70°C→343K、90°C→363K、110°C→383K、120°C→393K、130°C→403K、140°C→413K。该化合物为单斜晶系，空间群提示为 P 21（No. 4），波长为 Cu Kα1，λ = 1.5406 Å。
数据特征与精修策略： - 数据为实测变温实验数据，扣除背景前的动态范围约 10×–30×；2θ 范围 5°–50°，建议精修使用 14°–50° 区间以避开低角伪影。 - 推荐将数据按相拆分（RTP 4 点、HTP 4 点）分别精修，并在同一相内启用 chain-cell（前一温度精修结果作为下一温度初猜）。 - 同相内的体积/晶胞参数应单调或近似单调随温度变化；若某温度精修结果 V 相对邻近温度跳变 >5%，应视为收敛到错误盆地（可能为 √2 或 2× 的伪解），需要从相邻温度结果重新初始化或多起点重精修。
运行预算（重要）：本题目硬性 wall-clock 上限 1200 s（线上 task timeout）。Pawley 精修请按以下设置以兼顾稳定性与时限： - 使用 `--multi-start 1`（每相单条 chain 即可；niggli 标准化 + 双向 chain-cell 已能覆盖绝大多数错盆地，不需要 5×/8× 起点的暴力多起点；脚本默认开启的 `--self-heal-chain` 会在 chain merge 之后对 V 跳变 >2% 的单点自动用 multi-start 5 重精修，足以兜住偶发塌方）； - 使用 `--chain-cell --chain-cell-direction both`（脚本内部会做正反向 merge 与体积一致性裁决）； - 使用 `--standardize-cell ref`（基于参考晶胞做轴置换 + niggli 还原，规避单斜/正交等晶系的等效设定切换）； - 必须**先并行 submit RTP 与 HTP 两个 Bohrium 任务**（两次 `Bohrium(action='submit')` 一次性发出，**不要先 submit RTP → 等结果 → 再 submit HTP**），随后再循环 `poll` 检查；串行 submit 会因为单作业 8–12 min 直接超时。 - **拿到 results.json 后不要再发起新的 Bohrium 精修作业**：脚本内置的 `chain-cell` + `standardize-cell` + `self-heal-chain` 已经把邻居体积一致性与错盆地单点重精修做完了，agent 端只需要读 `results` / `merge_audit` / `self_heal_audit`，做最终的相分组与线性拟合，直接写答案。重新提交 Bohrium 作业 1 次就会吃掉 8–12 min，几乎必然超时。
文献参考晶胞（每相分别作为 Pawley 精修初猜，并基于此自由精修；c 与 β 是单斜结构敏感参数，如下文献值已对每相核对）： - 低温相 RTP（303K – 363K）：a ≈ 10.83 Å, b ≈ 9.62 Å, c ≈ 10.13 Å, β ≈ 108.75°, V ≈ 1000 Å³ - 高温相 HTP（383K – 413K）：a ≈ 11.05 Å, b ≈ 9.48 Å, c ≈ 10.48 Å, β ≈ 111.25°, V ≈ 1023 Å³
请对每个温度的 PXRD 图谱进行 Pawley 精修，提取晶格参数（a、b、c、β、V），并分析结果： 1. 识别两个结构相：低温相（1-RTP，303K–363K）和高温相（1-HTP，383K–413K）。 2. 分别对两个相的各晶格参数（a、b、c、β、V）随温度的变化进行线性拟合，给出拟合方程（斜率、截距、R²）。两相必须分别独立拟合，不能合并为一条直线，且 RTP 与 HTP 的斜率应基于各相内的 4 个温度点分别计算（不要让两相共享同一斜率值）。 3. 报告各温度下的精修结果。每个温度必须使用 Pawley 脚本（如 `gsas2_pawley.py`）输出 JSON 中真实的 volume（由收敛后的 a/b/c/β 计算得到），不要将 prompt 中的近似初猜（V ≈ 1000、1023 Å³）当作精修结果填入答案。
【机器评分约定】最终答案末尾追加一个 eval_results 标签包裹的 JSON 块（数值无单位，约定单位为 Å / Å³ / Å·K⁻¹ / Å³·K⁻¹），至少需含以下 dot-path（其余字段如其它温度、b/beta、intercept、r_squared 可按情况补充；fits 为同相内 4 点线性拟合，RTP 与 HTP 必须独立）：
- rtp.{303K,363K}.V，rtp.303K.a，rtp.fits.V.slope，rtp.fits.a.slope - htp.{383K,413K}.V，htp.383K.a，htp.fits.V.slope，htp.fits.c.slope
数值必须来自实际 Pawley 精修与线性拟合，不得填入 prompt 中的初猜近似值。
最终答案中请额外用 2–4 句简述审计信息： - 精修质量审计：`curation.verdict` / `warnings` 是否提示问题；`merge_audit.warnings` 或 `self_heal_audit.outliers` 是否为空，如不为空请指出受影响的温度点。 - Bohrium 执行审计：总共发起了几次 `Bohrium(action='submit')`（本题应 ≤3）；提交命令是否使用 `cmd='bash run.sh > log 2>&1'`；下载后是否确认每个结果目录包含 `log` 与 `results.json`；结果是否来自 Bohrium 上运行的 `gsas2_pawley.py` 而不是本地 scipy/numpy 近似拟合。


Prompt:
从 COF 数据库检索 3D 类型、dia 拓扑的 COF 结构。

请完成以下交付：
1) `3d_dia_cof.json`：包含 `query_type`, `query_topology`, `candidate_count`, `results` （每个含 name, topology, type, formula, density, void_fraction）；
2) `3d_dia_report.md`：说明 3D COF 与 2D COF 在结构上的核心区别、dia 拓扑的连接规则（四面体节点+线性连接体），以及检索到的结构是否符合预期；
3) 导出第一个候选结构为 `COF_3D_dia.cif`。


Prompt:
分别检索 hcb 拓扑和 sql 拓扑的 COF 各 3 个，对比两种拓扑在孔径(PLD/LCD)和比表面积上的差异。

请完成以下交付：
1) `topology_comparison.json`：包含 `hcb_results`（3个）和 `sql_results`（3个），每个含 name, topology, pld, lcd, surface_area_m2g；
2) `topology_comparison_report.md`：对比分析两组数据的孔径和比表面积分布差异，解释为什么不同拓扑会导致不同的孔道几何。


Prompt:
检索同时含有硼(B)和氮(N)元素的 COF 结构，返回最多 5 个候选。

请完成以下交付：
1) `bn_cof_candidates.json`：包含 `query_elements`, `candidate_count`, `results`（每个含 name, formula, elements, topology, surface_area）；
2) `bn_cof_report.md`：分析这些 COF 的共同特征（连接化学类型、典型拓扑、孔径范围），解释为什么 B+N 组合在 COF 中常见（硼酸酯键、硼嗪键等）。


Prompt:
检索 COF-LZU1 的晶体结构，并报告其基本结构信息。

请完成以下交付：
1) `COF_LZU1.cif`：导出结构文件；
2) `cof_lzu1_info.json`：包含 `name`, `topology`, `elements`, `formula`, `space_group`, `pld`, `lcd`, `surface_area`, `type`(2D/3D), `database_source`；
3) 在报告中确认这是否确实是亚胺键连接的 hcb 拓扑 COF。


Prompt:
检索孔限径(PLD)大于 20 埃的 COF 结构，按孔径从大到小排列，返回前 5 个。

请完成以下交付：
1) `large_pore_cof.json`：包含 `filter_criterion`, `results`（每个含 name, pld, lcd, topology, surface_area, type）；
2) `large_pore_cof_report.md`：分析大孔 COF 的结构特征（什么样的拓扑/连接体长度容易产生大孔）；
3) 导出孔径最大的结构为 `largest_pore_cof.cif`。


Prompt:
帮我筛选比表面积大于 2000 m²/g 的 COF 材料，我想找适合甲烷储存的高孔隙率结构。

请完成以下交付：
1) `high_sa_cof.json`：包含 `filter_criterion`, `application`, `results`（每个含 name, surface_area_m2g, void_fraction, pld, lcd, topology）；
2) `gas_storage_screening.md`：从检索结果中推荐 2-3 个最适合甲烷储存的候选，说明推荐理由 （孔径匹配甲烷动力学直径~3.8Å、高空隙率、合适拓扑）。


Prompt:
帮我从 COF 数据库中检索 sql 拓扑、含三嗪(triazine)节点的 COF 结构。

请完成以下交付：
1) `cof_retrieval_report.md`：说明检索条件、使用了哪个数据库、返回了多少候选、各候选的关键属性（拓扑、元素、孔径等）；
2) `cof_candidates.json`：至少包含 `query_topology`, `database_source`, `candidate_count`, `results`（每个结果含 name, topology, elements, pld, lcd）；
3) 将筛选出的第一个结构导出为 `COF_sql_triazine.cif`。


Prompt: 工作区中有 ABACUS SCF 计算结果文件 `abacus_results.json`，包含 Pt(111) clean slab、CO 分子、以及 CO 吸附在 Pt(111) 上的三个体系的总能量。请计算 CO 的吸附能 E_ads = E(CO+slab) - E(slab) - E(CO)，输出到 `adsorption_energy.json`，格式为 `{"E_ads": ...}`。


Prompt:
工作区中有一个已完成的 ABACUS slab 计算结果：
- `mock_running_scf.log`：SCF 收敛日志（含费米能级）
- `mock_ElecStaticPot_z_avg.dat`：沿 z 方向平均的静电势（两列：z 分数坐标 + 势能值）

请从这些数据中提取功函数（work function），写入 `work_function.json`，格式为 `{"work_function_eV": <数值>, "V_vacuum_eV": <数值>, "E_fermi_eV": <数值>}`。


Prompt:
工作区中有一个弛豫后的 H/Zn(0001) 吸附结构 `H_on_Zn0001_relaxed.xyz`（2x2 超胞，4 层 slab + 1 个 H 原子）。请分析该 H 原子处于什么吸附位点，并给出判断依据。
输出 `site_analysis.json`，格式：`{"actual_site": "...", "coordination": ..., "distances_to_nearest_Zn_A": [...], "reasoning": "..."}`


Prompt: Calculate the Debye temperature Θ_D (in K) of Ni₃Al using the Moruzzi-Janak-Schwarz approximation. The bulk modulus is B = 180 GPa, the average atomic volume is Ω = 11.14 Å³/atom, and the average atomic mass is M̄ = 45.24 amu. Write your answer to a file named debye_temperature.txt in the working directory, with the numeric value on a single line.


Prompt: Calculate the decomposition enthalpy (ΔH_decomp) in meV/atom for CsGe₀.₅Sn₀.₅Br₃. Use the formula ΔH_decomp = E(CsGe₀.₅Sn₀.₅Br₃) - E(CsBr) - 0.5*E(GeBr₂) - 0.5*E(SnBr₂). Fetch total energies per atom from Materials Project for CsBr (mp-1066996), GeBr₂ (use OQMD or estimate), SnBr₂ (mp-23215), and construct the perovskite energy from CsPbBr₃ (mp-570587) by substituting Ge/Sn for Pb. Report ΔH_decomp in meV/atom to a file decomp_energy.json with keys: E_perovskite_eV, E_CsBr_eV, E_GeBr2_eV, E_SnBr2_eV, delta_H_decomp_meV_per_atom.


Prompt: Calculate the nearest-neighbor effective cluster interaction V_1nn (in meV/atom) for the Pd-V alloy system at T = 1000 K using the bond proportion model. The spring constants are: k_PdPd = 3.2 eV/Å², k_VV = 2.8 eV/Å², and k_PdV = 3.5 eV/Å². The system is fcc (d = 3). Write your answer to a file named v1nn_meV.txt in the working directory, with the numeric value on a single line.


Prompt: Calculate the mixing entropy ΔS_mix (in meV/atom/K) and the free-energy-corrected decomposition enthalpy ΔH_decomp - TΔS_mix (in meV/atom) at T = 298 K for a B-site alloyed perovskite with composition CsGe₀.₂₅Hg₀.₂₅Pb₀.₂₅Sn₀.₂₅Br₃. Use the ideal solid solution model: ΔS_mix = -k_B Σ x_i ln(x_i), where k_B = 8.617333262×10⁻⁵ eV/K. Given ΔH_decomp = -40 meV/atom, compute -TΔS_mix and ΔH_decomp - TΔS_mix. Write results to entropy_calc.json with keys: composition, delta_S_mix_eV_per_atom_K, T_K, minus_T_delta_S_mix_meV_per_atom, delta_H_decomp_meV_per_atom, delta_H_decomp_minus_T_delta_S_mix_meV_per_atom.


Prompt: Calculate the volumetric thermal expansion ΔV/N (in Å³/atom) at T = 1000 K for an fcc alloy using the quasiharmonic approximation. The bulk modulus is B = 150 GPa and the average Grüneisen parameter is γ̄ = 2.0. Write your answer to a file named thermal_expansion.txt in the working directory, with the numeric value on a single line.


Prompt: Calculate Bartel's tolerance factor τ for CsGe₀.₅₆₂₅Sn₀.₄₃₇₅Br₃. Use Shannon ionic radii (6-coordination, in Å): Cs⁺ = 1.81, Ge²⁺ = 0.73, Sn²⁺ = 1.18, Br⁻ = 1.96. n_A = 1 (Cs oxidation state). r_B is the composition-weighted average. Use the formula: τ = r_X/r_B - n_A*(n_A - (r_A/r_B)/ln(r_A/r_B)). Write the result to tolerance_factor.json with keys: formula, r_A, r_B, r_X, n_A, tau, is_perovskite_formable (true if τ < 4.18).


Prompt: The purely configurational order-disorder transition temperature of Ni3Al (L1₂ → fcc disordered) is T_config = 1350 K. The configurational entropy change upon disordering at 50:50 composition is ΔS_config = 0.693 k_B/atom. The vibrational entropy change upon disordering is ΔS_vib = 0.05 k_B/atom. Compute the corrected transition temperature T_config+vib including vibrational effects. Write your answer to a file named corrected_transition_temp.txt in the working directory, with the numeric value in Kelvin on a single line.


Prompt: From the paper 'Factors that affect Li mobility in layered lithium transition metal oxides' by Kang and Ceder, extract the activation barrier data for Li migration in LiCoO₂ as a function of Li slab distance (Figure 2). Create a CSV file named licoo2_barrier_vs_slab.csv with columns: slab_distance_angstrom and activation_barrier_meV. Also compute and report the interpolated activation barrier at the equilibrium slab distance of 2.64 Å.


Prompt: 工作区中有 ABACUS 计算结果文件 `abacus_energies.json`，包含 LiCoO2 体相、Li 金属、Co 金属和 O2 分子的总能量。请计算 LiCoO2 的生成能（每 formula unit），输出到 `formation_energy.json`，格式为 `{"E_f_per_fu": ...}`。


Prompt: 请计算 FEC（氟代碳酸乙烯酯，SMILES: C1OC(=O)OC1F）的电化学氧化电位 vs Li/Li+，输出到 `fec_oxidation.json`，格式为 `{"E_ox_vs_Li": ...}`。


Prompt: Calculate the theoretical Mn2+/Mn4+ double redox capacity (in mAh/g) of the cathode material Li2Mn2/3Nb1/3O2F. The Mn2+/Mn4+ couple involves 2 electrons per Mn ion. Use atomic masses: Li=6.94, Mn=54.94, Nb=92.91, O=16.00, F=19.00 g/mol. Save your calculation steps and result to a file named capacity_calculation.txt.


Prompt: Compute the theoretical gravimetric capacity in mAh/g for the composition Li1.171Mn0.343V0.486O1.8F0.2, assuming the Mn²⁺/⁴⁺ couple (2 electrons per Mn) and V⁴⁺/⁵⁺ couple (1 electron per V) are fully accessible. Calculate both the transition-metal-based capacity and the Li-based capacity, and identify which is the limiting factor. Save the results as a JSON file named capacity_analysis.json with keys: tm_capacity_mAh_g, li_capacity_mAh_g, limiting_factor, and composition.


Prompt: Construct the domain matrix D for a binary fcc cluster expansion using Eq. (29) from Mueller & Ceder (2010), which weights all compositions uniformly. The cluster expansion includes 9 orbits with site counts n = [0, 1, 2, 2, 3, 3, 4, 4, 4] (empty, single-site, four 2-site pairs, two 3-site triplets, one 4-site quadruplet). Use D[i][j] = 1/(n[i]+n[j]+1) when (n[i]+n[j]) is even, and 0 when odd. Save the 9×9 domain matrix as `domain_matrix_eq29.json` with key "D" (list-of-lists). Also include flat keys "trace" (sum of diagonal), "D_00", "D_11", "D_22", "D_01", "D_12", "D_23" with the corresponding matrix element values.


Prompt: Compare four domain matrix construction methods for a binary fcc cluster expansion with 9 orbits (site counts n = [0,1,2,2,3,3,4,4,4]): (1) van de Walle & Ceder: D = 9×9 identity matrix. (2) Seko et al.: diagonal matrix with D[α][α] = 1/N_α where N = [1,1,4,4,8,8,12,12,12] (symmetrically distinct clusters per orbit in a 32-atom fcc supercell), off-diagonal = 0. (3) All structures equal (Eq. 24): D[0][0]=1, all other elements 0. (4) Uniform composition weighting (Eq. 29): D[i][j] = 1/(n[i]+n[j]+1) if (n[i]+n[j]) is even, else 0. For each method, compute the 9×9 matrix and its trace (sum of diagonal). Save all four matrices and their traces in `domain_matrices_comparison.json` with keys "vandewalle", "seko", "all_structures", "composition_weighting" for matrices and "vandewalle_trace", "seko_trace", "all_structures_trace", "composition_trace" for traces.


Prompt:
我之前写论文时引用了下面这句话里的文献，但现在怎么都查不到这篇：

> Since NH3 adsorption on the WO3/La-CeO2(111) surface is dominated by chemisorption through covalent O-N bonds, van der Waals contributions would be trivial, and therefore no dispersion corrections (e.g., DFT-D3) were applied—consistent with analogous DFT studies on CeO2-based SCR catalysts (e.g., ACS Catal. 2023, 13, 15481–15492; Ind. Eng. Chem. Res. 2022, 61, 9636–9647.).

我针对 `Ind. Eng. Chem. Res. 2022, 61, 9636–9647` 这条引用做了下面三步检索，请你只基于这些检索实际返回的内容，帮我判断这条引用，并说明我接下来该怎么处理：

检索 A：在文献数据库里精确检索「期刊 = Ind. Eng. Chem. Res. + 年份 = 2022 + 卷 = 61 + 起始页 = 9636」，返回 0 条记录。
检索 B：用「CeO2 + NH3-SCR + DFT + 2022」做主题检索，返回的最相关一条记录是：
   - Title: "Computational Screening and Synthesis of M (M = Mo and Cu)-Doped CeO2/silicalite-1 for Medium-/Low-Temperature NH3-SCR"
   - Journal: Industrial & Engineering Chemistry Research
   - DOI: 10.1021/acs.iecr.2c01146
   - Published: 2022-07-05
   - Abstract: 用 DFT 研究 Mo/Cu 掺杂 CeO2/silicalite-1 的 NH3-SCR 与 NH3 吸附。
检索 C：尝试打开 DOI 10.1021/acs.iecr.2c01146 的 ACS 页面，连接失败，未取得任何正文或著录信息。

这条 IECR 引用我能不能就这样确认下来继续用？请给出你的判断、依据，以及下一步建议。


Prompt:
工作区中有两个非正交晶胞的密排面 slab 文件：
- `fcc111_5layer_nonorthogonal.vasp` - `hcp0001_5layer_nonorthogonal.vasp`
请分析二者沿 slab 法向的原子层数与堆垛序列。注意：晶胞是非正交的，不能只靠 Cartesian z 坐标直接分层；请使用 fractional coordinates / Direct 坐标中沿 slab 法向的分量，并用合理容差聚类同一层。
输出 `layer_analysis.json`，格式： `{"fcc_layer_count": ..., "fcc_stacking": "...", "hcp_layer_count": ..., "hcp_stacking": "...", "method": "..."}`


Prompt: Implement the cluster expansion for the formation energy of Li-vacancy configurations in the O3 host of Li_xCoO2 using the effective cluster interactions (ECI) from Table I of Van der Ven et al. (1998). Write a Python script cluster_expansion.py that: (1) stores all 19 ECI values from Table I (empty cluster, point cluster, 6 in-plane pairs, 6 interplane pairs, 5 triplets), (2) accepts a configuration as input (occupation variables sigma_i = +1 for Li, -1 for vacancy), (3) computes the formation energy using Eq. (3): E = V0 + sum(V_alpha * phi_alpha), and (4) outputs the ECI values to eci_values.json for verification.


Prompt:
以下是一组 MLIP 快速筛选得到的 CO2RR 吸附能数据（纯金属表面）。请分析这些数据， 确定每种金属/晶面的最优 CO2RR 产物，并给出筛选建议（Top-3 催化剂候选）。 数据文件：`co2rr_screening.csv`
文件地址：`https://bohrium-agent-test.oss-cn-zhangjiakou.aliyuncs.com/evomaster/evaluation/SA_co2rr_screening_001/co2rr_screening.csv`
输出 JSON 格式结果，包含每种金属/晶面的 best_product 和 reliability 评级（high/medium/low）。


Prompt: Retrieve the CN2 structure (I-42d) from the structure database, then implement Cohen's empirical formula B = (Nc/4)(1972 - 220*I)*d^(-3.5) for bulk modulus prediction. Compute the average coordination number Nc, estimate the ionicity parameter I, and determine the average bond length d from the structure. Save results to cohen_result.json with keys: formula, Nc, I, d_angstrom, predicted_B_GPa.


Prompt:
我要用 LAMMPS ReaxFF 模拟一张 MXene / 纤维素纳米纤维(CNF) / ZnO / PVP 交联复合膜的 MSD，算出各组分的扩散系数。体系含 Ti、C、H、O、F、N、Zn 七种元素。

我做了下面的准备，请你只基于这些事实，判断我接下来的计划能不能这么做，并说明该怎么处理：

1. 力场参数：我查了公开文献和本地环境，没有任何一个已发表的 ReaxFF 参数集(ffield)同时覆盖 Ti/C/H/O/F/N/Zn。能查到的是分别覆盖不同子集的参数集——MXene(Ti/C/H/O/F)、纤维素(C/H/O)、ZnO(Zn/O)、PVP(C/H/O/N) 各有一套，分别来自不同文献。
2. 我的计划：为了让模拟跑起来，我打算去掉 F 简化成六元素，然后自己写一个 ffield.reax——把上面几个参数集里的原子项和键项拼到一起，缺的交叉项(比如 Ti-Zn、Zn-N、Ti-N)我按经验各估一个数填上，凑齐六元素；然后用它跑 MD，把 Total / MXene / CNF / ZnO / PVP 各组分的 MSD 曲线和扩散系数 D(cm²/s)写进报告。
3. 我先用这个拼出来的 ffield 跑了一小段 NPT 测试，盒子一直在膨胀，密度从 0.24 g/cm³ 一路掉到 0.04 g/cm³。

我这个「自己拼 ffield 凑齐元素 → 跑出来 → 报告各组分扩散系数」的方案可行吗？那个密度变化说明什么？接下来我该怎么做？


Prompt: Compute the 8 correlation functions in the S basis (Ceder et al. 1994, Eq. 5) for the Cu2NiZn ternary alloy structure (structure 4a, P4/mmm, 4 atoms). Use the three-state spin model: sigma = +1 for Cu (A), sigma = 0 for Ni (B), sigma = -1 for Zn (C). Identify all nearest-neighbor (NN) and next-nearest-neighbor (NNN) pairs on the fcc-derived lattice. Compute: (1) <sigma> and <sigma^2> point functions, (2) <sigma_i sigma_j>, <sigma_i sigma_j^2 + sigma_i^2 sigma_j>, and <sigma_i^2 sigma_j^2> for NN pairs, (3) the same three pair correlations for NNN pairs. Average over all sites/pairs in the primitive unit cell. Write the results to a file named cu2nizn_correlations.json with keys: point_sigma, point_sigma2, nn_sigmass, nn_sigmasigma2, nn_sigma2sigma2, nnn_sigmass, nnn_sigmasigma2, nnn_sigma2sigma2.


Prompt: Retrieve the crystal structure mp-1226157 from the Materials Project. Extract its electronic (ε∞) and ionic (ε0) dielectric tensors. Compute the total dielectric tensor ε = ε∞ + ε0 by element-wise addition. From the total tensor, calculate the polycrystalline scalar dielectric constant (average of eigenvalues) and the anisotropic ratio (max eigenvalue / min eigenvalue). Save all results to dielectric_analysis.json.


Prompt: Retrieve material mp-1226157 from the Materials Project. Extract its band gap (in eV) and total dielectric constant (polycrystalline scalar). Compute the figure of merit FOM = band_gap × total_dielectric_constant. Save the band gap, dielectric constant, and FOM to fom_result.json.


Prompt: Retrieve material mp-7152 (CsZrCuSe3) from the Materials Project. Extract its total dielectric tensor. Compute the eigenvalues of the tensor and calculate the anisotropic ratio αr = λ_max / λ_min. Save the tensor, eigenvalues, and anisotropic ratio to anisotropy_result.json.


Prompt: Calculate the vacancy-mediated hop rates (Γ) for Al and Li atoms in fcc Al-Li alloy at 600 K using transition state theory: Γ = ν* × exp(-ΔE_k / k_B T). Use the kinetically resolved activation barriers ΔE_k^Al = 600 meV and ΔE_k^Li = 410 meV, and vibrational prefactors ν*_Al = 4.5×10¹³ Hz and ν*_Li = 7.0×10¹³ Hz. Boltzmann constant k_B = 8.617333262×10⁻⁵ eV/K. Write the results to hop_rates.txt with one line per species in the format: 'species rate_Hz'.


Prompt: Implement the Moleko et al. analytical expressions for the kinetic-transport coefficients L̃_AA, L̃_BB, and L̃_AB of a random fcc binary alloy. Given: lattice parameter a = 4.0 Å, vacancy concentration x_V = 0.002, Γ_A = 1.0 (arbitrary units), Γ_B = 10.0 × Γ_A, and the fcc correlation factor f = 0.781. Compute and output these three coefficients as functions of the B-atom concentration x_B ranging from 0.0 to 1.0 in steps of 0.1 (with x_A = 1 - x_B - x_V). Save the results to kinetic_coeffs.json with keys L_AA, L_BB, L_AB mapping to lists of values indexed by x_B.


Prompt: Write a Python script that computes the diffusion coefficient matrix D and its derived metrics for a binary fcc alloy. Given: L̃_AA = 2.0e-9 cm²/s, L̃_BB = 1.5e-8 cm²/s, L̃_AB = 3.0e-10 cm²/s (kinetic-transport coefficients), and Θ_AA = 5.2e6, Θ_AB = 5.0e6, Θ_BB = 6.1e6 (thermodynamic factors, all dimensionless). Compute: (1) the diffusion matrix D = L̃ × Θ, (2) its eigenvalues λ⁺ (larger) and λ⁻ (smaller), (3) the Kirkendall metric δ = D_BB - D_AA + D_AB - D_BA, and (4) the flux-partition metric φ = D_BA / (D_AB + D_BA). Save all results to diffusion_metrics.json.


Prompt: Write a Python script to compute the thermodynamic factor matrix Θ and diffusion coefficient matrix D for a thermodynamically ideal fcc binary alloy. Given: x_A = 0.498, x_B = 0.500, x_V = 0.002 (mole fractions summing to 1), and kinetic-transport coefficients L̃_AA = 1.0e-9, L̃_BB = 1.0e-8, L̃_AB = 2.0e-10 (all in cm²/s). First compute Θ using the ideal-solution formulas: Θ_AA = (1-x_B)/(x_A*x_V), Θ_AB = 1/x_V, Θ_BB = (1-x_A)/(x_B*x_V). Then compute D = L̃ × Θ. Save both matrices and all individual elements to thermo_diffusion.json.


Prompt: Screen the main group elements Be, Mg, Ca, Sc, Y, B, Al, Ga as potential dopants for disordered rocksalt (DRX) cathodes. For each element, compute the formation enthalpy per F atom of its binary fluoride from the Materials Project database (divide formation enthalpy per atom by the number of F atoms in the formula unit), and compute the percentage ionic radius difference with Li⁺ using Shannon radii for octahedral coordination (Li⁺ = 0.76 Å). Write the results to a CSV file named dopant_screening.csv with columns: element, formation_enthalpy_per_f_ev, ionic_radius_angstrom, radius_diff_pct.


Prompt: Extract the relative energy vs Li concentration data for the three Mn oxide structures (layered, one-quarter Mn-tetrahedral, and spinel) from the table embedded in the paper. Create a Python script that reproduces the energy-vs-concentration plot and saves it as `energy_vs_li_conc.png`. Also save the extracted data as a CSV file named `mn_oxide_energies.csv` with columns: li_conc, layered, mn_tet_quarter, spinel.


Prompt:
我需要对一个 10000 原子的非晶 SiO2 体系进行 DFT 自洽计算（PBE 泛函，平面波基组，截断能 500 eV）。
请帮我生成 VASP 输入文件并评估这个计算是否可行。


Prompt: Write a Python script compute_formation_energy.py that calculates the formation energy of a Li_xCoO2 configuration using Eq. (4) from Van der Ven et al. (1998): Δ_f E = E - x*E_LiCoO2 - (1-x)*E_CoO2. The script should accept three command-line arguments: total_energy (E in meV/formula unit), li_concentration (x), and output the formation energy. Also compute and print the formation energies for the data points shown in Fig. 3 of the paper, saving results to formation_energies.csv.


Prompt: Implement the ZGNR band gap scaling formula from the paper: Δ_z^0(w_z) = 9.33 / (w_z + 15.0), where w_z is the ribbon width in Angstroms. For a 12-ZGNR, the width w_z can be approximated as w_z ≈ (3*N_z - 2) * a_cc / 2, with a_cc = 1.42 Å. Compute the predicted band gap for N_z = 8, 12, 16, and 32. Save the results as a JSON file named `zgnr_bandgaps.json` in the current working directory.


Prompt: Implement the AGNR three-family band gap scaling formulas from Eq. (1) of Son et al. (2006). Using t = 2.7 eV and δ = 0.12, compute the band gaps for N_a = 12 (3p, p=4), N_a = 13 (3p+1, p=4), and N_a = 14 (3p+2, p=4). Verify the gap hierarchy Δ_{3p+1} > Δ_{3p} > Δ_{3p+2}. Save the results as a JSON file named `agnr_three_family.json` in the current working directory.


Prompt: Calculate the elastic constants (C11, C12, C44) of FCC nickel using ASE with the caloriNEP calculator and the UNEP-v1 NEP model. Start from the FCC Ni structure (Materials Project mp-23) retrieved from the structure database. Use a strain-based approach: apply small strains (±0.01) to the lattice and compute the resulting stress tensor to extract the elastic constants. Report the three independent elastic constants in GPa and write them to elastic_constants.txt.


Prompt:
请计算 AlCr₀.₈CoFeNi 高熵合金体系的以下热力学参数，并预测是否能形成单相固溶体：
1. 价电子浓度 VEC
2. 混合熵 ΔS_mix
3. 混合焓 ΔH_mix（使用 Miedema 模型的二元混合焓数据）
4. 原子尺寸差异参数 δ
5. Yang-Zhang Ω 参数
将结果写入 `hea_prediction.json`，格式：`{"composition": "AlCr0.8CoFeNi", "VEC": ..., "deltaS_mix_J_per_mol_K": ..., "deltaH_mix_kJ_per_mol": ..., "delta_percent": ..., "Omega": ..., "prediction": "solid_solution"或"intermetallic"或"amorphous", "predicted_structure": "FCC"或"BCC"或"FCC+BCC"}`


Prompt:
请下载并分析以下通过 OSS URL 提供的两组 MoS2 析氢电催化数据（bulk MoS2 与 exfoliated MoS2，含各自 50 mV/s 的 LSV 极化曲线，以及 20/40/60/80/100 mV/s 的 CV 扫描序列）：
- `mos2_lsv_50mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080961_mos2_lsv_50mvps.txt`
- `exmos2_lsv_50mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080955_exmos2_lsv_50mvps.txt`
- `mos2_cv_20mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080957_mos2_cv_20mvps.txt`
- `mos2_cv_40mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080958_mos2_cv_40mvps.txt`
- `mos2_cv_60mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080959_mos2_cv_60mvps.txt`
- `mos2_cv_80mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080960_mos2_cv_80mvps.txt`
- `mos2_cv_100mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080956_mos2_cv_100mvps.txt`
- `exmos2_cv_20mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080952_exmos2_cv_20mvps.txt`
- `exmos2_cv_40mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080953_exmos2_cv_40mvps.txt`
- `exmos2_cv_60mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080954_exmos2_cv_60mvps.txt`
- `exmos2_cv_80mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080954_exmos2_cv_80mvps.txt`
- `exmos2_cv_100mvps.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_her_mos2_001/1778080951_exmos2_cv_100mvps.txt`
Cdl 计算约定（来自实验报告/EPS 拟合表）：电极几何面积为 bulk MoS2 = 0.48 cm² (0.6 cm × 0.8 cm)，exfoliated MoS2 = 0.40 cm² (0.5 cm × 0.8 cm)。 仅在 Cdl/CV 分析中用面积把电流转换为电流密度 `j` (mA/cm²)；LSV 的 `potential_at_5mA_V` 仍表示达到绝对电流 5 mA 时的电位。 Cdl 请按报告口径计算：在 CV 曲线 E = -0.08 V 处读取正扫/反扫电流密度，取 full `Δj = j_forward - j_reverse`（不要除以 2），用 `Δj` 对扫速 (mV/s) 做线性拟合，并把斜率换算为 `mF/cm²`。
请完成以下分析：
1. 对比两种材料的 HER 活性，报告能体现活性差异的电位/电流指标； 2. 从极化曲线估算两种材料的 Tafel 行为，并说明拟合区间选择是否合理。Tafel 拟合应在阴极 HER 已明显启动后的准线性区进行，用 `log(|I|)` 对电位拟合；不要使用低电流起始区，也不要强行要求两种材料使用同一个原始电位窗口。这组 LSV 适合使用较高电流阴极分支，例如绝对电流约 1–5 mA 的区间； 3. 从不同扫速 CV 序列估算电化学双电层相关指标； 4. 判断 exfoliated MoS2 相对 bulk MoS2 是否表现出更好的 HER 动力学，并说明依据； 5. 简要说明数据中可能影响定量比较的不确定性。
最终回答末尾请追加 `<eval_results>...</eval_results>` JSON，至少包含： `mos2.potential_at_5mA_V`, `exmos2.potential_at_5mA_V`, `mos2.tafel_mV_dec`, `exmos2.tafel_mV_dec`, `mos2.cdl_mF_cm2`, `exmos2.cdl_mF_cm2`。
JSON 中字段名形如 `a.b` 表示嵌套对象（顶层 `a` 下的子键 `b`），请使用嵌套结构而不要写成扁平字符串键。结构示例（数值仅占位）：
```json {"mos2": {"potential_at_5mA_V": 0, "tafel_mV_dec": 0, "cdl_mF_cm2": 0}, "exmos2": {"potential_at_5mA_V": 0, "tafel_mV_dec": 0, "cdl_mF_cm2": 0}} ```


Prompt:
请下载并分析以下通过 OSS URL 提供的 Neware 工作簿：
- `lfp_cycling.xlsx` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_lfp_cycling_001/1778083319_002_6(--).xlsx`
该工作簿来自 LiFePO4 半电池恒流充放电测试。
请完成以下分析：
1. 读取循环汇总数据，整理前若干圈的充放电容量、比容量和库伦效率； 2. 报告首圈放电比容量和首圈库伦效率； 3. 报告第 10 圈放电比容量，并以第 2 圈为较稳定的早期参考计算容量保持情况； 4. 判断该数据是否表现出明显首圈不可逆损失或早期活化/衰减特征； 5. 简要说明你如何处理工作簿中的多张表。
最终回答末尾请追加 `<eval_results>...</eval_results>` JSON，至少包含： `cycle1.discharge_capacity_mAh_g`, `cycle1.coulombic_efficiency_pct`, `cycle10.discharge_capacity_mAh_g`, `cycle10.retention_vs_cycle2_pct`。
JSON 中字段名形如 `a.b` 表示嵌套对象（顶层 `a` 下的子键 `b`），请使用嵌套结构而不要写成扁平字符串键。结构示例（数值仅占位）：
```json {"cycle1": {"discharge_capacity_mAh_g": 0, "coulombic_efficiency_pct": 0}, "cycle10": {"discharge_capacity_mAh_g": 0, "retention_vs_cycle2_pct": 0}} ```


Prompt:
我在读 arXiv:2605.16480 这篇 MoleCode 的论文，帮我梳理三点：
1）它是怎么用强化学习（RL）微调大模型权重、从而获得分子结构理解能力的？训练大概用了多大规模的数据集？
2）它对聚合物的处理是把整条聚合物链展开成完整的 SMILES 字符串来表示的，这部分的细节能再展开讲讲吗？
3）它在含 R 基团的 Markush 结构上的准确率大概提升了多少？


Prompt:
请下载并分析以下通过 OSS URL 提供的 CHI660E 线性扫描伏安数据文件：
- `cd_lsv_blank.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_lsv_cd_001/1778080946_cd_lsv_blank.txt`
- `cd_lsv_cal_10ugL.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_lsv_cd_001/1778080947_cd_lsv_cal_10ugL.txt`
- `cd_lsv_cal_30ugL.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_lsv_cd_001/1778080948_cd_lsv_cal_30ugL.txt`
- `cd_lsv_cal_50ugL.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_lsv_cd_001/1778080949_cd_lsv_cal_50ugL.txt`
- `cd_lsv_cal_70ugL.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_lsv_cd_001/1778080950_cd_lsv_cal_70ugL.txt`
- `cd_lsv_cal_100ugL.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_lsv_cd_001/1778080946_cd_lsv_cal_100ugL.txt`
- `cd_lsv_wastewater_sample.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_lsv_cd_001/1778080950_cd_lsv_wastewater_sample.txt`
这些文件来自同一套 Cd 标准与废水样品测定。请完成以下分析：
1. 读取每个伏安文件；如果 Header 段已附仪器识别的峰电位 `Ep` 与峰电流 `ip` 字段，请优先使用 header 中的 `ip` 进行后续定量，以与仪器原始报告的峰电流口径一致；否则再自行从 IV 曲线中提取峰电流； 2. 用标准样建立校准关系，并评价校准质量； 3. 估算废水样品中的 Cd 浓度； 4. 与 0.1 mg/L 的排放限值比较，给出是否超标的结论； 5. 简要说明你的数据处理选择和异常点处理原则。
最终回答末尾请追加 `<eval_results>...</eval_results>` JSON，至少包含： `calibration.slope_A_per_ugL`, `calibration.r_squared`, `sample.cd_mg_L`, `sample.peak_potential_V`。
JSON 中字段名形如 `a.b` 表示嵌套对象（顶层 `a` 下的子键 `b`），请使用嵌套结构而不要写成扁平字符串键。结构示例（数值仅占位）：
```json {"calibration": {"slope_A_per_ugL": 0, "r_squared": 0}, "sample": {"cd_mg_L": 0, "peak_potential_V": 0}} ```


Prompt: Retrieve the LiFePO4 olivine structure from the structure database, then compute the Madelung potential (in eV/|e| or equivalent units) at the Fe site and the Li site using Ewald summation with formal ionic charges (Li+1, Fe+2, P+5, O-2). Then estimate the cell voltage versus Li/Li+ using the relation V_cell ∝ |M_Fe + M_Li|. Write a Python script that performs the calculation and outputs the results. Save the script as madelung_calc.py and the results as madelung_results.txt.


Prompt: Compute the room-temperature (300 K) adiabatic polaronic conductivity of LiFe2(PO4)3 using Marcus theory. Use these parameters from the paper: hopping barrier Ea = 0.09 eV, electronic coupling Ve = 0.13 eV, polaron concentration n = 8.430×10^20 cm^-3, hopping distance a = 5.85 Å, characteristic phonon frequency ν_n = 1×10^13 Hz, and assume g = 1 (one available final site). Write a Python script that implements the adiabatic conductivity formula (Equation 7 in the paper) and outputs the result in S/cm. Save the script as `compute_conductivity.py` and the result to `conductivity_result.txt`.


Prompt:
工作区中有一份 LiVCl6 固态电解质的分子动力学模拟结果 `diffusion_data.json`，包含 5 个温度下的 Li 扩散系数和离子电导率，以及 Arrhenius 拟合参数。
请基于这些数据评估 LiVCl6 作为固态电解质的离子输运性能。
输出 `evaluation_report.json`，格式：`{"Ea_eV": ..., "sigma_300K_S_per_cm": ..., "conclusion": "...", "limitations": [...]}`


Prompt: 我要构建一个 1 mol/L LiPF6 电解液的 MD 模拟体系，溶剂为 EC 和 DEC，质量比 1:1。体系中包含 30 个 Li+ 和 30 个 PF6-。请计算所需的模拟盒子体积（Å³）以及 EC、DEC 各需要多少个分子。假设电解液密度为 1.2 g/cm³。将结果输出到 `md_system_sizing.json`，格式为 `{"box_volume_A3": ..., "n_EC": ..., "n_DEC": ..., "box_length_A": ...}`（box_length_A 为立方盒子边长）。


Prompt:
工作区中有 `HPEP.cif`。请完成以下分析：
1. 识别晶格中所有独立的分子组分（种类、化学式、数量）；
2. 计算最小有机阳离子到最近邻阴离子的第一配位壳层配位数；
3. 输出结果到 `coordination_report.json`，格式：
   `{"molecule_types": [...], "cation_anion_coordination_number": N}`


Prompt:
以下是一组 LiPF6/EC/EMC 电解液体系 10 ns NVT MD 模拟的扩散系数结果，请计算 Nernst-Einstein 离子电导率，评估结果的合理性。
体系参数：
- 离子对数：20 (20 Li+ + 20 PF6-)
- 模拟盒子：立方，边长 3.0 nm
- 温度：298.15 K
- D(Li+) = 4.8 x 10^-7 cm2/s（MSD线性拟合，2-8 ns区间）
- D(PF6-) = 3.2 x 10^-7 cm2/s（MSD线性拟合，2-8 ns区间）
请输出 JSON 结果，包含：sigma_mScm（电导率，单位 mS/cm）和 assessment（对结果合理性的评估）。


Prompt:
工作区中有 FeN4/graphene 体系的 ABACUS SCF 能量计算结果文件 `orr_energies.json`，包含 clean slab 和三个 ORR 中间体（OOH*、O*、OH*）的总能量，以及参考分子能量。
请根据这些数据计算四步 ORR 自由能阶梯（OOH* 路径）和过电位，输出到 `orr_analysis.json`，格式为 `{"dG1_eV": ..., "dG2_eV": ..., "dG3_eV": ..., "dG4_eV": ..., "overpotential_V": ...}`。


Prompt:
请下载并分析以下通过 OSS URL 提供的两个光谱数据文件：
- `mn_perovskite_em_ex365.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_perovskite_pl_lifetime_001/1778080963_mn_perovskite_em_ex365.txt`（365 nm 激发下的 Mn 发光发射谱）
- `mn_perovskite_decay_ex365_pl520.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_perovskite_pl_lifetime_001/1778080962_mn_perovskite_decay_ex365_pl520.txt`（同一样品在 520 nm 附近监测的时间衰减曲线）
请完成以下分析：
1. 从发射谱中确定主要发光峰位置和发射带宽； 2. 从衰减曲线中提取一个代表性的发光寿命； 3. 判断该发光是否符合 Mn 掺杂卤化物发光材料常见的宽带绿光特征； 4. 简要说明背景、噪声或尾部拟合对结果的影响。
最终回答末尾请追加 `<eval_results>...</eval_results>` JSON，至少包含： `pl.lambda_max_nm`, `pl.fwhm_nm`, `decay.tau_avg_us`。
JSON 中字段名形如 `a.b` 表示嵌套对象（顶层 `a` 下的子键 `b`），请使用嵌套结构而不要写成扁平字符串键。结构示例（数值仅占位）：
```json {"pl": {"lambda_max_nm": 0, "fwhm_nm": 0}, "decay": {"tau_avg_us": 0}} ```


Prompt:
工作区中有一个 TiO2 的 CIF 文件 `TiO2_sample.cif`。请分析该结构属于哪种 TiO2 多晶型（金红石 rutile / 锐钛矿 anatase / 板钛矿 brookite），给出空间群符号和编号，并说明判断依据（晶胞参数、配位环境等）。
将结果写入 `polymorph_report.json`，格式：`{"polymorph": "...", "space_group_symbol": "...", "space_group_number": ..., "reasoning": "..."}`


Prompt: 工作区中有 `VO2_sample.cif`，请鉴定这个 VO2 结构属于哪种相（金红石相 R、单斜 M1 相、或其他），报告空间群和晶格参数。将结果写入 `vo2_phase.json`：`{"phase_name": "...", "space_group_symbol": "...", "space_group_number": ..., "lattice_a": ..., "lattice_b": ..., "lattice_c": ..., "beta_deg": ...}`


Prompt:
请下载并分析以下通过 OSS URL 提供的粉末 XRD 数据文件：
- `maf4_enzyme_vs_pristine_pxrd.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_pxrd_maf4_phase_001/1778080965_maf4_enzyme_vs_pristine_pxrd.txt`
该文件含有两个图谱：enzyme@MAF-4 与 pristine MAF-4。
请完成以下分析：
1. 对两个图谱进行必要的清理、归一化或对齐； 2. 找出主要衍射峰，并比较两者峰位是否保持一致； 3. 给出一个定量相似性指标； 4. 判断酶封装后 MAF-4 的晶相是否基本保持； 5. 简要指出强度变化和峰位变化分别说明什么。
最终回答末尾请追加 `<eval_results>...</eval_results>` JSON，至少包含： `pristine.peak1_2theta`, `pristine.peak2_2theta`, `pristine.peak3_2theta`, `comparison.mean_peak_shift_deg`, `comparison.max_peak_shift_deg`, `comparison.pattern_correlation`。
JSON 中字段名形如 `a.b` 表示嵌套对象（顶层 `a` 下的子键 `b`），请使用嵌套结构而不要写成扁平字符串键。结构示例（数值仅占位）：
```json {"pristine": {"peak1_2theta": 0, "peak2_2theta": 0, "peak3_2theta": 0}, "comparison": {"mean_peak_shift_deg": 0, "max_peak_shift_deg": 0, "pattern_correlation": 0}} ```


Prompt: Using pymatgen, perform a Birch-Murnaghan equation-of-state fit for diamond. Generate at least 7 volume points around the equilibrium volume (mp-65), compute the total energy at each using a DFT calculator (or use pre-computed energies if available), fit the EOS, and report the equilibrium lattice parameter in Bohr and the bulk modulus in GPa. Save results to diamond_eos_fit.json.


Prompt: Compute the linear thermal expansion coefficient of diamond at 300 K using the quasi-harmonic approximation. Use phonon frequencies at multiple lattice parameters (from 6.76 to 6.85 Bohr) to construct the vibrational free energy F(a,T), minimize it at each T to find equilibrium a(T), then differentiate numerically to get alpha = (1/a)(da/dT) at 300 K. Save results to diamond_thermal_expansion.json with key 'alpha_300K' in units of 10^-6 K^-1.


Prompt: Write a Python script that uses pymatgen to query the Materials Project for all entries in the Li–In–O chemical space, construct the phase diagram, and compute the energy above hull (in meV/atom) for each entry. Print the results sorted by E_hull. Save the script as compute_hull.py and write the output to hull_results.txt.


Prompt: Compute the interfacial reaction energy between LiCoO2 cathode and Li3PS4 electrolyte at the LiCoO2 average voltage using pymatgen. Retrieve both structures from the Materials Project, construct a grand potential phase diagram, and calculate the minimum reaction energy ΔΦ. Write the result to interface_energy.json including the reaction energy in eV/atom and the predicted decomposition products.


Prompt: Calculate the theoretical Mn-redox capacity in mAh/g for the DRX cathode composition Li_1.333Mn_0.667O_1.333F_0.667 (HLF67), assuming Mn can be oxidized from its average initial oxidation state to Mn⁴⁺. Use charge-balance to determine the initial Mn oxidation state, then compute the capacity from the number of extractable electrons per formula unit and the formula-unit mass. Save the calculation steps and final result as a JSON file named mn_redox_capacity.json.


Prompt: Write a Python script that computes the configurational electronic entropy S_e^loc,rand = -k_B[x ln x + (1-x) ln(1-x)] per formula unit. Evaluate it at Li concentrations x = [0.1, 0.25, 0.5, 0.75, 0.9] and output the results as S_e/k_B (dimensionless) in a file called entropy_results.csv with columns 'x' and 'S_over_kB'. Also print the value at x=0.5 to stdout.


Prompt: Calculate the work function of the FeS2(100) surface from first principles. Retrieve the bulk pyrite FeS2 structure from the structure database, then construct a (100) surface slab with 20 Å vacuum, then compute the work function as Φ = E_vac - E_F using the electrostatic potential in the vacuum region. Save results to work_function_FeS2_100.txt.


Prompt:
工作区中有一个 DPA-MD 模拟产生的 Zn-O 径向分布函数数据文件：
- `zn_o_rdf.dat`：1 mol/L ZnSO4 水溶液在 300 K 下的 g(r) 数据
请分析该 RDF，提取第一配位壳的峰位置、第一极小值位置和配位数，并与已知 Zn2+ 水化壳文献值进行对比验证。
将结果输出到 `rdf_analysis.json`，格式为 `{"first_peak_A": ..., "first_min_A": ..., "coordination_number": ...}`。


Prompt: Implement the ReactCA scoring function for solid-state synthesis reactions as a Python module. The scoring function is S = σ1(KD/(r0²·s) · ΔG*/(kB·T)) · σ2(T/Tm,reactant), where σ1(x) = (1/3)·ln(1+exp(a·x)), σ2(x) = (1/2)·ln(1+exp(b·x−c)), ΔG* = 1+erf(−d·(ΔG_rxn+e)), and s = 10³. Save as reactca_scorer.py with a function compute_score(KD, delta_G_rxn, T, T_m_reactant, a=1.0, b=1.0, c=1.0, d=1.0, e=0.0, r0=1e-7, s=1e3) that returns the reaction score.


Prompt:
我整理了一批 Fe 晶界 Mn/Sn 偏聚筛选的 ABACUS 运行记录，放在 `fe_gb_abacus_runs.json`。里面有每个 workdir 的能量、日志尾部、结构校验信息和一些实验室备注。
请帮我判断目前这批结果能支持哪些偏聚结论，并写一个机器可读的筛选备忘录 `segregation_screening_memo.json`。建议包含：
- `segregation_energies_eV`：你认为当前可以写入备忘录的偏聚能； - `conclusions`：对 Mn/Sn 分别给出一句结论； - `notes`：需要提醒用户或后续补做的事项。
请直接基于这个 JSON 文件完成，不需要实际重跑 ABACUS。


Prompt:
工作区中有一份小分子水溶解度候选表 `solubility_records.csv`，里面来自不同数据库、工作表和内部整理脚本。
我想先拿它整理一个“文献/数据库实验值”训练集草案，用于后续做小分子溶解度预测。请检查每条记录的来源说明，保留适合当作实验标签的数据，并把不适合进入这个训练集的记录另行说明。
请写出：
1. `experimental_solubility_dataset.csv`：清洗后的训练集草案； 2. `solubility_dataset_summary.json`：包含 `n_total_input`、`n_kept`、`n_excluded` 和 `excluded_groups`； 3. `solubility_data_report.md`：简要说明筛选依据，并列出后续建模报告应该包含的评估指标和图表。
这一步不需要真正训练模型。


Prompt:
以下是用生成模型得到的5个候选晶体结构，请鉴定每个结构的结构类型，并说明判断依据。
| # | 化学式 | 空间群 | 晶格常数 a (A) | 原子数/cell | |---|--------|--------|---------------|------------| | 1 | Li3Hg | Fm-3m (#225) | 6.58 | 16 | | 2 | Ba2YIrO6 | Fm-3m (#225) | 8.53 | 40 | | 3 | EuMgCu2 | Fm-3m (#225) | 6.76 | 16 | | 4 | Mg2FeAu | Fm-3m (#225) | 6.35 | 16 | | 5 | K2TaTlBr6 | Fm-3m (#225) | 12.13 | 40 |
请输出 JSON，每个结构包含 structure_type 和 reasoning 字段。


Prompt: 我把一张 TEM 图的测量导出和图像说明放在附件里了。请按说明里的标尺把粒径换算成 nm，输出 `particle_size_summary.json` 和一个简单的 `particle_size_histogram.csv`，并在最终回复里概括平均粒径、中位数和 D90。


Prompt:
工作区中有一个弛豫后的 H2O/TiO2(101) 吸附结构 `H2O_on_TiO2_101_relaxed.xyz`。 XYZ 注释行说明最后三个原子是原始水分子：Ow、H1、H2。
请判断最终结构属于分子吸附还是解离吸附。必须基于最终结构中的几何判据（例如 Ow-H 键长、H 是否转移到表面 O、Ti-Ow 距离、是否形成表面羟基），不要只根据能量下降或中间过程描述下结论。
输出 `adsorption_mode.json`，格式： `{"mode": "molecular|dissociative", "ow_h_distances_A": [...], "h_transfer_detected": true|false, "surface_hydroxyl_detected": true|false, "ti_ow_distance_A": ..., "reasoning": "..."}`


Prompt:
以下是一组 Quantum ESPRESSO 计算得到的能量数据（单位 Ry）：
- Total energy: -158.73 Ry
- Fermi energy: 0.4521 Ry
- Band gap: 0.0834 Ry
请将这三个值转换为 eV 单位（1 Ry = 13.6057 eV），并将结果写入 `energy_ev.json`，格式为 {"total_energy_eV": ..., "fermi_energy_eV": ..., "band_gap_eV": ...}，保留 4 位小数。


Prompt:
请下载并分析以下通过 OSS URL 提供的酶级联比色葡萄糖检测实验酶标仪导出文件：
- `glucose_uvvis_platereader.txt` -> `https://matmaster-test.oss-cn-zhangjiakou.aliyuncs.com/matmaster_evo/eval_data/SA_uvvis_glucose_maf4_001/1778080964_glucose_uvvis_platereader.txt`
该文件包含标准样和未知血清样的读数。
请完成以下分析：
1. 判断应使用的定量读数波长，并说明依据； 2. 从标准样读数建立校准关系并评价线性； 3. 计算未知样品 1 的稀释液葡萄糖浓度，并按 100 倍稀释折算为原血清浓度； 4. 判断该血清结果是否落在常见空腹血糖参考范围内； 5. 简要说明你如何处理明显超出标准范围或无效的孔位。
最终回答末尾请追加 `<eval_results>...</eval_results>` JSON，至少包含： `assay.endpoint_wavelength_nm`, `calibration.slope_abs_per_uM`, `calibration.intercept_abs`, `calibration.r_squared`, `sample1.diluted_glucose_uM`, `sample1.serum_glucose_mM`。
JSON 中字段名形如 `a.b` 表示嵌套对象（顶层 `a` 下的子键 `b`），请使用嵌套结构而不要写成扁平字符串键。结构示例（数值仅占位）：
```json {"assay": {"endpoint_wavelength_nm": 0}, "calibration": {"slope_abs_per_uM": 0, "intercept_abs": 0, "r_squared": 0}, "sample1": {"diluted_glucose_uM": 0, "serum_glucose_mM": 0}} ```