ReaxFF教程:使用ReaxFF中巨正则蒙特卡洛(GCMC)方法模拟锂化放电过程的电压特性

Posted · Add Comment

英文教程

英文教程阅读说明

准备知识
  1. 该功能版本号≥2017.211,下载地址:https://www.scm.com/support/downloads/development-snapshots/,svnRevision大于r62841的都可以
  2. 熟悉ReaxFF建模、保存任务、提交任务等基本操作
  3. 熟悉巨正则系综的含义更好
  4. 知道如何在命令行里面运行命令(本教程需要在命令行运行一次分析命令,得到电压特性数据)
  5. 知道如何将数据从文本格式转移到Excel中,并生成曲线,参考:如何从图谱中导出数据到Excel表格
阅读说明

该教程非常详细,实际上初学者也可以照着说明一步一步做出来,当然结果会有差异,这是正常的。本教程主要包含如下几个小节:

  1. 从外部导入Alpha硫晶体的*.cif文件,分别采用NPT系综低温弛豫(ADFinput > ReaxFF > Main > Task: Molecular Dynamics)、能量最小化(ADFinput > ReaxFF > Main > Task: Geometry Optimization)来得到其平衡结构。需要说明的是,教程中能量最小化得到的能量是-8535.99 kcal/mol,实际我们得到的能量可能更低。
  2. 计算Li原子的化学势:计算体心立方锂晶体的8*8*8超胞(共512个Li原子)的能量,除以512,即得到Li原子的化学势-37.70kcal/mol
  3. 对Alpha硫晶体进行巨正则系综蒙特卡洛模拟
  4. 从上一步GCMC模拟生成的*.rxkf文件中,抽取电压特性等数据:将LiVoltageProfile.py文件和*.rxkf放到同一个文件夹内,在ADF环境变量对当前命令行生效的前提下,在命令行执行:startpython LiVoltageProfile.py gcmc_test.rxkf。将生成一个文件,名为voltage_profile.out。该文件实际上是一个文本文件,给出了4列数据,其中第二列是LixS中的x值,第三列是电压,第四列是体积变化。用户可以自行用这两列数据生成教程中的曲线

上面这几步的计算模拟文件(请点击此处下载),在阅读困难的情况下,可以辅助您阅读教程,但由于GCMC这一步很耗时间(对应文件gcmc_test.adf),可能需要一天左右才能完成,因此设置的参数数值放宽了条件,没有教程里面那么严格,因此这一步仅供参考。

 
  • 标签

  • 关于费米科技

    费米科技以促进工业级模拟与仿真的应用为宗旨,致力于推广基于原子级别模拟技术和基于图像模型的仿真技术,为学术和工业研究机构提供研发咨询、软件部署、技术攻关等全方位的服务。费米科技提供的模拟方案具有面向应用、模型新颖、功能丰富、计算高效、简单易用的特点,已经服务于众多的学术和工业用户。

    欢迎加入我们!(点击链接)

  • 最近更新

  • 联系方式

    • 留言板点击留言
    • 邮箱:sales_at_fermitech.com.cn
    • 电话:010-80393990
    • QQ: 1732167264
  • 订阅费米科技新闻

    • 邮件订阅:
      您可以使用常用的邮件地址接收费米科技定期发送的产品更新和新闻。
      点击这里马上订阅
    • 微信订阅:
      微信扫描右侧二维码。
  •