计算化学公社

标题: 关于NVT平衡后电场模拟失败的主题 [打印本页]

作者
Author: 19821838673 时间: 2025-11-24 09:57
标题: 关于NVT平衡后电场模拟失败的主题
各位老师好~我做的是一个电池体系，正极材料是碘单质和双层石墨烯，电池充放电过程中碘单质会和碘离子反应生成碘三离子，对应阳抗衡离子是锌离子，溶液环境是2M硫酸锌溶液。我的体系中还存在一种两性离子（即阴阳离子同时存在一个结构中），如图1所示。我将两性离子吸附在了电池正极石墨烯表面，并对石墨烯和两性离子层都做了固定。目的是想看碘三离子在电池环境的电场作用下被吸附到两性离子层的N+离子上的过程，如图2。这里想请问一下大家，我做了NVT平衡之后，我准备再在Z负方向上施加一个0.5V/nm的电场，但是却出现如图3中的报错，请问会是什么原因呢？我检查了轨迹没有问题，能量最小化也收敛得很好（emtol=100)，并且NVT平衡过程我用的退火升温到常温，也排除原子速度过大导致系统崩溃的原因。此前尝试过多种调整方式，有试过把电场和冻结去掉，再次进行能量最小化，但是结果没有任何变化，显示体系的势能依旧是2.99407e+10。实在找不到原因，特来向各位前辈请教！
最后在附件里提供源文件，老师们有空的话可以顺便检查下！

作者
Author: sobereva 时间: 2025-11-25 01:55
我没时间具体看文件，如果只是加了电场后模拟才崩溃，先尝试带着电场做em之后再MD。如果还不行，把电场从很小一点点往大加，反复对比、测试找原因

作者
Author: 19821838673 时间: 2025-11-25 09:09

sobereva 发表于 2025-11-25 01:55
我没时间具体看文件，如果只是加了电场后模拟才崩溃，先尝试带着电场做em之后再MD。如果还不行，把电场从很 ...

好的卢老师，我尝试下把电场改小，再来跑模拟，感谢您的答复

作者
Author: 19821838673 时间: 2025-11-25 11:10

sobereva 发表于 2025-11-25 01:55
我没时间具体看文件，如果只是加了电场后模拟才崩溃，先尝试带着电场做em之后再MD。如果还不行，把电场从很 ...

卢老师，我改小了电场，而且像您说的，不去电场的条件下跑了能量最小化，总能量收敛至<100后跑的field，但还是出现一模一样的结果。我怀疑是我的field.mdp参数设置的不合理，您能抽空帮我看一下这个文件吗？

作者
Author: sobereva 时间: 2025-11-25 11:27

19821838673 发表于 2025-11-25 11:10
卢老师，我改小了电场，而且像您说的，不去电场的条件下跑了能量最小化，总能量收敛至

乱七八糟的设置先都去掉，诸如wall、冻结
二维周期性体系不要用DispCorr = EnerPres
甭写完全多余的fourierspacing = 0.16
rvdw设得明显偏小
尝试用三维周期性

作者
Author: 19821838673 时间: 2025-11-25 13:48

sobereva 发表于 2025-11-25 11:27
乱七八糟的设置先都去掉，诸如wall、冻结
二维周期性体系不要用DispCorr = EnerPres
甭 ...

好的，感谢卢老师，我按照您说的试一下

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