使用 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)。修改上述代码,输出ASE晶胞优化的结构为cif
时间: 2024-03-07 15:46:44 浏览: 185
ASE12.5.1_资料_性能和调优
要将优化后的晶体结构输出为cif文件,可以使用ASE中的`write()`函数。具体来说,可以在代码最后添加以下代码:
```
from ase.io import write
# 输出优化后的结构为cif文件
write('optimized_structure.cif', box, format='cif')
```
其中,`write()`函数的第一个参数是要输出的文件名,第二个参数是要输出的结构,第三个参数是输出格式,这里选择了cif格式。这样,优化后的晶体结构就会被保存到名为`optimized_structure.cif`的文件中。
阅读全文