计算化学公社

标题: 求助：水系锌离子电解液模拟的退火问题以及平衡相到产生相的延续问题 [打印本页]

作者
Author: AnnieXu 时间: yesterday 12:11
标题: 求助：水系锌离子电解液模拟的退火问题以及平衡相到产生相的延续问题
老师们好，目前研究各种添加剂对水系锌离子电解液的改善情况，模拟时遇到以下疑惑：
1. 到底需要进行退火操作吗，我目前需要研究不同温度的体系，如果需要进行退火，这样设置可行吗
Tcoupl  = V-rescale
tau_t = 0.5
tc_grps = system
ref_t = 298.15
annealing = single
annealing-npoints = 2
annealing-time = 0 1000000
annealing-temp = 0 298.15
;
Pcoupl    = Berendsen
pcoupltype = isotropic
tau_p = 5.0
ref_p = 1.0
compressibility = 4.5e-5
;
gen_vel  = no
gen_temp = 298.15
gen_seed = -1

2. 预平衡用V-rescale以及Berendsen，产生相用V-rescale以及PR，相应的 tau_t 和 tau_p 改为0.2和2，预平衡与产生相均设置 gen_vel  = no 。
我查看了论坛上的帖子http://bbs.keinsci.com/forum.php ... 5&fromuid=83569，此时平衡相切换到产生相是加上 -t npt.cpt 更好吗？
若进行退火操作，则直接将其与预平衡合成一步？

作者
Author: yuzc 时间: yesterday 14:47
1. 这种退火是不可取的，你的目的应该是充分松弛盒子，并且让SSIP, CIP, AGG的比例达到一个合理的水平。对于松弛盒子的目的，你应该做高温到低温的退火。同时应该关注盒子体积变化率，避免出现真空区。
2. 原则上退火最后个温度平台区间可以用来做产生相。现在有C-rescale，可以不用PR了，除非你要各向异性的盒子变化。

作者
Author: AnnieXu 时间: yesterday 15:45

yuzc 发表于 2026-1-22 14:47
1. 这种退火是不可取的，你的目的应该是充分松弛盒子，并且让SSIP, CIP, AGG的比例达到一个合理的水平。对 ...

感谢您的回复，请问高温到低温的退火是从0K-高温-模拟实际温度吗，这个高温的选择有依据吗？

作者
Author: sobereva 时间: 13 hour ago
不是必然非得先拉到高温再降温，只不过这样可能能比直接从低温升到目标温度用更短的模拟时间达到平衡，具体视体系而定。
对于你的情况，应当在较短时间内把参考温度拉到目标温度（中途不是必须经过高温阶段），然后维持足够长的时间使得体系达到平衡

续跑时一律接着之前的cpt文件

作者
Author: AnnieXu 时间: 5 hour ago

sobereva 发表于 2026-1-23 01:57
不是必然非得先拉到高温再降温，只不过这样可能能比直接从低温升到目标温度用更短的模拟时间达到平衡，具体 ...

谢谢老师

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