Q1：朱老师，请问我做AIMD的体系跑MOF材料，在NVT下设置的初始和末态温度为100K，0.5fs的步长，跑200步，OSZICAR中的T就从100开始掉到50左右，请问我应该怎么改进达到设置的温度啊？
A：一开始会下降，过0.5ps会平衡；
温控没真正起作用，耦合太弱 ， 体系在前几十步把动能转成了势能；
步数太少导致你看到的几乎全是瞬时波动。
Q2：朱老师，电荷密度差图的等值面数值是不是（0.0004 e/ų）？
A：对，截图里 Isosurface level = 0.0004 就是画等值面的“阈值”；
单位取决于加载的“标量场”是什么物理量：
CHGCAR / 差分 CHGCAR：一般就是 e/ų
LOCPOT：电势（eV）
ELF：无量纲（0–1）
Q3：朱老师，这个一维线性平均和一维宏观平均分别什么含义，有什么区别？
A：这个是平面平均值：
一维线性平均：把沿 z 方向的分布趋势提取出来；
一维宏观平均：在线性平均的基础上，把原子尺度的高频振荡抹掉；
Q4：朱老师，柔性具有呼吸效应的MOF材料结构优化采用PE方法，晶格畸变超过6%是不是优化失败了？
A：不一定，可能是真的发生了呼吸相变;
也可能是 方法或者设置把结构拉偏了，导致结构发生假坍塌，假膨胀；
结构是否化学合理：
金属节点配位有没有断键；
有没有出现明显不合理的短接触；
有机配体有没有扭曲变形。
Q5：朱老师，我在做分子动力学时，前面都是正常运行，但是在128步开始结构就突然崩溃了（前一个127步还很正常），温度也只能显示为*，这一般是什么原因啊？
A：某一步开始原子距离变得极端不合理 → 力或速度瞬间爆炸 → 温度大到超输出格式上限
AIMD 每一步都需要自洽电子结构来给出靠谱的力；
对含 H、柔性键角/扭转很多的 MOF 来说，2 fs 的 AIMD 很容易数值不稳定。
