cp2k做AIMD前预平衡或者结构优化是否为必需？

wwp · 发表于 Post on 2022-2-5 21:35:54

各位老师好，我想请教一下在进行AIMD模拟模拟前，体系的预平衡或者是预优化是否为必需？或者说是初始结构可以自行进行调整然后直接进行AIMD
模拟？？？感谢各位老师！！！！

sobereva · 发表于 Post on 2022-2-6 02:44:42

是否需要结构优化看模型的合理性
如果本身就没有严重不合理接触、过强的斥力，比如直接用较准确的晶体结构当初始结构跑，可以直接跑MD而不需要先优化
预平衡看你的目的。如果你要区分平衡相（预平衡）和产生相，在产生相部分收集数据，预平衡总是要做的

wwp · 发表于 Post on 2022-2-6 08:36:14

sobereva 发表于 2022-2-6 02:44
是否需要结构优化看模型的合理性
如果本身就没有严重不合理接触、过强的斥力，比如直接用较准确的晶体结构 ...

sob老师您好，谢谢您的回答！！！我目前在研究有机物在铂表面的解离过程，主要是想请教一下在添加有机物分子到铂表面上时，分子的摆放位置是否可以自己选择？
按您的答复，我的理解就是只要结构相对合理，就是可以自己选择摆放吗？
谢谢老师，明白了预平衡和采集数据的区别。

sobereva · 发表于 Post on 2022-2-6 12:31:11

wwp 发表于 2022-2-6 08:36
sob老师您好，谢谢您的回答！！！我目前在研究有机物在铂表面的解离过程，主要是想请教一下在添加有机物 ...

我不清楚你说的摆放位置具体是怎么个选择法
不知道你是要找过渡态还是优化吸附构型什么的

wwp · 发表于 Post on 2022-2-6 18:50:19

sobereva 发表于 2022-2-6 12:31
我不清楚你说的摆放位置具体是怎么个选择法
不知道你是要找过渡态还是优化吸附构型什么的

sob老师您好，可能是我表达没有到位，麻烦您再看一下。类似于将一个水分子人为地放到铂的表面（并没有进行整个体系的结构优化），是选择将氢原子朝向铂表面，还是将氧原子朝向铂表面。这两种初始结构都可以作为MD的初始结构吗？
只是说两种结构通过模拟之后最终得到的结果不一样吗？（我的理解）
谢谢老师您的指点。

Eudaimonia · 发表于 Post on 2022-2-6 19:38:31

wwp 发表于 2022-2-6 18:50
sob老师您好，可能是我表达没有到位，麻烦您再看一下。类似于将一个水分子人为地放到铂的表面（并没有进 ...

对DFT计算而言优化出的结果可能会有区别，但是对MD可能区别就比较小了，而且就你截图这么简单的一个体系你直接跑一下就好了啊，很快的；
一般问问题建议先自己尝试再说，不可能其他人帮你确定大多数参数；

wwp · 发表于 Post on 2022-2-6 19:59:27

Eudaimonia 发表于 2022-2-6 19:38
对DFT计算而言优化出的结果可能会有区别，但是对MD可能区别就比较小了，而且就你截图这么简单的一个体系 ...

好的，我只是用着这个体系打个比方。我想表达的就是对于md计算而言，初始结构中水的位置并没有严格的要求吗？因为看到有的文献中，选取的初始结构都是通过DFT结构优化得到的最稳定的结构。谢谢老师的指点！！！

Eudaimonia · 发表于 Post on 2022-2-6 22:12:28

wwp 发表于 2022-2-6 19:59
好的，我只是用着这个体系打个比方。我想表达的就是对于md计算而言，初始结构中水的位置并没有严格的要求 ...

我基本没有怎么接触过MD，不太好说；
但是仅我对DFT的了解而言，一旦你的有机物存在多个官能团，DFT大多只能在局部得到最优解，具体以何种方式吸附在催化剂表面可能需要大量的初始构型测试，比较优化后的吸附能才能得到相对靠谱的最稳定结构；（但是带来的工作量会非常大）

期待其他老师的回答；

wwp · 发表于 Post on 2022-2-6 22:23:56

Eudaimonia 发表于 2022-2-6 22:12
我基本没有怎么接触过MD，不太好说；
但是仅我对DFT的了解而言，一旦你的有机物存在多个官能团，DFT大多 ...

好的谢谢老师。“DFT大多只能在局部得到最优解”说得很有道理

丁越 · 发表于 Post on 2022-2-6 22:26:19

wwp 发表于 2022-2-6 19:59
好的，我只是用着这个体系打个比方。我想表达的就是对于md计算而言，初始结构中水的位置并没有严格的要求 ...

吸附分子的盒子可以由好多方式构建，如packmol搭建一个水盒子，然后将slab和水盒子对摞起来就行；也可以用ABCluster搭建，M$也可以搭建，再用build layer摞起来。
摞起来后整体仍然需要用DFT优化释放斥力，否则原子因为初速度过大而跑飞。

wwp · 发表于 Post on 2022-2-6 22:34:18

丁越发表于 2022-2-6 22:26
吸附分子的盒子可以由好多方式构建，如packmol搭建一个水盒子，然后将slab和水盒子对摞起来就行；也可以 ...

丁越老师您好，

您所说的DFT优化释放斥力，是否可以理解为初始结构保证各个分子之间相距不要太近（以免将结构跑散）即可？随后就通过MD进行模拟，体系达到平衡即可进行采样分析？

谢谢老师指点！！！

丁越 · 发表于 Post on 2022-2-6 22:39:33

wwp 发表于 2022-2-6 22:34
丁越老师您好，

您所说的DFT优化释放斥力，是否可以理解为初始结构保证各个分子之间相距不要太近（以 ...

就是做结构优化，如果不优化导致分子受力很大，结构一开始就会跑飞

wwp · 发表于 Post on 2022-2-6 22:47:26

丁越发表于 2022-2-6 22:39
就是做结构优化，如果不优化导致分子受力很大，结构一开始就会跑飞

丁越老师您好

好的老师，知道了。那我还想请问一下，如果没有做结构，MD也能正常跑起来，初始结构也没有跑散，体系最后也达到平衡了。这样是否也是可行的？

谢谢老师指点！！！

丁越 · 发表于 Post on 2022-2-6 22:50:56

wwp 发表于 2022-2-6 22:47
丁越老师您好

好的老师，知道了。那我还想请问一下，如果没有做结构，MD也能正常跑起来，初始结构也 ...

可以

wwp · 发表于 Post on 2022-2-7 08:19:34

丁越发表于 2022-2-6 22:50
可以

丁越老师您好

谢谢您的指点！！！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] cp2k做AIMD前预平衡或者结构优化是否为必需？