计算化学公社

标题: cp2k使用md模拟聚乙烯热裂解的温度问题 [打印本页]

作者
Author: xcgotk 时间: 2026-1-8 16:37
标题: cp2k使用md模拟聚乙烯热裂解的温度问题
请教各位老师，我在做聚乙烯的热裂解，实际裂解温度是900度，使用cp2k模拟md，在900度下模拟10ps，分子没有显著的断裂趋势，升高温度至3000度，则很快出现断裂，这样的模拟是否有意义，如果想尽量贴合实际实验条件，需要如何调整
第二个问题，聚乙烯断裂的话，键长到多少可以视为断裂呢

作者
Author: UW_0728. 时间: 2026-1-8 22:07
温度单位是什么说清楚，摄氏度还是K？
一般来说，动力学中设置更高温度有助于更快观察到期望的现象；想做高温动力学又想尽量贴合实验条件的话，个人认为优化好合适的晶格参数、然后高温动力学时用NVT系综即可
至于键断裂的判据，我认为按VMD轨迹显示的来看应该就可以（可以看博文http://sobereva.com/545，里面有如何在播放的时候实时自动更新连接关系的说明）

作者
Author: sobereva 时间: 2026-1-9 00:11
研究这种问题不可能做到完全贴合不是很高温的实验条件，否则模拟需要的时间尺度根本算不动。用比实际更高的温度加快现象的出现是普遍做法。

作者
Author: xcgotk 时间: 2026-1-9 12:46

UW_0728. 发表于 2026-1-8 22:07
温度单位是什么说清楚，摄氏度还是K？
一般来说，动力学中设置更高温度有助于更快观察到期望的现象；想做 ...

好的，谢谢

作者
Author: xcgotk 时间: 2026-1-9 12:46

sobereva 发表于 2026-1-9 00:11
研究这种问题不可能做到完全贴合不是很高温的实验条件，否则模拟需要的时间尺度根本算不动。用比实际更高的 ...

谢谢老师解惑

作者
Author: xcgotk 时间: 2026-1-11 16:33

UW_0728. 发表于 2026-1-8 22:07
温度单位是什么说清楚，摄氏度还是K？
一般来说，动力学中设置更高温度有助于更快观察到期望的现象；想做 ...

想请教您，关于优化合适的晶格参数如何操作，微塑料我觉得是一个孤立体系，这个晶格参数挺随意的

作者
Author: xcgotk 时间: 2026-1-11 16:37

sobereva 发表于 2026-1-9 00:11
研究这种问题不可能做到完全贴合不是很高温的实验条件，否则模拟需要的时间尺度根本算不动。用比实际更高的 ...

请问老师，这种普遍做法是否有相关的合理性支撑呢，例如做聚苯乙烯，过高的温度如果使得苯环过早地开环，这相对于实验过程就会是不合理断键，或者说在我描述的后者情况下，苯环可能开环但是一定会比碳链断裂晚，这样也是合理的

作者
Author: sobereva 时间: 2026-1-12 03:54

xcgotk 发表于 2026-1-11 16:37
请问老师，这种普遍做法是否有相关的合理性支撑呢，例如做聚苯乙烯，过高的温度如果使得苯环过早地开环， ...

如果你对于断键先后顺序都有要求的话，应该算不同位置断键的反应速率常数进行对比

高温下跑动力学，跑的次数足够多的话，对轨迹做统计分析，难断键的键比起易断键的键，断键现象平均出现的时间长度也应该更晚

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