这是本文档旧的修订版!
S 晶体中,随机插入一定数量的 Li 原子,形成 Li$_{0.4}S,然后通过退火操作形成非晶结构,并基于该结构进行 1600K 恒温分子动力学模拟,从该过程中计算 Li 原子扩散系数
S 晶体 *.cif 文件下载(后解压得到*.cif)。AMSinput → File → Import Coordinates读取 *.cif 文件,总计128 个 S 原子:
另开一个AMSinput并创建一个Li原子,File → Export Coordinates保存为 Li.xyz。Edit → Builder,添加 51 个 Li原子,注意最小原子间距从默认 2.5Å 改为 1Å,否则Li原子插不进去:
优化过程这里不再详细展开,用户自行进行即可,Task选择 Geometry Optimization,力场选择LiS.ff,并在Details → Geometry Optimization窗口中,勾选Optimize lattice,保存作业并运行。收敛后,Movie的最后一帧用于后续的退火操作。
基于AMSinput中的新结构,(为了看起来更舒服,可以Bonds→Guess bonds重新猜测键级),设置退火参数:
退火结果:SCM → Movie → MD Properties → Temperature:
基于最后一帧结构,再次如前优化结构与晶格常数。