请教无定形交联聚合物如何建模

tjuptz · 发表于 Post on 2019-8-31 21:42:53

苏玖染发表于 2019-8-31 21:00
分叉的可以搞，可以用specbond.dat形成特殊键，我也是才试成功的

谢谢提示，我尝试一下^_^

sobereva · 发表于 Post on 2019-8-31 21:53:08

苏玖染发表于 2019-8-31 21:00
分叉的可以搞，可以用specbond.dat形成特殊键，我也是才试成功的

有例子没？就像http://bbs.keinsci.com/thread-14423-1-1.html里文章中的体系

苏玖染 · 发表于 Post on 2019-9-1 16:20:17

本帖最后由苏玖染于 2019-9-1 16:23 编辑

没有教程额，我只是试出来的。

其实使用Amber是可以比较容易的做支化型聚合物的，因为可以使用bond连接功能，所以制作amber的非标准残基，去处理就行了。但是这只能做gaff力场的，如果要做opls或者其它力场的，那么可以将amber这一步作为生成构型的一步。而拓扑文件，如果要使用gmx跑模拟的话，则根据自己所需力场去做gmx的非标准残基。

我打个比方哈，比如我们要生成ABABABAB交替共聚物，在B上面又连接CCCCC，我那么按照正常的流程先把A、B、C做成gmx非标准残基，由于我们有了构型，直接用pdb2gmx就行了，但是需要注意两点，1. 就是做B的rtp文件的时候不需要表明它与C连接；2.做C的rtp文件的时候需要做两个，第一个是连接B的rtp（不需要标明与B成键），第二个是后面和自身C相连的rtp（需要标明与自身形成键）。后面就可以使用specbond.dat文件中添加键的方式连接B和C

specbond.dat存在于$GMXLIB/share/gromacs/top文件夹下，我们要添加以下一行，并把第一行的条目数改成9（因为原来有8个条目）。

B B的连接原子 1 C C的连接原子 1 连接原子之间的键长 B C'

这里的C是与自己相连的残基名；C‘是与B连接的残基名

后面使用下面命令就行了
$ gmx pdb2gmx -f DDA.pdb -merge all

pdb2gmx会加载oplsaa.ff文件夹下的rtp文件以及其上层目录下的specbond.dat文件，也会自动识别pdb文件中的残基，-merge all是将所有残合并成一个分子，因为我们在创建pdb的过程中有支链存在，所以pdb文件中会有TER这样的标识。另外，由于pdb2gmx不能识别异常二面角，所以在生成的拓扑里面一定要将异常二面角的函数类型改成4。

苏玖染 · 发表于 Post on 2019-9-1 16:32:13

由于高分子一般体积都很大，如果做模拟的话，还需要粗粒化，如果使用Martini的话，不知道可不可以把非标准残基做成Martni形式的。如果可以的话，那么做粗粒化也没问题。但是Martini的参数化有点麻烦，需要多尝试。

naoki · 发表于 Post on 2019-9-5 16:31:17

tjuptz 发表于 2019-8-31 09:16
老师，这种具有分支结构的聚合物pdb2gmx就不适用了吧

想请教一下您这种聚合物的拓扑用什么工具获取比较好呢

tjuptz · 发表于 Post on 2019-9-5 16:39:01

naoki 发表于 2019-9-5 16:31
想请教一下您这种聚合物的拓扑用什么工具获取比较好呢

本来想只能用x2top了，不过我觉得18#的方法是可以的

tjuptz · 发表于 Post on 2019-9-5 16:44:51

苏玖染发表于 2019-9-1 16:20
没有教程额，我只是试出来的。

其实使用Amber是可以比较容易的做支化型聚合物的，因为可以使用bond连接 ...

您好，假如交联发生在~~~ABABAB~~~的中间B和另外一个链ABABAB~~~的端基A这种写rtp有没有注意的呢？

苏玖染 · 发表于 Post on 2019-9-6 10:56:52

