计算化学公社

 找回密码 Forget password
 注册 Register
Views: 52022|回复 Reply: 35
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

512

帖子

0

威望

4058

eV
积分
4570

Level 6 (一方通行)

truffle

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 naoki 于 2019-7-24 15:30 编辑

想请教下有没有依托GROMACS的代码或者软件,可以实现交联聚酰胺初始结构的构建。LAMMPS可以用Polymatic,但是还想问问GROMACS有没有类似的操作。没有的话我就去学习LAMMPS了



评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
luoc + 3

查看全部评分 View all ratings

No problem is insoluble in all conceivable circumstances.

236

帖子

0

威望

1229

eV
积分
1465

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

36#
发表于 Post on 2023-2-19 22:08:21 | 只看该作者 Only view this author
本帖最后由 JCenter 于 2023-2-19 22:09 编辑
naoki 发表于 2023-1-12 15:02
同行啊 不过我用的polymatic 没用过pysimm

哈哈哈,刚看到消息,我放弃pysimm,改用polymatic了。主要是polymatic有manual。另外,我刚开始做,想加下大佬的qq,有不懂的想请教大佬。我的qq是1714463007。
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

512

帖子

0

威望

4058

eV
积分
4570

Level 6 (一方通行)

truffle

35#
 楼主 Author| 发表于 Post on 2023-1-12 15:02:51 | 只看该作者 Only view this author
JCenter 发表于 2023-1-10 22:59
大佬,我最近在学习如何进行哌嗪和均苯三甲酰氯形成的高分子聚合物建模,看了公社内大量的帖子,发现需要 ...

同行啊 不过我用的polymatic 没用过pysimm
No problem is insoluble in all conceivable circumstances.

236

帖子

0

威望

1229

eV
积分
1465

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

34#
发表于 Post on 2023-1-10 22:59:46 | 只看该作者 Only view this author
naoki 发表于 2019-8-7 12:41
谢谢!polymatic我差不多知道怎么用了,就是结构单元的LAMMPS data文件还没搞定

大佬,我最近在学习如何进行哌嗪和均苯三甲酰氯形成的高分子聚合物建模,看了公社内大量的帖子,发现需要用lammps集成polymatic(由于停更,好像是用pysimm包了,搜了pysimm的使用,包括包里的示例,硬是没看懂如何操作,不知道您是否使用过pysimm进行建模,是的话,有几个问题想问下您。
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

1

帖子

0

威望

34

eV
积分
35

Level 2 能力者

33#
发表于 Post on 2022-4-29 22:42:38 | 只看该作者 Only view this author
naoki 你好,我是lammps 和 MS的初学者,最近也是在学习TMC和MPD交联聚合成聚酰胺,仔细看了你帖子,但是由于网上的资料有限,仍有许多不会操作的地方,想向你请教一下

28

帖子

0

威望

182

eV
积分
210

Level 3 能力者

32#
发表于 Post on 2021-9-23 21:43:46 | 只看该作者 Only view this author
EdwardLimit 发表于 2019-8-6 17:30
贴上一些相关的附件,我也在研究这个

这些文件怎么样才能下载呢

515

帖子

1

威望

4697

eV
积分
5232

Level 6 (一方通行)

31#
发表于 Post on 2021-5-25 20:26:26 | 只看该作者 Only view this author
Lacrimosa 发表于 2021-5-24 11:57
或许可以尝试一下Moltemplate
http://www.moltemplate.org/examples.html

谢谢,我了解下,换到lammps要学习的工具有点多了

365

帖子

5

威望

3868

eV
积分
4333

Level 6 (一方通行)

Nerv

30#
发表于 Post on 2021-5-24 11:57:49 | 只看该作者 Only view this author
tjuptz 发表于 2021-5-23 21:06
确实有时间没做研究了,最近得空又想起这个想搞明白。谢谢解惑,另外发现polymatic如果不用ms搞cvff力场 ...

或许可以尝试一下Moltemplate
http://www.moltemplate.org/examples.html
God's in his heaven,all is right with the world

515

帖子

1

威望

4697

eV
积分
5232

Level 6 (一方通行)

29#
发表于 Post on 2021-5-23 21:06:55 | 只看该作者 Only view this author
naoki 发表于 2021-5-23 14:15
好久不见ptz兄,我通过polymatic聚合时通常设置交联度为100%,实际模拟聚合时由于空间原因,总有部分TMC ...

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

512

帖子

0

威望

4058

eV
积分
4570

Level 6 (一方通行)

truffle

28#
 楼主 Author| 发表于 Post on 2021-5-23 14:15:31 | 只看该作者 Only view this author
本帖最后由 naoki 于 2021-5-23 14:17 编辑
tjuptz 发表于 2021-5-22 21:29
楼主好,最近看了下polymatic的手册,也看了几篇相关文章,对于加聚反应(均聚或共聚)还比较方便;但对于 ...

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

515

帖子

1

威望

4697

eV
积分
5232

Level 6 (一方通行)

27#
发表于 Post on 2021-5-22 21:29:57 | 只看该作者 Only view this author
楼主好,最近看了下polymatic的手册,也看了几篇相关文章,对于加聚反应(均聚或共聚)还比较方便;但对于缩聚反应我看文章里好多都是用残基来做聚合的,但是用残基的话,得到的结构最后比如像TMC就会有一个或者两个酰羰基。我看楼主在另外一个帖子里给的结构最后酰羰基部分都用羟基封端了,是如何做的呢?
ps 我看有的文献用H或者OH再做一次聚合以封端

515

帖子

1

威望

4697

eV
积分
5232

Level 6 (一方通行)

26#
发表于 Post on 2019-9-7 09:41:34 | 只看该作者 Only view this author
苏玖染 发表于 2019-9-7 09:17
我的意思是,直接把端基A交联的那条链看作支链,也就是说,端基A看成是支链的头,那么端基A的rtp肯定要处 ...

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

124

帖子

0

威望

2904

eV
积分
3028

Level 5 (御坂)

25#
发表于 Post on 2019-9-7 09:17:02 | 只看该作者 Only view this author
tjuptz 发表于 2019-9-6 11:27
您做的是如18#所述的分支的咯?
恕我愚钝,看做分支的话就不用specbond.dat而是直接在rtp定义里搞?还没 ...

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

515

帖子

1

威望

4697

eV
积分
5232

Level 6 (一方通行)

24#
发表于 Post on 2019-9-6 11:27:19 | 只看该作者 Only view this author
苏玖染 发表于 2019-9-6 10:56
这个需要试,我没有做过交联的部分,端基A如果要交联,rtp肯定不一样;
其实你可以把端基连接的那个链看 ...

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

124

帖子

0

威望

2904

eV
积分
3028

Level 5 (御坂)

23#
发表于 Post on 2019-9-6 10:56:52 | 只看该作者 Only view this author
tjuptz 发表于 2019-9-5 16:44
您好,假如交联发生在~~~ABABAB~~~的中间B和另外一个链ABABAB~~~的端基A这种写rtp有没有注意的呢?

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

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 21:28 , Processed in 0.807460 second(s), 29 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list