计算化学公社

标题: AIMD预测材料不稳定该怎么给审稿人回复稳定性对其应用的影响 [打印本页]

作者
Author: 2422759918 时间: 2024-4-14 18:21
标题: AIMD预测材料不稳定该怎么给审稿人回复稳定性对其应用的影响
想问下大家，返修一篇论文，预测四层的二维材料用于太阳能电池，但是审稿人问到预测的材料热力学稳定性怎么样，稳定性如何影响器件的性能？我尝试跑了一下分子动力学，但是不太稳定，比如有的键变长，个别原子丢失，这是不是说明材料不稳定，那我应该怎么和审稿人解释啊，是不是不稳定就不能发了，有没有大佬支个招啊，真的想把它发出去。

作者
Author: sobereva 时间: 2024-4-15 04:40
“原子丢失” 指代不明

如果常温下几十ps内结构特征都维持不住，而且确信动力学模拟设置无误、计算过程无异常（如没有出现SCF基本没收敛），说明当前结构就根本没有实际意义。

作者
Author: 2422759918 时间: 2024-4-15 10:21

sobereva 发表于 2024-4-15 04:40
“原子丢失” 指代不明

如果常温下几十ps内结构特征都维持不住，而且确信动力学模拟设置无误、计算过程 ...

IMCAR参数如下,麻烦老师看下有没有什么问题

ISTART =  1
ISPIN  =  1
LREAL  = Auto
PREC =  A
LWAVE  = .F.
LCHARG = .F.
ADDGRID= .TRUE.
ISMEAR =  0
SIGMA  =  0.05
EDIFF  =  1E-05
ISIF = 2
ALGO = N
IBRION =  0          (Activate MD)
NSW =  5000       (Max ionic steps)
EDIFFG = -1E-02       (Ionic convergence; eV/A)
POTIM  =  1          (Timestep in fs)
TEBEG  =  300    (Start temperature K)
TEEND  =  300    (Final temperature K)
MDALGO =  1
ANDERSEN_PROB = 0.5       (Andersen Thermostat)
ISYM =  0       (Switch symmetry off)
Nwrite =  1       (For long MD-runs use Nwrite=0 or Nwrite=1)

作者
Author: logzzz 时间: 2024-4-15 12:46
只有INCAR能说明什么呢。你有用合理的超胞去跑么？如何设置的K点。你分子动力学最后得到的能量与时间的变化图谱是什么。不行，你就找找网上的教程，有很多。

作者
Author: 2422759918 时间: 2024-4-15 14:08
本帖最后由 2422759918 于 2024-4-15 14:09 编辑

logzzz 发表于 2024-4-15 12:46
只有INCAR能说明什么呢。你有用合理的超胞去跑么？如何设置的K点。你分子动力学最后得到的能量与时间的变化 ...

好的，k点是Gmma点111，超胞大概140多个原子，能量随时间的变化基本稳定在300k左右，上下跳动小于4k，感觉是达到平衡了，就是这个结构，有部分原子之间的键变长了，用vesta打开都不成键了

作者
Author: logzzz 时间: 2024-4-15 16:46

2422759918 发表于 2024-4-15 14:08
好的，k点是Gmma点111，超胞大概140多个原子，能量随时间的变化基本稳定在300k左右，上下跳动小于4k，感 ...

如果不稳定的话，会有比较大的跳点。你需要看到是体系的能量，不是温度。你设置了300到300的温度，它肯定在300附近波动啊。

作者
Author: 2422759918 时间: 2024-4-15 19:57
本帖最后由 2422759918 于 2024-4-15 20:00 编辑

你好，跑了3000步，能量变化如下图，麻烦您帮忙看看，是不是跑的步数太小了，只跑了3000步（3000fs）。
(, 下载次数 Times of downloads: 18)

欢迎光临计算化学公社 (http://bbs.keinsci.com/)