计算化学公社

标题: 求助 MSST冲击模拟发生时间问题 [打印本页]

作者
Author: joke122 时间: 2023-5-24 10:57
标题: 求助 MSST冲击模拟发生时间问题
本帖最后由 joke122 于 2023-5-24 10:59 编辑

各位老师好，我在使用CP2K进行MSST冲击模拟时，在0.1时间步长下，模拟了1.8ps，没有发生冲击（压力温度体积并没有很大变化），是我的输入文件问题还是只是因为体系原因导致冲击发生时间晚？

作者
Author: sobereva 时间: 2023-5-26 03:12
盒子太小，X方向扩成2倍，Y方向扩成3倍

如果还不行，减小CMASS再试

作者
Author: joke122 时间: 2023-5-26 09:23

sobereva 发表于 2023-5-26 03:12
盒子太小，X方向扩成2倍，Y方向扩成3倍

如果还不行，减小CMASS再试

好的老师我去试试，老师，我的输入文件的其他设置是没有问题的是吧

作者
Author: sobereva 时间: 2023-5-26 09:33

joke122 发表于 2023-5-26 09:23
好的老师我去试试，老师，我的输入文件的其他设置是没有问题的是吧

没明显问题

作者
Author: joke122 时间: 2023-5-31 22:45

sobereva 发表于 2023-5-26 09:33
没明显问题

老师，我还有一个问题，体系太小对于冲击有什么影响

作者
Author: sobereva 时间: 2023-6-1 03:50

joke122 发表于 2023-5-31 22:45
老师，我还有一个问题，体系太小对于冲击有什么影响

模拟不出真实的现象，相当于人为做了不同区域里原子之间的等价性约束

作者
Author: NGC626 时间: 2023-10-7 11:33

sobereva 发表于 2023-6-1 03:50
模拟不出真实的现象，相当于人为做了不同区域里原子之间的等价性约束

sob老师好，请问在CP2K讲义中，您在计算用MSST方法模拟石墨转变为金刚石这个算例大约耗费了多长时间呢？

作者
Author: sobereva 时间: 2023-10-7 21:44

NGC626 发表于 2023-10-7 11:33
sob老师好，请问在CP2K讲义中，您在计算用MSST方法模拟石墨转变为金刚石这个算例大约耗费了多长时间呢？

你看一下file文件包里相应out文件末尾的耗时统计便知

作者
Author: NGC626 时间: 2023-10-8 11:03

sobereva 发表于 2023-10-7 21:44
你看一下file文件包里相应out文件末尾的耗时统计便知

谢谢sob老师

作者
Author: ddddgc 时间: 2025-2-24 14:19
本帖最后由 ddddgc 于 2025-2-24 15:03 编辑

sobereva 发表于 2023-5-26 03:12
盒子太小，X方向扩成2倍，Y方向扩成3倍

如果还不行，减小CMASS再试

sob老师好，我在进行MSST模拟时也出现楼主类似的情况（PBE做的），然后想测试一下尺寸对结果的影响（但是换了xTB以更快的测试）
结果发现，相同尺寸，相同冲击设置（MD部分设置一致），xTB模拟的结果发生了压缩，但PBE的不行。
想向老师请教一下原因？

两个不同方法模拟的输入文件（均使用Multiwfn生成并修改）中，差别主要在DFT设置部分，
其中，PBE是按讲义操作设置的（开smearing），
xTB是将计算方法设为GFN1-xTB，其余默认（用了OT方法，即使依然用smearing，也会发生压缩，但第一次SCF较难收敛）

作者
Author: sobereva 时间: 2025-2-24 15:36

ddddgc 发表于 2025-2-24 14:19
sob老师好，我在进行MSST模拟时也出现楼主类似的情况（PBE做的），然后想测试一下尺寸对结果的影响（但是 ...

对势能面的描述的差异。两个方法本身就有很大差异，这没什么新鲜的

作者
Author: ddddgc 时间: 2025-2-24 16:29

sobereva 发表于 2025-2-24 15:36
对势能面的描述的差异。两个方法本身就有很大差异，这没什么新鲜的

明白了

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