计算化学公社

 找回密码 Forget password
 注册 Register

用ztop构建大分子/聚合物的拓扑文件

查看数: 31892 | 评论数: 46 | 收藏 Add to favorites 61
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-7-14 21:43

正文摘要:

本帖最后由 ggdh 于 2021-8-23 17:28 编辑 这个帖子介绍如何用ztop构建大分子的拓扑文件。ztop的安装以及基本用法参见:构建及优化[小大]分子力场参数脚本ztop.py 用ztop拟合二面角可以参考 使用ztop拟合二面角 ...

回复 Reply

boqiang 发表于 Post on 2025-6-26 16:13:49
boqiang 发表于 2025-6-26 16:01
不太理解为什么要切C-F键,尝试了一下,得不到想要的结构,反而是A1H与A1T有连接,我的目的是两端的C-H按 ...

按照这个方法,先定义头和尾,中间定义两端C-H,可以实现,感谢。
boqiang 发表于 Post on 2025-6-26 16:01:15
ggdh 发表于 2025-6-13 21:17
因为A1 只能和A1 连, A2 只能和A2 连,所以会连成你这样
你这个利益,应该用首先定义头位基团(H和T), ...

不太理解为什么要切C-F键,尝试了一下,得不到想要的结构,反而是A1H与A1T有连接,我的目的是两端的C-H按顺序连接,跟C-F没有关系, 即A1-A2-A1-A2-A1-A2...这样的顺序,而不是A1-A2-A2-A1-A1-A2....
萌兔子酱 发表于 Post on 2025-6-22 11:13:11

老师,您好,我想构造如图所示的支化结构,其中数字表示个数,A,B分别为两种基团。我在构造的过程中,总会报错。这个是我的输入指令,“ztop.py --loadlib -b S[KIGIGIGIGd]3[KIGIGIGJMS[KIGIGIGIGd]2[KIGIGIGIQNS[KIGIGIGIGd]3]]1 -o Sp3.top,Sp3.pdb”,因为支化中结构含有支化结构,这就意味着括号包含括号,是否指令不能这么进行使用或者不能进行这么构造?
ggdh 发表于 Post on 2025-6-13 21:17:14
boqiang 发表于 2025-6-9 16:27
关于接口随机匹配问题,是否可以指定连接规则。
1. 比如,结构: CH3CHF2(如图1),接口分别为A1,A2,
2. ...

因为A1 只能和A1 连, A2 只能和A2 连,所以会连成你这样
你这个利益,应该用首先定义头位基团(H和T),
头基团F那边切断(切C-F键),位点定义位A1T
中间基团,位点定义位A1H (切C-H键),A1T(切C-F)
尾部定义一个基团A1H (切C-H键),
这种情况下A1T只能和A1H连,(规则逻辑是避免位点最后一位字母重复,比如H和H不能连)
最后链接的操作是-b H[P]8T
注意去gitee更新最新版本的ztop
boqiang 发表于 Post on 2025-6-9 16:27:31
关于接口随机匹配问题,是否可以指定连接规则。
1. 比如,结构: CH3CHF2(如图1),接口分别为A1,A2,
2. 把CH2(看成A),C2F2(看成B),构建聚合度为10的聚合物时,连接方式是A-B-B-A-A-B..方式,不是按顺序A-B-A-B-A-B..的方式进行(如图2);


处理命令如下:
ztop.py -f "P;p=AAA.top;x=AAA.gro;site=A1:3-4,A2:5-6;resname=ABC;comment=test" --savelib
ztop.py --loadlib -b [P]10

202506091623422039..png (194.12 KB, 下载次数 Times of downloads: 131)

图2

图2

202506091619089955..png (62.72 KB, 下载次数 Times of downloads: 123)

图1

图1
Hyperion 发表于 Post on 2024-12-7 23:47:21
消除净电荷 记得加--total_charge 0  


-----------发自我的手机
821533922 发表于 Post on 2024-5-16 23:28:51
期待老师能出一个含随机函数的,方便构建无规共聚
一条君 发表于 Post on 2024-3-23 17:16:37
本帖最后由 一条君 于 2025-1-15 22:54 编辑

老师,10.接口参数处理,以及电荷处理中  提到【这个电荷调整会输出在屏幕上】,我的好像看不到,是bug还是调整了呢
site-pair to connect is H1-P1, Left site: P2_1
site-pair to connect is P2_1-A21, Left site: A22_2
....
site-pair to connect is D11_40-P1, Left site: P2_41
site-pair to connect is P2_41-T2, Left site:
HP[APDP]10T build susscessful
_________________
答:看帮助ztop.py --help
有个选项--total_charge 0 可以把净电荷调为0



一条君 发表于 Post on 2024-3-23 13:56:39
本帖最后由 一条君 于 2024-3-25 15:06 编辑
neocc 发表于 2022-12-8 01:29
老师您的例子里面切出的非头尾单元都是中心或者镜面对称的结构吧。如果想做聚乙烯醇长链的话需要怎么办呢? ...

您好,分享一下我的想法,不知可行否
【4.连接位点命名和连接规则 a.第一个字符通常和片段的标签一致,不用来匹配连接位点】文中说通常一样,但是好像也可以不一样,比如【d.比如一个片段上有位点A1,A2,另外一个片段上有位点A1,B2,】,出现了两个不同字母,所以
非对称结构,切片后直接命名为A1,B1,根据【第一个字符不同的优先度更高】,程序运行时应该会形成A1B1A1B1结构吧
————————————更新
实测可行
ggdh 发表于 Post on 2023-12-14 16:03:48
821533922 发表于 2023-11-30 21:45
按照例子处理,却一直报错,不知道是哪里安装出了问题

目前重新下载应该没有我问题了,还有问题直接加我qq说
Azrael 发表于 Post on 2023-12-7 16:45:34
钟老师,我按照您的步骤安装和跑例子报错了,是不是更新的问题

截图 2023-12-07 16-44-03.png (113.9 KB, 下载次数 Times of downloads: 150)

截图 2023-12-07 16-44-03.png
821533922 发表于 Post on 2023-12-3 21:51:24
huifeideyu 发表于 2023-12-2 16:01
同样的问题,你咋解决的?

还没解决,以前用过老版本没有问题,不知道是不是更新原因
huifeideyu 发表于 Post on 2023-12-2 16:01:29
821533922 发表于 2023-11-30 21:45
按照例子处理,却一直报错,不知道是哪里安装出了问题

同样的问题,你咋解决的?
821533922 发表于 Post on 2023-11-30 21:45:15
按照例子处理,却一直报错,不知道是哪里安装出了问题

202311302144229411..jpg (72.59 KB, 下载次数 Times of downloads: 147)

202311302144229411..jpg

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

GMT+8, 2025-8-14 01:59 , Processed in 0.245168 second(s), 32 queries , Gzip On.

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