这里会显示出您选择的修订版和当前版本之间的差别。
后一修订版 | 前一修订版 | ||
adf:ereaxff [2021/08/09 22:38] – 创建 liu.jun | adf:ereaxff [2024/02/05 13:06] (当前版本) – liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ======eReaxFF:锂离子电池中碳酸乙烯酯(EC)分解的模拟====== | + | ======eReaxFF-锂离子电池中碳酸乙烯酯(EC)分解的模拟====== |
- | 版本要求2021。 | + | 版本要求2021。eReaxFF是将电子与原子一样,视为一个普通粒子(当然质量、电荷与原子有所不同),进行分子动力学模拟。 |
=====一、模型===== | =====一、模型===== | ||
首先创建电子、EC、Li原子(实际上是阳离子,但是在模型参数设置的时候,才进行电荷设置),然后将三者生成均匀混合物。 | 首先创建电子、EC、Li原子(实际上是阳离子,但是在模型参数设置的时候,才进行电荷设置),然后将三者生成均匀混合物。 | ||
====1,创建电子==== | ====1,创建电子==== | ||
- | 在AMSinput窗口任意添加一个原子,然后File - Export Coordinates - .xyz,保持文件名例如为El。打开生成的El.xyz文件,将内容中的元素符号名字改为El,并保存文件。ReaxFF中将El识别为电子,但在计算中仍然需要用户指定El的电荷为-1(详见后文)。 | + | 在AMSinput窗口任意添加一个原子,然后File - Export Coordinates - .xyz,保持文件名例如为El。打开生成的El.xyz文件,将内容中的元素符号改为El,并保存文件。ReaxFF中将El识别为电子,但在计算中仍然需要用户指定El的电荷为-1(详见后文)。 |
====2,创建Li、EC==== | ====2,创建Li、EC==== | ||
在AMSinput窗口创建一个Li原子,然后类似如上导出Li.xyz文件;创建EC分子,并使用窗口右下角齿轮优化分子结构,然后导出EC.xyz。 | 在AMSinput窗口创建一个Li原子,然后类似如上导出Li.xyz文件;创建EC分子,并使用窗口右下角齿轮优化分子结构,然后导出EC.xyz。 | ||
行 11: | 行 10: | ||
{{ : | {{ : | ||
- | 现在有三个xyz文件备用。 | + | 现在有了三个xyz文件备用。 |
====3,创建EC、Li、电子混合物==== | ====3,创建EC、Li、电子混合物==== | ||
行 42: | 行 41: | ||
{{ : | {{ : | ||
- | 保存并提交作业。初次提交,有可能出现“Unreasonably large charges produced by ACKS2”的报错,这时候可以考虑把El和Li放到一起,甚至在空间上略微重叠,这样就可以避免这个问题了。 | + | 保存并提交作业。 |
+ | <color blue> | ||
=====结果查看===== | =====结果查看===== | ||
SCM - Movie可以看到电子与分子中的原子的运动轨迹: | SCM - Movie可以看到电子与分子中的原子的运动轨迹: | ||
{{ : | {{ : | ||
+ | |||
+ | 经过几百ps后,可以观察到乙烯、二氧化碳等产物。 | ||
+ | ====键级==== | ||
+ | 在AMSinput中的键级,未经计算,因此其键级是图形窗口根据简单的价电子规则而猜测出来的。在Movie中是经过计算的,因此鼠标选中两个原子,窗口左下角则会显示两个原子之间的键级。 | ||
+ | =====参考文献===== | ||
+ | Md Mahbubul Islam and Adri C. T. van Duin, //Reductive Decomposition Reactions of Ethylene Carbonate by Explicit Electron Transfer from Lithium: An eReaxFF Molecular Dynamics Study//, J. Phys. Chem. C 2016, 120, 48, 27128–27134 | ||
+ |