修改代码:要求输出ASE晶胞优化后的结构。mol_ini = Molecule.from_file("Fe32C9/60_1/POSCAR") mol_opt = Molecule.from_file("Fe32C9/60_1/CONTCAR") mol_refined, chemenv = refine(mol_ini, submols) # 使用 ASE 来完成晶胞及原子结构的进一步优化 mol = mol_ini.clone() p1 = spdkit.to_ase_atoms(mol) calc_elastic = ElasticModel(mol, chemenv) p1.set_calculator(calc_elastic) box = Mushybox(p1) # dyn = FIRE(p1, maxstep=0.1) dyn = FIRE(box, maxstep=0.1) dyn.run(fmax=0.001, steps=100)
时间: 2024-03-06 22:52:15 浏览: 155
您需要将代码中的最后一行修改为以下内容:
```
ase.io.write("Fe32C9/60_1/optimized_POSCAR", box.to_atoms())
```
这将输出ASE优化后的结构到名为“optimized_POSCAR”的文件中。如果需要输出其他格式的文件,可以将`optimized_POSCAR`替换为所需格式的文件名,如`optimized.xyz`或`optimized.cif`。
阅读全文
相关推荐














