计算化学公社

标题: 问MOF与高分子链跑动力学时，MOF形状如何不改变的情况到平衡 [打印本页]

作者
Author: zaqxs234 时间: 2021-2-8 21:22
标题: 问MOF与高分子链跑动力学时，MOF形状如何不改变的情况到平衡
(, 下载次数 Times of downloads: 19) 如图，红色的为CIF导入的MOF，蓝色的为两条高分子链，在跑动力学NPT的时候，MOF会被挤压变形，如果把MOF通过Modify→Constraints冻结结构，跑NPT又跑不了，提示无法控压。我想研究的是MOF与高分子链平衡后的位置关系，分析RDF曲线，现在就是得不到平衡态，不知道有没有大神可以给小弟解答一下，如何才能保持MOF形状不改变还能跑到平衡态

作者
Author: EdwardLimit 时间: 2021-2-9 16:14
你这模型有问题，怎么跑都不可能跑出对的结果。Polymer和MOF的相溶性不是分子层面的问题，是聚集体与聚集体的问题，现在MD研究的程度也就研究两者的界面相互作用，做不到Polymer把MOF完全包起来。如果你想研究Polymer和MOF的相互作用，用MS做基本上是死路，MS的力场不能混用，UFF或者Dreding是MOF常用的力场，但是这些力场不适用于Polymer. 力场不对，MD的结果就没有意义。所以我还是建议你学lammps或者gromacs, MOF的力场一般都是有大组开发的，你可以依据自己的需求去查找文献下载，然后获取里面的参数，polymer的力场你一般用opls-aa,charm或者amber都达到发文章的要求了。

作者
Author: zaqxs234 时间: 2021-2-11 01:21

EdwardLimit 发表于 2021-2-9 16:14
你这模型有问题，怎么跑都不可能跑出对的结果。Polymer和MOF的相溶性不是分子层面的问题，是聚集体与聚集体 ...

首先感谢您的解答，但是我的MOF里含有Cu金属元素，UFF和Dreding都优化不了（报错：The model contains forcefield types that are not part of the selected forcefield.），因此我的MOF和Polymer都用的COMPASSⅡ，但是MOF结构变化较大，我导师说实验中MOF是不改变形貌的，我不知道能否继续做下去，我的初步想法就是直接用COMPASSⅡ力场跑出动力学（不考虑MOF结构变化），然后分析Polymer和MOF的相互作用（RDF），请问是否可行？（刚研一很多东西都不是很懂，问的比较多，再次感谢！）

作者
Author: tjuptz 时间: 2021-2-12 19:49

zaqxs234 发表于 2021-2-11 01:21
首先感谢您的解答，但是我的MOF里含有Cu金属元素，UFF和Dreding都优化不了（报错：The model contains fo ...

感觉不行，模型不太合理

作者
Author: EdwardLimit 时间: 2021-2-13 17:40

zaqxs234 发表于 2021-2-11 01:21
首先感谢您的解答，但是我的MOF里含有Cu金属元素，UFF和Dreding都优化不了（报错：The model contains fo ...

你COMPASS能优化是因为某些原子分配了错误的力场导致可以优化，MOF结构变化大说明力场不合适!

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