AIMD预测材料不稳定该怎么给审稿人回复稳定性对其应用的影响

2422759918

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

2422759918

你好，跑了3000步，能量变化如下图，麻烦您帮忙看看，是不是跑的步数太小了，只跑了3000步（3000fs）。

logzzz

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

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

2422759918

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

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

logzzz

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

2422759918

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)

sobereva

“原子丢失” 指代不明

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