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

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

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

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

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

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)

“原子丢失” 指代不明

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