adf:voltageprofileofdischarging2020
使用ReaxFF中巨正则蒙特卡洛(GCMC)方法模拟锂硫电池锂化放电过程的电压特性
英文教程
英文教程导读
准备知识
该功能适用于版本号大于AMS2020
熟悉ReaxFF建模、保存任务、提交任务等基本操作
熟悉巨正则系综的含义更好
知道如何在命令行里面运行命令(本教程需要在命令行运行一次分析命令,得到电压特性数据)
-
阅读说明
该教程非常详细,实际上初学者也可以照着说明一步一步做出来,当然结果会有差异,这是正常的。本教程主要包含如下几个小节:
从外部导入Alpha硫晶体的*.cif文件,分别采用NPT系综低温弛豫(AMSinput > ReaxFF > Main > Task: Molecular Dynamics)、能量最小化(AMSinput > ReaxFF > Main > Task: Geometry Optimization)来得到其平衡结构。需要说明的是,教程中能量最小化得到的能量是-8540 kcal/mol,实际我们得到的能量可能更低。
计算Li原子的化学势:计算体心立方锂晶体的8*8*8超胞(共512个Li原子)的能量,除以512,即得到Li原子的化学势-37.70kcal/mol
对Alpha硫晶体进行巨正则系综蒙特卡洛模拟
从上一步GCMC模拟生成的*.results/ams.rkf文件中,抽取电压特性等数据:将教程中可以下载的LiVoltageProfile.py文件放到该任务所在文件夹内,在AMS环境变量对当前命令行生效的前提下,在命令行执行:
$AMSBIN/plams LiVoltageProfile.py -v resultsdir=*.results
注意:
对Windows来说,在命令行执行的具体方式为:执行AMS202*.*\ams_command_line.bat,将自动读取环境变量。输入sh回车,之后即可像Linux一样,进行命令行工作
将上面命令行中的*替换为具体作业的名字,例如$AMSBIN/plams LiVoltageProfile.py -v resultsdir=GCMC.results
输入该命令后回车,过几分钟后,会输出4列数据。其中第一列是Li原子的个数,第二列是LixS的x值,第三列是电压,第四列是体积(Å3),如下图所示:
用户可以自行用这中间两列数据生成教程中的曲线。
adf/voltageprofileofdischarging2020.txt · 最后更改: 2024/06/14 07:46 由 liu.jun