tjuptz 发表于 2019-9-5 16:44
您好，假如交联发生在~~~ABABAB~~~的中间B和另外一个链ABABAB~~~的端基A这种写rtp有没有注意的呢？

这个需要试，我没有做过交联的部分，端基A如果要交联，rtp肯定不一样；
其实你可以把端基连接的那个链看作支链处理

tjuptz · 发表于 Post on 2019-9-6 11:27:19

苏玖染发表于 2019-9-6 10:56
这个需要试，我没有做过交联的部分，端基A如果要交联，rtp肯定不一样；
其实你可以把端基连接的那个链看 ...

您做的是如18#所述的分支的咯？
恕我愚钝，看做分支的话就不用specbond.dat而是直接在rtp定义里搞？还没想清楚

苏玖染 · 发表于 Post on 2019-9-7 09:17:02

tjuptz 发表于 2019-9-6 11:27
您做的是如18#所述的分支的咯？
恕我愚钝，看做分支的话就不用specbond.dat而是直接在rtp定义里搞？还没 ...

我的意思是，直接把端基A交联的那条链看作支链，也就是说，端基A看成是支链的头，那么端基A的rtp肯定要处理一下，然后只用特殊键连接

tjuptz · 发表于 Post on 2019-9-7 09:41:34

苏玖染发表于 2019-9-7 09:17
我的意思是，直接把端基A交联的那条链看作支链，也就是说，端基A看成是支链的头，那么端基A的rtp肯定要处 ...

好的，我再研究一下怎么处理

tjuptz · 发表于 Post on 2021-5-22 21:29:57

楼主好，最近看了下polymatic的手册，也看了几篇相关文章，对于加聚反应（均聚或共聚）还比较方便；但对于缩聚反应我看文章里好多都是用残基来做聚合的，但是用残基的话，得到的结构最后比如像TMC就会有一个或者两个酰羰基。我看楼主在另外一个帖子里给的结构最后酰羰基部分都用羟基封端了，是如何做的呢？
ps 我看有的文献用H或者OH再做一次聚合以封端

naoki · 发表于 Post on 2021-5-23 14:15:31

本帖最后由 naoki 于 2021-5-23 14:17 编辑

tjuptz 发表于 2021-5-22 21:29
楼主好，最近看了下polymatic的手册，也看了几篇相关文章，对于加聚反应（均聚或共聚）还比较方便；但对于 ...

好久不见ptz兄，我通过polymatic聚合时通常设置交联度为100%，实际模拟聚合时由于空间原因，总有部分TMC酰羰基没有连上N原子，但是数量很少，最多十来个，我是手动补羟基，以表示未交联的部分会水解变为羧基（虽然这对我所研究的问题没啥影响）。之后再跑lammps优化。使用H或OH再跑交联也可以，对我来说还是手动补几个来得更快。

tjuptz · 发表于 Post on 2021-5-23 21:06:55

naoki 发表于 2021-5-23 14:15
好久不见ptz兄，我通过polymatic聚合时通常设置交联度为100%，实际模拟聚合时由于空间原因，总有部分TMC ...

确实有时间没做研究了，最近得空又想起这个想搞明白。谢谢解惑，另外发现polymatic如果不用ms搞cvff力场的话，比如我用opls手动把完整单体改残基data，那个先列出类型再和序号对应起来很麻烦，而且成键后atom type必须改一下，gaff更没有直接能生成单体data 的。现在考虑转到ms🤔

Lacrimosa · 发表于 Post on 2021-5-24 11:57:49

tjuptz 发表于 2021-5-23 21:06
确实有时间没做研究了，最近得空又想起这个想搞明白。谢谢解惑，另外发现polymatic如果不用ms搞cvff力场 ...

或许可以尝试一下Moltemplate
http://www.moltemplate.org/examples.html

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

[GROMACS] 请教无定形交联聚合物如何建模

评分 Rate