# 使用 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晶胞优化的结构

可以在上述代码的最后加上以下代码来输出优化后的结构：

# 输出优化后的结构
print("优化后的晶格参数：", box.get_cell())
print("优化后的原子坐标：")
for atom in box:
    print(atom.symbol, atom.position)

其中，`box.get_cell()`可以获取优化后的晶格参数，`for atom in box:`可以遍历优化后的原子坐标并输出。