计算化学公社

标题: CP2K跑AIMD温度为350K和500K，体系300-400个原子，怎么设time constant? [打印本页]

作者
Author: Oliviaw 时间: 2024-6-30 21:06
标题: CP2K跑AIMD温度为350K和500K，体系300-400个原子，怎么设time constant?
一个300-400个的host-guest体系，想讨论热稳定性，设置温度分别为350K和500K，跑10ps，步长1fs，PBE-D3BJ/DZVP-MOLOPT-SR-GTH，CSVR，WAVELET，跑AIMD。这种情况需要设置time constant多少合适呢？

作者
Author: sobereva 时间: 2024-7-1 18:14
一般可以先尝试200 fs，看控温效果，如果温度始终控得不那么接近期望温度，可以再用更小的使热浴效果更显著

作者
Author: Oliviaw 时间: 2024-7-3 22:35

sobereva 发表于 2024-7-1 18:14
一般可以先尝试200 fs，看控温效果，如果温度始终控得不那么接近期望温度，可以再用更小的使热浴效果更显著

我试的时候是用的time constrant 100，设定温度为350K，跑了快4 ps，温度差不多稳定在320 K。请问这个控温是不是效果可以？

作者
Author: sobereva 时间: 2024-7-5 08:38

Oliviaw 发表于 2024-7-3 22:35
我试的时候是用的time constrant 100，设定温度为350K，跑了快4 ps，温度差不多稳定在320 K。请问这个控 ...

还行。想更快平衡在350附近也可以尝试50

作者
Author: ljh123 时间: 2024-11-30 13:56

sobereva 发表于 2024-7-5 08:38
还行。想更快平衡在350附近也可以尝试50

sob老师，请问这个time constant可以设置很小跑产生相吗，因为我设置100仍然偏离设定温度40K，但是设置为20可以稳定在设定温度

作者
Author: sobereva 时间: 2024-12-2 04:59

ljh123 发表于 2024-11-30 13:56
sob老师，请问这个time constant可以设置很小跑产生相吗，因为我设置100仍然偏离设定温度40K，但是设置为 ...

可以。但建议先尝试50

作者
Author: wilely 时间: 2025-6-29 16:23
你这个是按照孤立体系来算的把？

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